Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,207:1 +G 100% noncoding (80/87 nt) leuP ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,2071.G88.8% 287.0 / 34.5 89noncoding (80/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (5/5);  new base G (48/31);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAA  >  NC_000913/4606068‑4606337
                                                                                                                                             |                                                                                                                                     
gTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTgaga                                                                                                                                        >  1:646046/1‑140 (MQ=255)
gTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCATCGATTGGTGCGa                                                                                                                                        >  1:451641/1‑140 (MQ=255)
      aTAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGAAAGGACACTAACACCTGCAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTGGATATCGCCTCGTTTGGTGCGA‑‑ggggggg                                                                                                                             <  2:358/140‑1 (MQ=21)
      aTAATTGGGGCGCGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑ggggggg                                                                                                                             <  2:211/140‑1 (MQ=255)
        aaTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCACGTTTGGTGCGA‑‑gggggggg                                                                                                                            <  2:1316227/140‑1 (MQ=255)
        aaTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAATCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑ggggggggg                                                                                                                           >  2:942785/1‑139 (MQ=255)
          ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAGCTTACTTTTTTTGATATCGAATCGTTTGGTGCGA‑‑ggggggggga                                                                                                                          >  2:113795/1‑138 (MQ=255)
              tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTg                                                                                                                       >  1:463163/1‑140 (MQ=255)
               gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGa                                                                                                                      <  1:1255699/140‑1 (MQ=255)
               gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCATCGTTTGGTGCGAGGGGGGGGGACTTGa                                                                                                                      <  1:549739/140‑1 (MQ=255)
                  aGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAAcc                                                                                                                   >  1:1306728/1‑140 (MQ=255)
                          gACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTc                                                                                                           <  1:1155440/140‑1 (MQ=255)
                          gACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTc                                                                                                           <  1:226352/140‑1 (MQ=255)
                                aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     <  1:709233/140‑1 (MQ=255)
                                aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     <  1:709208/140‑1 (MQ=255)
                                aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAACTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     <  1:709456/140‑1 (MQ=255)
                                 aCCCCCACGTCCGTAAGGACACTAACACCTGATGCTAGCGCGTCTACCATTTCCGCCACCTTCTCACAGTCATCTTACTTTTTTTGATATCGCCTCTTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                                                                                                   <  1:84631/140‑1 (MQ=255)
                                 aCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTATACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                                                                                                   <  1:85030/140‑1 (MQ=255)
                                     ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                >  1:1382862/1‑140 (MQ=255)
                                     ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                >  1:1383275/1‑140 (MQ=255)
                                     ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTAGATAACGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGACCGTAAGAacac                                                                                                >  1:3080/1‑140 (MQ=25)
                                     ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGAAATATTACTTTTTTTGATATCGAATCGATTGGTGCGAGGGGGGGGGACTTGAACACCCACGACAGTAATAacac                                                                                                >  1:1382767/1‑140 (MQ=18)
                                      cACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGACATCTTACTTTTTTTGATATCGCCACCTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                               <  1:455112/140‑1 (MQ=37)
                                       aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                                                              >  2:13698/1‑140 (MQ=255)
                                       aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTAATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGTCTTCAACCCCCACGTCCTTCAGAACACTa                                                                                              >  2:13670/1‑140 (MQ=25)
                                        cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                             <  1:1495769/140‑1 (MQ=255)
                                        cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                             <  1:430275/140‑1 (MQ=255)
                                          tCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCGCGTGTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAacac                                                                                          >  1:341006/1‑140 (MQ=255)
                                            cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGAGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                         <  2:188981/140‑1 (MQ=255)
                                            cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCATCGTATGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTCACAcc                                                                                         >  2:1515897/1‑140 (MQ=37)
                                              tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGATCCCCCACGTCCGTAAGAACACTAACACCTg                                                                                       >  2:1356705/1‑140 (MQ=255)
                                              tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                         >  1:915520/1‑138 (MQ=255)
                                              tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                         <  2:915520/138‑1 (MQ=255)
                                              tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGACCGTAAGAACACTAAAACCTg                                                                                       >  2:493183/1‑140 (MQ=255)
                                                       tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACAAACACCTGAAGCTAgc                                                                               >  2:1464506/1‑140 (MQ=255)
                                                       tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                              >  1:1782488/1‑140 (MQ=255)
                                                       tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGACATCTTACTTTTTTTGATATCGAAGCGTTTGGTGGGA‑‑GTGGGGGGGACTTGAACACCCACTTCCGAAAGAAGAATAACACCTGATGCTAgag                                                                              >  2:1749171/1‑138 (MQ=25)
                                                       tAACACCTGAAGCTAGCGCGTCTACCAAGTCCGCCTCCTTCGCACAGGCATCTAGGTTTTTTTTCTTTCGGGGCGTTTGGTGCGAGGGGGGGGGACTTGTACCCCCACGTCCGTATGTACAGTGAGACCTGAAGGTAgcg                                                                              >  2:1209292/1‑140 (MQ=16)
                                                        aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAgcgc                                                                             <  1:1680885/140‑1 (MQ=255)
                                                           acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                                                          <  2:1689925/140‑1 (MQ=255)
                                                            ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                         >  2:1814289/1‑140 (MQ=255)
                                                            ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTGTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCAGTCTGTACACTAACACCAGAAGCTAGCGCGTCt                                                                         >  2:1923777/1‑140 (MQ=37)
                                                             cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                                                        >  1:1899683/1‑140 (MQ=255)
                                                             cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGAATCGTTTGGTGCGAGGGGGGGGGCCTTGAACACCCACGGGCGTAATAAAACTAAAACCAGAAGCTATACCGTCTa                                                                        >  1:1412735/1‑140 (MQ=25)
                                                               gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                       >  1:389631/1‑140 (MQ=255)
                                                               gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                      <  2:938619/140‑1 (MQ=255)
                                                               gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                      >  2:1666768/1‑140 (MQ=255)
                                                               gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAca                                                                      >  1:704058/1‑139 (MQ=255)
                                                               gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGACTTCTTACTTTTTTTGATATCGCATCGGTTGGTGGGA‑‑GTGGGGGGGAATTGAACAACAACTTCAGTATGATCACTCTCTCCTGAAGCTAGCGCGTCTAcc                                                                      >  2:1757071/1‑140 (MQ=19)
                                                                    tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                 >  1:1508121/1‑140 (MQ=255)
                                                                     aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                >  2:1442913/1‑140 (MQ=255)
                                                                     aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                >  1:1374191/1‑140 (MQ=255)
                                                                       ctagTCTACCTTTTCCGGCTCCTTCGCACATTGATCTTAATTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCAACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              <  2:834425/137‑1 (MQ=18)
                                                                       cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCAACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              <  2:837403/140‑1 (MQ=255)
                                                                       cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCAACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              <  2:1477577/140‑1 (MQ=255)
                                                                           tcgaCCAATTCCGCCACCTTCGCACAGGCATCTTACTATTTTTGTTATCGCCTCGATAGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                          >  2:1837784/4‑140 (MQ=21)
                                                                           tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                          <  1:144360/140‑1 (MQ=255)
                                                                           tCTACCAATTCCGCCACCTTCGCACAGTCAGCTTACATTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACTTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATACCGCCAcc                                                          >  2:1838609/1‑140 (MQ=21)
                                                                                 aaTTCCGCCACTTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                    <  1:1546441/140‑1 (MQ=255)
                                                                                 aaTTCCGCCACTTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                    <  1:1375840/140‑1 (MQ=255)
                                                                                 aaTTCCGCCACCTTCGCACAGTCAACTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                    <  1:1338469/140‑1 (MQ=255)
                                                                                     ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                <  2:1374191/140‑1 (MQ=255)
                                                                                      cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                               <  1:516108/140‑1 (MQ=255)
                                                                                        ccACCTTCGCACATTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                             <  2:801797/140‑1 (MQ=255)
                                                                                        ccACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCGGTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGGGc                                              >  1:559870/1‑140 (MQ=21)
                                                                                          aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                           <  1:536616/140‑1 (MQ=255)
                                                                                          aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                           <  1:824962/140‑1 (MQ=255)
                                                                                                      tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                               <  2:664697/76‑1 (MQ=255)
                                                                                                      tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                               >  1:664697/1‑76 (MQ=255)
                                                                                                      tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                               >  2:1601066/1‑140 (MQ=255)
                                                                                                      tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAAACCCCACGTCCGTAAGAAAACTAACACCTGAAGCTATCGCGGCAACCAATTCCGCCACCTTCGCCCAGTGCGAGCAAATTCAACgg                               >  1:1481887/1‑139 (MQ=21)
                                                                                                          cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                           >  1:445397/1‑140 (MQ=255)
                                                                                                           ttACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                         <  2:1624084/140‑1 (MQ=255)
                                                                                                              cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCACTATCAACGTGGTTTTTg                       >  2:356367/1‑140 (MQ=255)
                                                                                                              cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                       >  2:1558093/1‑140 (MQ=255)
                                                                                                              cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                       >  1:446267/1‑140 (MQ=255)
                                                                                                              cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                       >  1:446095/1‑140 (MQ=255)
                                                                                                              cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                       >  2:492519/1‑140 (MQ=255)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCACTATCAACGTGGTCTATGGt                     >  2:1425579/1‑140 (MQ=21)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                     >  2:1876129/1‑140 (MQ=255)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                     >  2:1876326/1‑140 (MQ=255)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                     >  2:1503829/1‑140 (MQ=255)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGAGAGCAATATCAAAGTGGTTTTTGGt                     >  2:236114/1‑140 (MQ=255)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCGAGCAATATCAACGTGGTTTTTGGt                     >  2:1590828/1‑140 (MQ=255)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAATCCGCCACCTTCGCCCAGTGCGAGCAATATAAACGTGGTTTTTGGt                     >  2:1085036/1‑140 (MQ=255)
                                                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAAAACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAACAACGTGGTGTTTGGt                     >  2:1401134/1‑140 (MQ=21)
                                                                                                                tattttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGAAAGAACACTAACACCTGAAGCTAGCGCGACTACCAATTCCGCCACCTTCGCCCAGAGCGAGCAAAATCAACGAGGTTTTAGGt                     >  2:1590834/3‑140 (MQ=18)
                                                                                                                    ttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACTCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                 >  2:165628/1‑140 (MQ=255)
                                                                                                                    ttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                 >  2:661577/1‑140 (MQ=255)
                                                                                                                      gATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTACCACCTGAAGCTCGCGCGTATACCAATTCCGACACATTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgg               <  1:1099222/140‑1 (MQ=21)
                                                                                                                         atCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggg            <  1:436623/140‑1 (MQ=255)
                                                                                                                         atCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTGGTGCGAgggggg           <  2:526107/140‑1 (MQ=255)
                                                                                                                          tCGCCTCGTTTGGTGCGTG‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTACCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg           <  1:609491/140‑1 (MQ=255)
                                                                                                                               tCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACtt     >  1:104231/1‑140 (MQ=255)
                                                                                                                                  tttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGa   <  1:661577/140‑1 (MQ=255)
                                                                                                                                  tgtGGTGCGA‑‑GGGGGGGGACTTCAACCCCCACGGCCGGAAGAACACTAACACCAGAAGCTAGCGCGTCTACCAATTCCGCCCCCAGCCCCCACAGCCACCAATATCAACGGGGGTGTTGGTTCGAGGGGCGGGACTAGaa  >  2:601974/3‑140 (MQ=9)
                                                                                                                                             |                                                                                                                                     
GTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAA  >  NC_000913/4606068‑4606337

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