Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,752,112 Δ1 bp 5.8% coding (121/975 nt) yiaU → putative LysR‑type DNA‑binding transcriptional regulator YiaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,752,1120C. 5.8% 346.2 / 13.6 103coding (121/975 nt)yiaUputative LysR‑type DNA‑binding transcriptional regulator YiaU
Reads supporting (aligned to +/- strand):  ref base C (43/54);  new base . (3/3);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 105 total aligned reads displayed.

ATAAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTATTGCAGATTAT  >  NC_000913/3751963‑3752261
                                                                                                                                                      |                                                                                                                                                     
ataAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGACATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATc                                                                                                                                                       >  2:1438296/1‑150 (MQ=255)
 taAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGCCTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCt                                                                                                                                                      >  2:707278/1‑150 (MQ=255)
 taAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCt                                                                                                                                                      <  1:201595/150‑1 (MQ=255)
    aaTTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGc                                                                                                                                                   <  2:1185651/150‑1 (MQ=255)
           tataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCGGCGACTGTACTTGGCATCGCGCAGGCCAACGTCAGCAAATATCTTGCTGAtttt                                                                                                                                            >  2:851574/1‑150 (MQ=255)
           tataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgtt                                                                                                                                            <  2:92902/150‑5 (MQ=255)
           tataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgtt                                                                                                                                            <  1:946674/150‑5 (MQ=255)
           tataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAtttt                                                                                                                                            <  1:954995/150‑1 (MQ=255)
           tataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAtttt                                                                                                                                            >  1:1129867/1‑150 (MQ=255)
            ataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                                                                                                                           <  1:47291/150‑1 (MQ=255)
            ataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                                                                                                                           >  1:53643/1‑150 (MQ=255)
             taATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGa                                                                                                                                          >  2:1455377/1‑150 (MQ=255)
             taATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGa                                                                                                                                          >  2:1573080/1‑150 (MQ=255)
             taATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCAAGTGAAAATATCAGACATGCCGCGACTGTACTTGGCATCGCAAAGGACAACGTCAGAAAATATCTTGCTGATTTTGa                                                                                                                                          >  2:267519/1‑150 (MQ=255)
               aTTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                                                                                                                                        <  2:767191/150‑1 (MQ=255)
               aTTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                                                                                                                                        <  1:1066008/150‑1 (MQ=255)
               aTTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                                                                                                                                        <  1:626481/150‑1 (MQ=255)
                      tGTTGCGATGACAAAATTACAACTTAAGTACCGGTAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGt                                                                                                                                 <  1:183306/150‑1 (MQ=255)
                          gCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGt                                                                                                                             >  2:928703/1‑150 (MQ=255)
                           cGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGtt                                                                                                                            <  1:680730/150‑1 (MQ=255)
                             aTGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGACAGTGAAAATATCAGCCATGCCTCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTTAATCAAAAGTGGGTTTa                                                                                                                          <  1:267519/150‑1 (MQ=255)
                               gACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                                                                        <  2:825614/150‑1 (MQ=255)
                               gACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                                                                        <  1:657904/150‑1 (MQ=255)
                                 cAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAg                                                                                                                      >  2:439612/1‑150 (MQ=255)
                                 cAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAg                                                                                                                      >  2:745231/1‑150 (MQ=255)
                                 cAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAg                                                                                                                      >  2:318184/1‑150 (MQ=255)
                                     aTTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                                                                                                                  <  2:1162404/150‑1 (MQ=255)
                                     aTTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                                                                                                                  >  1:1162404/1‑150 (MQ=255)
                                      ttACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCttt                                                                                                                 >  1:150826/1‑150 (MQ=255)
                                       tACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGGCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                                                                                >  2:1591278/1‑150 (MQ=255)
                                          aaCTTAAGTACCGGGAGTTAAAAATTATCTCTTTAATCGCTGGCAGTGAAAATATCAGACATGCCTCGACTGTACTTGGCATCGCACAGGCCAACGTCTGCAAATATCTTGCTTATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc                                                                                                             <  1:1003274/150‑1 (MQ=255)
                                          aaCTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc                                                                                                             >  2:331256/1‑150 (MQ=255)
                                               aaGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacc                                                                                                        <  1:1238640/150‑1 (MQ=255)
                                               aaGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacc                                                                                                        <  1:1627042/150‑1 (MQ=255)
                                               aaGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacc                                                                                                        <  1:1621588/150‑1 (MQ=255)
                                                 gTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAaccac                                                                                                      >  1:991188/1‑150 (MQ=255)
                                                      gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGc                                                                                                 >  2:829408/1‑150 (MQ=255)
                                                      gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGc                                                                                                 <  1:1434937/150‑1 (MQ=255)
                                                           ttAAAAATTATCTCGGTAATCGCTGCGAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTg                                                                                            <  2:1521420/150‑1 (MQ=255)
                                                             aaaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGAt                                                                                          >  1:510019/1‑150 (MQ=255)
                                                                 aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                                                                                      <  2:16256/150‑1 (MQ=255)
                                                                   tATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTcac                                                                                    <  1:491380/150‑1 (MQ=255)
                                                                           tAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTg                                                                            <  1:1438296/150‑1 (MQ=255)
                                                                             aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                                                                          <  1:389446/150‑1 (MQ=255)
                                                                             aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                                                                          <  2:1142846/150‑1 (MQ=255)
                                                                             aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                                                                          <  2:1606954/150‑1 (MQ=255)
                                                                             aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                                                                          <  1:395380/150‑1 (MQ=255)
                                                                                gCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCAcc                                                                       <  2:53643/150‑1 (MQ=255)
                                                                                  tGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgc                                                                     >  1:761115/1‑150 (MQ=255)
                                                                                   gCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGc                                                                                        <  2:605663/130‑1 (MQ=255)
                                                                                   gCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGc                                                                                        >  1:605663/1‑130 (MQ=255)
                                                                                   gCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgcg                                                                    >  1:800830/1‑150 (MQ=255)
                                                                                   gCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgcc                                                                    >  1:440783/1‑149 (MQ=255)
                                                                                   gCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGAGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgcg                                                                    <  1:1499529/150‑1 (MQ=255)
                                                                                      aGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTg                                                                 <  1:928703/150‑1 (MQ=255)
                                                                                       gTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGt                                                                >  1:1105511/1‑150 (MQ=255)
                                                                                         gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTa                                                              <  2:505362/150‑1 (MQ=255)
                                                                                             atatCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGt                                                          <  1:1573080/150‑1 (MQ=255)
                                                                                              tatCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGACCCACCCGGCAGCTGATGCTCACACCTTTTGGCCCCGCGCTGTTACCGTa                                                         >  1:1092117/1‑150 (MQ=255)
                                                                                              tatCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTa                                                         <  2:417199/150‑1 (MQ=255)
                                                                                               atCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTAc                                                        <  1:1241825/150‑1 (MQ=255)
                                                                                                tCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACa                                                       <  1:331256/150‑1 (MQ=255)
                                                                                                tCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACa                                                       <  2:1504793/150‑1 (MQ=255)
                                                                                                  aGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTTCATc                                                     >  2:210377/1‑150 (MQ=255)
                                                                                                  aGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATc                                                     >  1:65980/1‑150 (MQ=255)
                                                                                                  aGCCATGCCGCGACTGTACTCGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATc                                                     <  1:1622612/150‑1 (MQ=255)
                                                                                                     cATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAAc                                                  <  1:1138477/150‑1 (MQ=255)
                                                                                                      aTGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACg                                                 <  1:318184/150‑1 (MQ=255)
                                                                                                       tGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACAcc                                                                                 <  2:290679/117‑1 (MQ=255)
                                                                                                       tGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACAcc                                                                                 >  1:290679/1‑117 (MQ=255)
                                                                                                       tGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTAGGCACCGCGCTGTTACCGTACATCAACGa                                                >  1:1110906/1‑150 (MQ=255)
                                                                                                              aCTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCCCCGCGCTGTTACCGTACATCAACGACATGCTg                                         >  2:1217476/1‑150 (MQ=255)
                                                                                                                   aCTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAg                                    >  2:1501264/1‑150 (MQ=255)
                                                                                                                   aCTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAg                                    >  2:1119401/1‑150 (MQ=255)
                                                                                                                   aCTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAg                                    >  2:310151/1‑150 (MQ=255)
                                                                                                                     ttGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGaa                                  <  2:510019/150‑1 (MQ=255)
                                                                                                                         cATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGa                              <  2:501433/150‑1 (MQ=255)
                                                                                                                          aTCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAg                             >  2:788661/1‑150 (MQ=255)
                                                                                                                           tCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGc                            <  1:829408/150‑1 (MQ=255)
                                                                                                                            cGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCa                           <  1:1330290/150‑1 (MQ=255)
                                                                                                                            cGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTAGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCa                           <  2:1110906/150‑1 (MQ=255)
                                                                                                                             gCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCaa                          <  1:1055786/150‑1 (MQ=255)
                                                                                                                             gCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCaa                          <  1:1119401/150‑1 (MQ=255)
                                                                                                                             gCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCaa                          >  1:559804/1‑150 (MQ=255)
                                                                                                                              cacaGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGCCCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAAc                         >  1:1159482/1‑150 (MQ=255)
                                                                                                                              cacaGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAAc                         <  2:791728/150‑1 (MQ=255)
                                                                                                                                    ccAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTtaata                   <  1:707278/150‑1 (MQ=255)
                                                                                                                                      aaaatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAg                             <  2:858543/134‑1 (MQ=255)
                                                                                                                                      aaaatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTaataat                 >  1:1060813/5‑150 (MQ=255)
                                                                                                                                      aaaatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCAGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAg                             >  1:858543/5‑138 (MQ=255)
                                                                                                                                      aaCGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTaataat                 >  2:1517727/1‑150 (MQ=255)
                                                                                                                                        aatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAAttt               >  1:476778/3‑150 (MQ=255)
                                                                                                                                           cAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGACCCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTAt            >  2:1422446/1‑150 (MQ=255)
                                                                                                                                           cAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCCCCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGCCATGCTGGACAGAAACGAGCAACTTAATAATTTTAt            >  2:1383089/1‑150 (MQ=255)
                                                                                                                                             gCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTATTg          <  2:800830/150‑1 (MQ=255)
                                                                                                                                             gCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTATTg          >  2:833571/1‑150 (MQ=255)
                                                                                                                                              cAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGAc                                      >  1:337797/1‑121 (MQ=255)
                                                                                                                                               aaaTATCTTGCAGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTCGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTATTGCa        <  1:851574/150‑1 (MQ=255)
                                                                                                                                                  aaatCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAAACACCCGGCAGCTGATGCTCAAACCTTTTGGCACAGCGATGTTACAGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTATTGCAGa      >  1:1247546/3‑150 (MQ=255)
                                                                                                                                                      cTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTATTGCAGATtat  >  1:286819/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATAAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTACATCAACGACATGCTGGACAGAAACGAGCAACTTAATAATTTTATTGCAGATTAT  >  NC_000913/3751963‑3752261

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: