Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,660,613 A→C 12.7% Y307D (TAC→GAC) ‡ 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,6130AC12.7% 232.2 / 22.1 98Y307D (TAC→GAC) ‡iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (45/40);  new base C (3/10);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 7.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

AAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGAT  >  NC_000913/2660468‑2660760
                                                                                                                                                 |                                                                                                                                                   
taaCGGATAGAGCTATGTTCGAGCTCGTCGTTCAGCCCCAGCGCGCGCTGCACGTTGGACGGTTCGAGGCTTGCTGACGTACAGGCGTAACCTGAAGAAAGTGCGAGGTCTTTCAGCTTCTTATTCAGCGACTCTCCTTCAACGTAGTTg                                                                                                                                                 <  1:637248/149‑1 (MQ=255)
aaaCGGATCGAGCTATGTGCCAGCTCGTCGTGCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGAGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTg                                                                                                                                                 <  1:99554/150‑1 (MQ=255)
aaaCGGATAGCGCTATGTGCCAGCTCGTCGTTCAGCCCCTGCGCGCGCAGCACGTTTGACGGTTCGAGGCTTGCTGACGTACAGGCGTAACCTTAAGAAAGTGCGAGGTCTTTCAGCGCGATAATCAGCGACTCTCCTTCAACGTAGTTg                                                                                                                                                 <  2:548870/150‑1 (MQ=255)
aaaCGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTg                                                                                                                                                 <  2:100092/150‑1 (MQ=255)
aaaCGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTg                                                                                                                                                 >  1:1186707/1‑150 (MQ=255)
    ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                             <  1:792846/150‑2 (MQ=255)
    ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                             <  2:1467086/150‑2 (MQ=255)
     gATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAggt                                                                                                                                            <  1:997326/150‑3 (MQ=255)
       tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                          <  1:304325/150‑1 (MQ=255)
       tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                          <  1:1527307/150‑1 (MQ=255)
       tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                          <  2:1310550/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                                                                                                                                >  1:627095/1‑150 (MQ=255)
                 tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAAt                                                                                                                                >  2:704655/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                                                                                   <  1:677111/150‑1 (MQ=255)
                                                              ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                   <  2:1186707/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                                                               <  2:1539166/150‑1 (MQ=255)
                                                                                  aGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTcc                                                               <  1:1517999/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              >  1:408633/1‑150 (MQ=255)
                                                                                                                                   cTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtctc              >  2:130193/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)
                                                                                                                                                 |                                                                                                                                                   
AAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGAT  >  NC_000913/2660468‑2660760

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: