Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,325,054 G→A W103* (TGG→TGA)  yciQ → DUF2207 domain‑containing protein YciQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,325,0540GA97.5% 109.9 / NA 40W103* (TGG→TGA) yciQDUF2207 domain‑containing protein YciQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/21);  minor base T (0/1);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTTAGCGAAATCGATCTTTTT  >  NC_000913/1324807‑1325300
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
ggtATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTTTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACAt                                                                                                                                                                                                                                                     <  1:43925/249‑1 (MQ=255)
            acacGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTGGGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGTCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGGTGGCGTCACCCGTGATGGGATTCCTGAACCCTGACATCTTGACTGGCat                                                                                                                                                                                                                                         <  2:320696/251‑1 (MQ=255)
                ggAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAACATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGCTGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCg                                                                                                                                                                                                                                     >  1:402897/1‑251 (MQ=255)
                      cccTTCTTATGAAGAAGGGATCTCGCTCTTTGAGGTGGAAGCCACGCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTa                                                                                                                                                                                                                                <  1:333828/250‑1 (MQ=255)
                                gaagaaGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGt                                                                                                                                                                                                                     >  1:193139/1‑251 (MQ=255)
                                        gATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCg                                                                                                                                                                                                             >  1:495909/1‑251 (MQ=255)
                                                  tttgacgctGAATCCATTCTGCACCCAGATGGGGTGCGCGTCATGAAAGAAAATATTTACTTTCAGGCGCGAAATCAGCTGATTAAGCCCCGCTTTTATCGTCATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGACATATCATATTGTGGGCGTCACCGGTGTTGGTATTCCTGAACACTGACATCTTGACTGTCATATCGGGTTAATGAGTATTGCCGTGGGCGATAaacaac                                                                                                                                                                                                    <  2:381876/242‑1 (MQ=255)
                                                   ttGATGTTGACGCCTCTCTGCAACCAGATTGGGTGCTCGACATCAAAGAAAAAATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCCCGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGTCGCTGCACTGCTGATCTATCATATTGTTGGCGTCACCCTTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAaacaa                                                                                                                                                                                                     <  1:291315/248‑1 (MQ=255)
                                                                    ctGCAACCAGATGGGGTGCTCGACATCAAAGAAAATACTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCTCTGCTGAACTAGCATATTGTTTGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACGGGCATATCGGGTTAATGAGTATTGTCGTGGTCGATAAACAACGTTTCTTGCCTCAAGGc                                                                                                                                                                                   <  1:58921/249‑1 (MQ=255)
                                                                         aCCAGATGGGGTGCTCGACAGCAAAGAAAATATTCATTGTCAGGAGCGAAATCAGCAGATTAAGCACGTCTTTTATCTTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATc                                                                                                                                                                            <  2:413074/251‑1 (MQ=255)
                                                                           cAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCAt                                                                                                                                                                          <  2:107171/251‑1 (MQ=255)
                                                                           cAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGAGACTATCAt                                                                                                                                                                          >  2:433195/1‑251 (MQ=255)
                                                                                ggggTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAATCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTAATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCATGAACCCTGACATCTTGACTGGCATATCGGGTTAAAGAGTATTGTCGTGGGCGATAAACAACGTTTATTGCCTCAAGGCGAATATCATTATc                                                                                                                                                                      >  2:296605/1‑250 (MQ=255)
                                                                                     gCTCGACATCAAAGAAAATATTCATTTTCAGGCGAGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAAAAACGTTTCTTGCCTCAAGGCGACTATCATTATAAAATTc                                                                                                                                                                >  1:90838/1‑251 (MQ=255)
                                                                                         gACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGGGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCGGCACTGCTGAACTATCATATTGTTGACGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGATTCTTGCCTCAAGGCGACTATCATTATCAAATTCAtt                                                                                                                                                             >  1:163642/1‑250 (MQ=255)
                                                                                             tCAAAGAACATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACACTGACATCTTGACTGGCATATCGGGTTAATGGGTATTGTCGTGGGCGATAAACAACGTTTTTTGCCTCAAGGCGACTATCATTATCAAATTCATTATc                                                                                                                                                          >  2:341530/1‑249 (MQ=255)
                                                                                              cAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAgg                                                                                                                                                       <  2:404372/251‑1 (MQ=255)
                                                                                                aaGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGtt                                                                                                                                                     >  1:401244/1‑251 (MQ=255)
                                                                                                           cATTGTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTATCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTTATGGTATTCCTGAACCCTGACATCTTGACTGGCATAGCGGGTTAATGAGTATTGTCGTGTGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCtt                                                                                                                                           <  2:464224/250‑1 (MQ=255)
                                                                                                            aTTTTCAGGCGCGAAATCAGCAGATTAATCACGGCTTTTATCGTGAGTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTc                                                                                                                                         <  1:340971/251‑1 (MQ=255)
                                                                                                                     cgcgAAATCAGCAGATTAAGCACGTCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAg                                                                                                                                <  2:101165/251‑1 (MQ=255)
                                                                                                                       cgAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAggg                                                                                                                              >  1:290036/1‑251 (MQ=255)
                                                                                                                                       aGCACGGCTTTTATCGTGATTTACCAAGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCt                                                                                                               >  1:344308/1‑250 (MQ=255)
                                                                                                                                             gCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTgg                                                                                                        <  2:452355/251‑1 (MQ=255)
                                                                                                                                              cTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTgg                                                                                                        >  1:505370/1‑250 (MQ=255)
                                                                                                                                                     cGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCAATGCTGAACTAACATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGACCGTGa                                                                                                 <  2:163642/250‑1 (MQ=255)
                                                                                                                                                                    cTATGGATGCAGCCTGATGGGGACGCTGCAATGCTGAACTATCATATTGTTAGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTAGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTACTGCGTGAGGGGGATTCTGGTCTGCTAATCTGGAACGTGACCGGTAACCACTGGcc                                                                                 >  1:205866/1‑251 (MQ=255)
                                                                                                                                                                                   gATGGGGACGCTGCACTGCTGAACTATCATTTTGTTGTCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGAACTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATaa                                                                  <  2:361612/251‑1 (MQ=255)
                                                                                                                                                                                           cGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGtt                                                          <  2:10299/251‑1 (MQ=255)
                                                                                                                                                                                           cGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAAGTTATAAGACCCGtt                                                          <  2:483855/251‑1 (MQ=255)
                                                                                                                                                                                              tGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCTGttttt                                                       >  2:150381/1‑251 (MQ=255)
                                                                                                                                                                                                            aTCATATTTTTGGCGTCACCCGTTATGGTTTTCCTGAACCCTGACCACTTGACTGGCTTATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTTCGTGAGGGGGATTCTGATCTGCTAATCTGTAACGTTACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTGCTCTCCAGTtctc                                          <  2:288202/250‑1 (MQ=255)
                                                                                                                                                                                                                tatTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAat                                      >  1:235538/1‑250 (MQ=255)
                                                                                                                                                                                                                  ttgttgGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATAtt                                   >  1:313473/1‑251 (MQ=255)
                                                                                                                                                                                                                    gttgGCGTCACCCGTGATGGTATTCATGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTAAAGACGACTATCATTATCACATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTCATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAATAACCGTTTTTCTCTCCAGTTCTCTAATATTgg                                 >  2:368756/1‑250 (MQ=255)
                                                                                                                                                                                                                            cACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATc                         <  1:482719/251‑1 (MQ=255)
                                                                                                                                                                                                                            cACCCGTGATGGGATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCGCTAATATTGCGGGTAATc                         <  2:235538/251‑1 (MQ=255)
                                                                                                                                                                                                                             aCCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATcc                        >  1:290055/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTTa                   <  2:137929/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccTGTCATCTTGACTGGCATATGGGGTTAATGAGTATTGTCGTTGGCGATAAACAACGTTTCTTGCCTCAAGGCGATTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATGGCGGGTAAGCCATTTAGCGAAATCGATCttttt  <  2:136040/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
GTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTTAGCGAAATCGATCTTTTT  >  NC_000913/1324807‑1325300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑