Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,660,621 C→G 8.4% S304T (AGC→ACC)  iscS ← cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,660,6210CG 8.4% 294.2 / 15.1 95S304T (AGC→ACC) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base G (3/5);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

TAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCC  >  NC_000913/2660475‑2660770
                                                                                                                                                  |                                                                                                                                                     
tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                                    <  1:304325/150‑1 (MQ=255)
tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                                    <  2:1310550/150‑1 (MQ=255)
tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                                    <  1:1527307/150‑1 (MQ=255)
 agagCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGag                                                                                                                                                   <  2:684671/150‑2 (MQ=255)
  gagCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCCGCACGAAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGATGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACg                                                                                                                                                  >  1:1277075/1‑150 (MQ=255)
     cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGGCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTg                                                                                                                                               <  1:164984/150‑1 (MQ=255)
      tATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCAGCTTCAACTACGTTGAAGGTGagtcgct                                                                                                                                              <  2:196931/150‑7 (MQ=255)
       aTGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTGCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTgag                                                                                                                                             <  2:1116035/150‑1 (MQ=255)
       aTGTGCCAGCTCGTCGATCCGCCCCAGCGCGCGCAGCACGTAGGACGGTGCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTgag                                                                                                                                             <  2:848308/150‑1 (MQ=255)
        tgtgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTgaga                                                                                                                                            <  2:1300885/150‑1 (MQ=255)
          tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAAt                                                                                                                                          >  2:704655/1‑150 (MQ=255)
          tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAAt                                                                                                                                          >  1:627095/1‑150 (MQ=255)
           gCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATg                                                                                                                                         >  2:965391/1‑150 (MQ=255)
              cgCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTc                                                                                                                                      <  1:1137050/149‑1 (MQ=255)
              aGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGCAGAAACTGCGAGGTCTTTCAGCGCCATACTCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGGTGAGAAGGTTc                                                                                                                                      >  2:1542143/1‑150 (MQ=255)
              aGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTGGAGAATGTTc                                                                                                                                      >  1:1581054/1‑150 (MQ=255)
               gCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCg                                                                                                                                     <  1:1512661/150‑1 (MQ=255)
                cTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCgg                                                                                                                                    <  1:981250/150‑1 (MQ=255)
                  cgtcgtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGgcg                                                                                                                                  >  1:729842/1‑150 (MQ=255)
                      gtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCAcc                                                                                                                              >  2:508504/1‑150 (MQ=255)
                       ttCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCg                                                                                                                             <  2:540165/150‑1 (MQ=255)
                         catccccAGCGCGCGCAGCACTTAGGAGGGTGCGAGGCGTGCTGACGGAATGGCGGAACCTGAAGTTACTGCGTGGGCTTTCAGCGCCATAATCAGCGACTCGCCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCgtg                                                                                                                           <  1:833914/147‑1 (MQ=255)
                           gCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCAATAATCAGCGCCTCACCTTCACCGTAGTTGAAGCTGACGTTGAGAATGTTCGCCGCACAGTGtt                                                                                                                         >  1:554418/1‑150 (MQ=255)
                             ccaaccgcgcgcTGCACGTAGGACGGTTCGAGGCTTGGTTACGAACAGGCGGAACCTGATGAAACTGCGAGGTCTTTCAGCGCCATAATCATCGAGTCACCTTCAACGTAGTTGAAGCTGACTTTGAGAATGTTCGGCGCACCGTGTTcc                                                                                                                       <  2:1462673/145‑1 (MQ=255)
                              ccAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCa                                                                                                                      >  1:1398299/1‑150 (MQ=255)
                              ccAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCa                                                                                                                      <  1:586415/150‑1 (MQ=255)
                               cAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAg                                                                                                                     <  2:1478265/150‑1 (MQ=255)
                                cgcgcgcgcAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAgg                                                                                                                    <  1:910622/149‑1 (MQ=255)
                                aGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAgg                                                                                                                    <  2:651565/150‑1 (MQ=255)
                                           caCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGAGTTCCAGGTCACCGATCac                                                                                                         >  1:761170/1‑149 (MQ=255)
                                                 ggACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAc                                                                                                   >  1:833751/1‑150 (MQ=255)
                                                  gACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAct                                                                                                  >  2:859063/1‑150 (MQ=255)
                                                  gACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAct                                                                                                  >  2:38056/1‑150 (MQ=255)
                                                   aCGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAActt                                                                                                 >  1:1539166/1‑150 (MQ=255)
                                                       ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                             <  2:1186707/150‑1 (MQ=255)
                                                       ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                             <  1:677111/150‑1 (MQ=255)
                                                         ggAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGa                                                                                           <  1:14953/149‑1 (MQ=255)
                                                           aGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGata                                                                                         <  1:214910/150‑1 (MQ=255)
                                                             gCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCCCCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATc                                                                                       >  2:706372/1‑150 (MQ=255)
                                                              cTTGCTGCCGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCt                                                                                      <  2:136462/150‑1 (MQ=255)
                                                              cTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCt                                                                                      >  1:1413872/1‑150 (MQ=255)
                                                                 gCTGTCGTACAGGCGGAACCTGTAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCACCGTCGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTg                                                                                   <  2:849220/150‑1 (MQ=255)
                                                                 gCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCGGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTg                                                                                   <  1:675171/150‑1 (MQ=255)
                                                                 gCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTg                                                                                   <  2:627049/150‑1 (MQ=255)
                                                                   tGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGAt                                                                                 >  1:1685586/1‑150 (MQ=255)
                                                                   agACGTACACGCGGAACCTGAAGAAACTGAGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACAGTTCAGGTAAACTTCAACGATATATTTGAt                                                                                 >  2:1667833/2‑150 (MQ=255)
                                                                    gACGTACAGGCGGAACCTGAAGAAACTGCGATGTCTTTAAGCGCCATAATCAGCGACTCACTTTCAACGTAGTTGAAGCTGACGTTGTGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATg                                                                                <  2:1642075/150‑1 (MQ=255)
                                                                           aGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTcc                                                                         <  1:1517999/150‑1 (MQ=255)
                                                                           aGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTcc                                                                         <  2:1539166/150‑1 (MQ=255)
                                                                           aGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTcc                                                                         <  1:218886/150‑1 (MQ=255)
                                                                            ggCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTccc                                                                        >  1:1154448/1‑149 (MQ=255)
                                                                            ggCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCca                                                                        >  2:1216489/1‑150 (MQ=255)
                                                                              cGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCcaca                                                                      >  1:1247165/1‑150 (MQ=255)
                                                                                  aCCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACATACg                                                                  >  2:1530523/1‑150 (MQ=255)
                                                                                     tGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGtt                                                               >  1:474034/1‑150 (MQ=255)
                                                                                      gaagaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATTTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATCTCTTTGATGACGTTCCACAGACGGTTa                                                              >  2:270401/1‑150 (MQ=255)
                                                                                      gaagaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTa                                                              >  2:1303692/1‑150 (MQ=255)
                                                                                      gaagaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTAAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTa                                                              >  1:75937/1‑150 (MQ=255)
                                                                                       aagaaACTGCGAGGTCTTTCAGCGACATAATCAGCGACTCACCTTCATCGTCTTTGAAGCTGACGTTTAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTAc                                                             <  1:443240/150‑1 (MQ=255)
                                                                                           aaCTGCCAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCTACGTAGTTGAAGCTGACGTTGAGATAGGTCGGCGCACCGTGTTCCAGGGCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAg                                                         <  1:288507/150‑1 (MQ=255)
                                                                                               gCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCcg                                                     >  2:1266377/1‑150 (MQ=255)
                                                                                                         tccccgcCAAAATCAGCGACCCACCTTCCACGGAGGTGAATCTGACGGTGATAATGTTCGGCGCACCGGGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGtt                                           <  1:1428551/146‑1 (MQ=255)
                                                                                                           aGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTAGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTcc                                         <  1:100013/150‑1 (MQ=255)
                                                                                                              gcCATAATCAGCGACTCACCTTCAACGTAGTTGCAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTTATGCCTTTCCACAGACGGTTACGCAGGCAGCGCCGACGTTCCAtc                                      >  2:34041/1‑150 (MQ=255)
                                                                                                              gcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtc                                      >  2:1022213/1‑150 (MQ=255)
                                                                                                                  tAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCgg                                  >  2:1424076/1‑150 (MQ=255)
                                                                                                                  tAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCgg                                  <  1:41741/150‑1 (MQ=255)
                                                                                                                  tAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTACACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCgg                                  >  1:877145/1‑150 (MQ=255)
                                                                                                                   aaTCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCATGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGt                                 >  2:1239498/1‑150 (MQ=255)
                                                                                                                     tCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCg                               <  1:1303692/150‑1 (MQ=255)
                                                                                                                       aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGcc                             >  1:855240/1‑150 (MQ=255)
                                                                                                                       aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGcc                             >  2:720514/1‑150 (MQ=255)
                                                                                                                       aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGcc                             >  1:166722/1‑150 (MQ=255)
                                                                                                                        tcGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCa                            <  1:1259660/149‑1 (MQ=255)
                                                                                                                         cGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAt                           <  1:859063/150‑1 (MQ=255)
                                                                                                                          acgtCAGCTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtc                          >  1:1051653/4‑150 (MQ=255)
                                                                                                                            gtCAGCTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATTCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtctc                        >  2:1076730/2‑150 (MQ=255)
                                                                                                                            cTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtctc                        >  2:130193/1‑150 (MQ=255)
                                                                                                                            cTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtctc                        >  1:408633/1‑150 (MQ=255)
                                                                                                                              cagcTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGGTACGCCGGCCGCGCAGGCGTTCCCTCTCGGTCGCCATcgctg                      >  2:696010/4‑149 (MQ=255)
                                                                                                                              cACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtctt                      >  2:1470466/1‑150 (MQ=255)
                                                                                                                              cACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtctt                      >  2:28450/1‑150 (MQ=255)
                                                                                                                                 cTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTcttctt                   >  1:1364742/1‑150 (MQ=255)
                                                                                                                                 cTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTGCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCGctgctt                   >  1:548634/1‑150 (MQ=255)
                                                                                                                                   tCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCATGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCtttt                 <  1:846802/150‑1 (MQ=255)
                                                                                                                                    cAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtg                <  2:669364/150‑1 (MQ=255)
                                                                                                                                     aaCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCCGGCCGCGCAGACTTGCCATCTCGGACGCCATCTCATCTTTtgc               >  2:282597/1‑150 (MQ=255)
                                                                                                                                      aCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcg              >  1:1605427/1‑150 (MQ=255)
                                                                                                                                       cGTAGTTGAAGCTGACGTTGAGACTGTGCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATAGCTTTGATGCCGTTCCACAGCCGGTTACGCCGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcga             <  1:53216/150‑1 (MQ=255)
                                                                                                                                        gTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcgat            <  1:1667833/150‑1 (MQ=255)
                                                                                                                                           gTTGAAGCTGACGTTGAGAATGTTCGGTGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGCCGTTCCATCTCGGTCGCCATCTCTTCTTTtgcgatgcg         >  2:45742/1‑150 (MQ=255)
                                                                                                                                           gTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcgatgcg         <  1:965391/150‑1 (MQ=255)
                                                                                                                                           gTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcgatgcg         <  2:410009/150‑1 (MQ=255)
                                                                                                                                            ttGACGCTGACGGTGAGAATGTTCGGCGCACCGTGGTCCAGGTCACCGGTCAGGTAAACTTCATCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcgatgcga        <  1:270401/150‑1 (MQ=255)
                                                                                                                                               aaGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAg     <  2:919797/150‑1 (MQ=255)
                                                                                                                                                  cTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGcc  <  1:508504/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCC  >  NC_000913/2660475‑2660770

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