Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,694 A→G 6.7% T92A (ACG→GCG)  ygiC → putative enzyme with ATPase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,6940AG 6.7% 314.2 / 9.6 108T92A (ACG→GCG) ygiCputative enzyme with ATPase activity
Reads supporting (aligned to +/- strand):  ref base A (49/50);  new base G (4/3);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

Only 100 of 109 total aligned reads displayed.

TCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAAC  >  NC_000913/3180545‑3180840
                                                                                                                                                     |                                                                                                                                                  
tCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGGGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGTAGTTTTGTGCGCCAGTCATGGCTGa                                                                                                                                                    <  2:470691/150‑1 (MQ=255)
   cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGc                                                                                                                                                 >  1:1115177/1‑150 (MQ=255)
      aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCAcc                                                                                                                                              <  1:752045/150‑1 (MQ=255)
         ttGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGc                                                                                                                                           <  2:771396/150‑1 (MQ=255)
           gAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCa                                                                                                                                         >  1:465260/1‑150 (MQ=255)
             aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtg                                                                                                                                       <  1:516519/150‑2 (MQ=255)
               aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtgac                                                                                                                                     <  1:477578/150‑4 (MQ=255)
                  tGGAAGAAGTCACTGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTACAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCttt                                                                                                                                  >  2:11216/1‑150 (MQ=255)
                   ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTa                                                                                                                                 <  1:860864/150‑1 (MQ=255)
                    gaagaaGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAt                                                                                                                                <  1:919903/150‑1 (MQ=255)
                       gaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcg                                                                                                                             <  1:1181769/150‑1 (MQ=255)
                        aagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgc                                                                                                                            >  1:1471002/1‑150 (MQ=255)
                         agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                                                           >  2:743065/1‑150 (MQ=255)
                           tCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTc                                                                                                                         <  1:1048810/150‑1 (MQ=255)
                            cACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCt                                                                                                                        >  1:1280553/1‑150 (MQ=255)
                            cACCGCAGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGCCGCACCAGCCATCGCTTTATTCGCGTCt                                                                                                                        >  1:528961/1‑150 (MQ=255)
                               cgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTTGCGTCTTGa                                                                                                                     <  1:968271/150‑1 (MQ=255)
                                gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGAt                                                                                                                    <  1:59553/150‑1 (MQ=255)
                                gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGAt                                                                                                                    >  2:325051/1‑150 (MQ=255)
                                 ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                   <  1:446875/150‑1 (MQ=255)
                                      cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                              <  2:1085226/150‑1 (MQ=255)
                                      cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                              >  1:784488/1‑150 (MQ=255)
                                      cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                              >  2:1187750/1‑150 (MQ=255)
                                        gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGTTGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGGCATGGCTGTCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                            <  1:1231681/150‑1 (MQ=255)
                                        gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                            <  1:1526971/150‑1 (MQ=255)
                                        gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                            >  1:53580/1‑150 (MQ=255)
                                        gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                            <  1:345949/150‑1 (MQ=255)
                                        gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                            >  2:1431001/1‑150 (MQ=255)
                                        gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATAGCCAGCGATGAGCTGATGACCAAATTACGCATTCCAAAACAAAACTGGAGTTTTGAGCGCCAGTCATGGCTGACGCACCAGCCATAGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                            >  2:835405/1‑150 (MQ=255)
                                         cACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTgg                                                                                                           <  1:230658/150‑1 (MQ=255)
                                           ccAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAATCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGa                                                                                                         <  2:1172447/150‑1 (MQ=255)
                                             aGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATg                                                                                                       <  1:870978/150‑1 (MQ=255)
                                                 gtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCAc                                                                                                   >  2:1359723/1‑150 (MQ=255)
                                                   gCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTg                                                                                                 >  2:476232/1‑150 (MQ=255)
                                                        aaaGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGaa                                                                                            <  2:1380057/150‑1 (MQ=255)
                                                         aaGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAAc                                                                                           >  2:777283/1‑150 (MQ=255)
                                                             ggtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTaa                                                                                       >  2:1175435/1‑150 (MQ=255)
                                                             ggtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTaa                                                                                       >  2:333772/1‑150 (MQ=255)
                                                                ggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCTCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACt                                                                                    <  1:1098267/150‑1 (MQ=255)
                                                                ggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACt                                                                                    >  1:1151299/1‑150 (MQ=255)
                                                                   aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                 >  1:1499024/1‑150 (MQ=255)
                                                                        tGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                              >  2:685150/1‑116 (MQ=255)
                                                                        tGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                              <  1:685150/116‑1 (MQ=255)
                                                                         gATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaaa                                                                           >  1:54307/1‑150 (MQ=255)
                                                                          aTCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAAt                                                                          >  1:1476024/1‑150 (MQ=255)
                                                                             gCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAAc                                                                       >  2:633454/1‑150 (MQ=255)
                                                                             gCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAAc                                                                       >  1:1325061/1‑150 (MQ=255)
                                                                                aGCGATGAGCTGATGACCAATTTCCGCATTCCACAACACACCTGGAGTTTTGTGCGCGAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                    <  2:815385/150‑1 (MQ=255)
                                                                                aGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                    <  2:804992/150‑1 (MQ=255)
                                                                                aGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                    <  2:450974/150‑1 (MQ=255)
                                                                                aGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGATGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                    <  1:1381415/150‑1 (MQ=255)
                                                                                 gCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTg                                                                                              >  2:865970/1‑123 (MQ=255)
                                                                                 gCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTg                                                                                              <  1:865970/123‑1 (MQ=255)
                                                                                  cGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGa                                                                  >  1:1226223/1‑150 (MQ=255)
                                                                                   gATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCAGGAAAATAACGCCGAt                                                                 >  1:1223884/1‑150 (MQ=255)
                                                                                         cTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCa                                                           >  2:243789/1‑150 (MQ=255)
                                                                                         cTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCa                                                           <  1:476232/150‑1 (MQ=255)
                                                                                         cTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCa                                                           >  1:1329865/1‑150 (MQ=255)
                                                                                           gatgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAAc                                                         >  1:688456/1‑150 (MQ=255)
                                                                                           gatgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAAc                                                         <  2:1030391/150‑1 (MQ=255)
                                                                                           gatgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAAc                                                         >  1:273275/1‑150 (MQ=255)
                                                                                            atgaCCAAAATCCGCCTTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGGCTTGATCTGGGGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACg                                                        <  2:783175/150‑1 (MQ=255)
                                                                                             tgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGt                                                       >  1:122824/1‑150 (MQ=255)
                                                                                               aCCAAATTCCGCATTCCAAAACACCCCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCa                                                     >  2:451101/1‑150 (MQ=255)
                                                                                                ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                    >  1:1114565/1‑150 (MQ=255)
                                                                                                ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                    >  1:493229/1‑150 (MQ=255)
                                                                                                ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                    <  1:317192/150‑1 (MQ=255)
                                                                                                  aaaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGTCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta                                                  >  1:766214/1‑150 (MQ=255)
                                                                                                  aaaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta                                                  <  1:605916/150‑1 (MQ=255)
                                                                                                  aaaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta                                                  >  2:1202912/1‑150 (MQ=255)
                                                                                                   aaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtat                                                 >  2:1421087/1‑150 (MQ=255)
                                                                                                    aTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtata                                                <  2:493229/150‑1 (MQ=255)
                                                                                                         gCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAgg                                           >  1:1027744/1‑150 (MQ=255)
                                                                                                          cATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggc                                          <  2:226929/150‑1 (MQ=255)
                                                                                                               cAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGt                                     <  1:333772/150‑1 (MQ=255)
                                                                                                               cAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGt                                     <  1:743065/150‑1 (MQ=255)
                                                                                                                 aaaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTc                                   >  1:1344303/1‑150 (MQ=255)
                                                                                                                      cacCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCa                              >  2:512797/1‑150 (MQ=255)
                                                                                                                        ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAAcc                                                                                          >  1:1005885/1‑88 (MQ=255)
                                                                                                                        ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAAcc                                                                                          <  2:1005885/88‑1 (MQ=255)
                                                                                                                        ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                    <  1:134325/126‑1 (MQ=255)
                                                                                                                        ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                    >  2:134325/1‑126 (MQ=255)
                                                                                                                        ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGt                            <  1:11216/150‑1 (MQ=255)
                                                                                                                         cTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTg                           <  1:1266957/150‑1 (MQ=255)
                                                                                                                          tGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTgg                          <  2:59535/150‑1 (MQ=255)
                                                                                                                           ggAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGGACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGa                         >  1:851645/1‑150 (MQ=255)
                                                                                                                             aGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCTCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATc                       <  2:1115177/150‑1 (MQ=255)
                                                                                                                               ttttGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctg                     <  2:465260/150‑1 (MQ=255)
                                                                                                                               gctggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTAGCGGCTTGATCTGGCGTGGGATGGCACTGGTGAACCTGAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctg                     >  2:264511/5‑150 (MQ=255)
                                                                                                                                 ttGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggc                   <  1:1187750/150‑1 (MQ=255)
                                                                                                                                 ttGCGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggc                   <  1:835405/150‑1 (MQ=255)
                                                                                                                                 tggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggc                   >  1:394303/3‑150 (MQ=255)
                                                                                                                                    tgCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctgg                <  1:1175435/150‑1 (MQ=255)
                                                                                                                                       gcCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAg             <  2:1344303/150‑1 (MQ=255)
                                                                                                                                         cAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGAt           >  2:1507261/1‑150 (MQ=255)
                                                                                                                                             gaTGGCTGGTGCGTCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGc       >  1:279471/2‑150 (MQ=255)
                                                                                                                                              aTGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCt      <  2:1347529/150‑1 (MQ=255)
                                                                                                                                              aTGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCt      <  1:1134362/150‑1 (MQ=255)
                                                                                                                                                ggCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTa    >  2:1492519/1‑150 (MQ=255)
                                                                                                                                                  cTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAAc  <  2:394303/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAAC  >  NC_000913/3180545‑3180840

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