Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,660,613 A→C 12.5% Y307D (TAC→GAC) ‡ iscS ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,660,6130AC12.5% 294.2 / 35.7 112Y307D (TAC→GAC) ‡iscScysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (45/53);  new base C (3/11);  total (48/64)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 112 total aligned reads displayed.

GAAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATG  >  NC_000913/2660467‑2660761
                                                                                                                                                  |                                                                                                                                                    
gaAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGtt                                                                                                                                                   >  1:1589330/1‑150 (MQ=255)
 aaaCGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTg                                                                                                                                                  >  2:1428218/1‑150 (MQ=255)
    cGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGTGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAg                                                                                                                                               <  1:724923/150‑1 (MQ=255)
     ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGATACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                              <  1:175459/150‑2 (MQ=255)
     ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                              <  1:618591/150‑2 (MQ=255)
     ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                              <  2:1522618/150‑2 (MQ=255)
     ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGc                                                                                                                                              >  2:907599/1‑150 (MQ=255)
     ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGc                                                                                                                                              >  1:27916/1‑150 (MQ=255)
      gATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAggt                                                                                                                                             <  1:8243/150‑3 (MQ=255)
      gATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCt                                                                                                                                             >  1:871330/1‑150 (MQ=255)
       aTAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTg                                                                                                                                            <  1:841733/150‑1 (MQ=255)
       aTAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTg                                                                                                                                            <  1:834111/150‑1 (MQ=255)
       aTAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTg                                                                                                                                            <  1:1286586/150‑1 (MQ=255)
       aTAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTg                                                                                                                                            <  2:592690/150‑1 (MQ=255)
        tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                           <  2:1610556/150‑1 (MQ=255)
         agagCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGAc                                                                                                                                          <  2:1559105/150‑1 (MQ=255)
           agCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGagtc                                                                                                                                        <  1:1298654/150‑4 (MQ=255)
            gCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGACTTGCTGACGTACAGGCGGAACCTGAAGAAACTCCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGtt                                                                                                                                       >  1:931747/1‑150 (MQ=255)
             cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCAGCTTCAACTACGTTGAAGGTGagtcgc                                                                                                                                      <  2:251258/150‑6 (MQ=255)
             cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTg                                                                                                                                      >  1:336932/1‑150 (MQ=255)
             cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTg                                                                                                                                      <  2:905463/150‑1 (MQ=255)
             cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTg                                                                                                                                      >  2:572848/1‑150 (MQ=255)
              aaTGTCCCAGCTCGTCGTTCACCCCCAGCGCGCGCAGCACGTAGGCCGGCTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGGCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTga                                                                                                                                     <  2:940567/149‑1 (MQ=255)
                 gtgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGaa                                                                                                                                  >  2:1487947/1‑150 (MQ=255)
                  tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGa                                                                                                                                           <  1:1049902/140‑1 (MQ=255)
                  tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGa                                                                                                                                           >  2:1049902/1‑140 (MQ=255)
                  tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAAt                                                                                                                                 >  1:830983/1‑150 (MQ=255)
                    ccAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGt                                                                                                                               >  2:111368/1‑150 (MQ=255)
                     cAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGtt                                                                                                                              >  2:1596036/1‑150 (MQ=255)
                     cAGCACGTCGTTCACCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTCACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGtt                                                                                                                              >  2:890835/1‑150 (MQ=255)
                      aGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTc                                                                                                                             >  2:596794/1‑150 (MQ=255)
                         tcgtcgTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGgc                                                                                                                          <  1:111368/150‑1 (MQ=255)
                         tcgtcgTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGgc                                                                                                                          >  2:907871/1‑150 (MQ=255)
                              gtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCAcc                                                                                                                     <  1:1428218/150‑1 (MQ=255)
                                   gCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCgg                                                                                                                           >  2:684253/1‑139 (MQ=255)
                                   gCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGtt                                                                                                                <  2:1449988/150‑1 (MQ=255)
                                        aGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACATACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAgg                                                                                                           <  1:720201/150‑1 (MQ=255)
                                          cgcgcgcAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTc                                                                                                         >  2:906184/1‑150 (MQ=255)
                                                cagcaCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGtt                                                                                                   >  2:1275643/1‑150 (MQ=255)
                                                  gcaCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCa                                                                                                 <  2:886536/150‑1 (MQ=255)
                                                   caCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAg                                                                                                <  1:1196906/150‑1 (MQ=255)
                                                    aCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAgg                                                                                               <  1:1464150/150‑1 (MQ=255)
                                                     cGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGt                                                                                              <  1:1218867/150‑1 (MQ=255)
                                                      gTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTa                                                                                             >  1:937863/1‑150 (MQ=255)
                                                        aGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTaaa                                                                                           <  1:1411027/150‑1 (MQ=255)
                                                        aGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTaaa                                                                                           <  1:1416596/150‑1 (MQ=255)
                                                         ggACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTGGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAc                                                                                          <  1:1488704/150‑1 (MQ=255)
                                                          gACGGTTCGAGGCTTGCTGACGTACAGGCGGAACGTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAct                                                                                         <  1:800259/150‑1 (MQ=255)
                                                          gACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAct                                                                                         <  1:1472996/150‑1 (MQ=255)
                                                             ggTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAActtct                                                                                      <  1:596794/150‑1 (MQ=255)
                                                               ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                    >  2:35814/1‑150 (MQ=255)
                                                               ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                    <  1:1596036/150‑1 (MQ=255)
                                                                tCGAGGCTTGCTGACGTACATGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCg                                                                                   <  1:1619362/150‑1 (MQ=255)
                                                                tCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCg                                                                                   <  2:1589330/150‑1 (MQ=255)
                                                                 cGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGGGAACCGTGTTCCAGTTCACCGTTCAGGTAAACTTCTTCGa                                                                                  <  2:336932/150‑1 (MQ=255)
                                                                 cGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGa                                                                                  <  1:683516/150‑1 (MQ=255)
                                                                   aGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGata                                                                                <  2:830983/150‑1 (MQ=255)
                                                                    ggCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGatat                                                                               >  1:343748/1‑150 (MQ=255)
                                                                      cTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCt                                                                             >  1:1061092/1‑150 (MQ=255)
                                                                              cGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGcc                                                                     >  1:188893/1‑150 (MQ=255)
                                                                                   aGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTcc                                                                <  1:572848/150‑1 (MQ=255)
                                                                                    ggCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCca                                                               >  2:1294485/1‑150 (MQ=255)
                                                                                    ggCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCca                                                               >  1:213453/1‑150 (MQ=255)
                                                                                     gCGGAACCTGCAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCcac                                                              <  2:1044478/150‑1 (MQ=255)
                                                                                     gCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCcac                                                              <  1:1487947/150‑1 (MQ=255)
                                                                                      cGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCcaca                                                             <  1:405069/150‑1 (MQ=255)
                                                                                           ccTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACgg                                                        >  2:856952/1‑150 (MQ=255)
                                                                                             tGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                           >  2:4580/1‑129 (MQ=255)
                                                                                             tGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGtt                                                      >  1:394765/1‑150 (MQ=255)
                                                                                               aagaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTAc                                                    <  2:188893/150‑1 (MQ=255)
                                                                                                agaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACg                                                   >  2:1179612/1‑150 (MQ=255)
                                                                                                agaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACg                                                   >  2:292277/1‑150 (MQ=255)
                                                                                                  aaaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCa                                                 <  1:906184/150‑1 (MQ=255)
                                                                                                  aaaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCa                                                 <  1:1294485/150‑1 (MQ=255)
                                                                                                   aaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAg                                                <  2:369856/150‑1 (MQ=255)
                                                                                                    aCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAgg                                               <  1:33046/150‑1 (MQ=255)
                                                                                                     cTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGc                                              >  2:593341/1‑150 (MQ=255)
                                                                                                      tGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGcc                                             >  1:510242/1‑150 (MQ=255)
                                                                                                          cggTCTTTCAGCGCCATGATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGGCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCgcgc                                         <  2:475224/149‑1 (MQ=255)
                                                                                                                   aGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTcc                                >  2:33123/1‑150 (MQ=255)
                                                                                                                    gcgcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCa                               >  2:1555659/1‑150 (MQ=255)
                                                                                                                     cgcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGGTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAt                              <  1:907871/150‑1 (MQ=255)
                                                                                                                      gcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtc                             >  2:30593/1‑150 (MQ=255)
                                                                                                                      gcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtc                             >  1:614997/1‑150 (MQ=255)
                                                                                                                       ccATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtct                            <  1:37836/150‑1 (MQ=255)
                                                                                                                       ccATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtct                            <  2:1387873/150‑1 (MQ=255)
                                                                                                                       ccATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtct                            <  1:43886/150‑1 (MQ=255)
                                                                                                                            aTCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTc                       <  2:614997/150‑1 (MQ=255)
                                                                                                                               aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGcc                    <  2:1453653/150‑1 (MQ=255)
                                                                                                                               aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGcc                    <  2:1201017/150‑1 (MQ=255)
                                                                                                                                 cGACTCACCTTCAACGTAGTTGAAGCTGCCGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCaa                  >  1:1614712/1‑149 (MQ=255)
                                                                                                                                  acgtCAGCTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtc                 >  1:172013/4‑150 (MQ=255)
                                                                                                                                   cgtCAGCTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtct                >  1:543458/3‑150 (MQ=255)
                                                                                                                                      cACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtctt             >  2:1554893/1‑150 (MQ=255)
                                                                                                                                       agcTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtcttc            >  1:140594/3‑150 (MQ=255)
                                                                                                                                        ccTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtcttct           <  1:593341/150‑1 (MQ=255)
                                                                                                                                         cTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTcttctt          >  1:260911/1‑150 (MQ=255)
                                                                                                                                            cAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtg       <  2:937863/150‑1 (MQ=255)
                                                                                                                                             aaCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc      <  2:440837/150‑1 (MQ=255)
                                                                                                                                             aaCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc      >  1:1278285/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATG  >  NC_000913/2660467‑2660761

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