Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 6.2% N33D (AAT→GAT)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4370TC 6.2% 296.6 / 6.7 110N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (52/47);  new base C (4/3);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

Only 100 of 111 total aligned reads displayed.

TTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAG  >  NC_000913/4472296‑4472578
                                                                                                                                              |                                                                                                                                             
ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTg                                                                                                                                              >  1:224542/1‑143 (MQ=255)
  gCTTTCCGTTGCGCCTCTGTCGTCAGCACCAGATTAAGGTGATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAg                                                                                                                                            <  1:656798/143‑1 (MQ=255)
    tttCAGTTTCGCCGCTGTCGCCAGCACCAGCGTAAGGTCTTAGCGGGCGATTGCTGATCGGGAAATGATATGTTTCTGATATTGCGGATTAGGTATTGCTATTGAAGTGACGCATGGGCACAAAACAGCCCATCGATTGAggg                                                                                                                                          >  2:172023/1‑143 (MQ=255)
    tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                                                                          >  2:684866/1‑143 (MQ=255)
     ttCACTTTCGCCGCTGTCTCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgggg                                                                                                                                         <  1:890404/143‑1 (MQ=255)
        aGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTg                                                                                                                                      <  2:611498/143‑1 (MQ=255)
        aGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTg                                                                                                                                      <  2:39336/143‑1 (MQ=255)
         gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                     <  1:551156/143‑1 (MQ=255)
         gTTTCGCCGATGTCGCCAGGAGCAGATTAAGCTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                     <  2:325641/143‑1 (MQ=255)
          tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                    >  1:712194/1‑143 (MQ=255)
            gcgccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCtttttt                                                                                                                                  <  1:481021/142‑6 (MQ=255)
             cgccgcTGTCGCCAGCACCAGAGTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCttttttt                                                                                                                                 <  1:922153/143‑7 (MQ=255)
              gccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGGCGGGCTgggccgg                                                                                                                                >  2:138242/1‑139 (MQ=255)
                  cTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCGTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                                                                            >  2:349751/1‑143 (MQ=255)
                        ccAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAAc                                                                                                                      >  1:14003/1‑143 (MQ=255)
                         caggagcAGATGAAGGGCATCGCAAATTAAAACGATGATGGACATGATACGTTGCTGATATAACCAATTAGCAATCTTTTATCTGCTGAAAAAAAAAAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                     <  2:260290/137‑1 (MQ=255)
                         cAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                     >  2:274586/1‑143 (MQ=255)
                          aGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGTAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACgg                                                                                                                    <  1:889121/143‑1 (MQ=255)
                          aGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCTATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACgg                                                                                                                    <  1:690654/143‑1 (MQ=255)
                           gcaataGATTAAGGTGATCGCGAGTAAGGTCGTTTATGGAAATGATCTGTTTGTGATATAGCGGATCAGTCATCTTTTATCGCCTGACGCCTGGGAAATAAAAAGCCCTTCGATTGAGGGGCTGGGAATGGGTGATCAACggg                                                                                                                   <  2:427563/138‑1 (MQ=255)
                                aGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaagaa                                                                                                              <  2:762140/143‑1 (MQ=255)
                                   ttAAGGTCATCGCGACTATGTTCGTTTATGTAAATTTTATGTTTTTTATATATTGTATTAGACATCTTTTATCTCCTGACGCCGGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                           <  1:368442/143‑1 (MQ=255)
                                    tAAGGGCATAGAGAATAAGGTGGTCTATGGACATGATATGTTTCTGGTATATCGGATTAGCCATCTTTTATTTCCTGACGACTAGGAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAc                                                                                                          <  2:631790/143‑1 (MQ=255)
                                       ggTCAGCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                                                                       <  1:506284/143‑1 (MQ=255)
                                        gTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                                                                                      >  1:361505/1‑143 (MQ=255)
                                           aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAACAAACGGAAGAAAAACGGCAGGc                                                                                                   >  2:737012/1‑143 (MQ=255)
                                            tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGACTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGcc                                                                                                  >  1:300270/1‑143 (MQ=255)
                                                   tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCAAGCGTCt                                                                                           >  1:309498/1‑143 (MQ=255)
                                                     aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttt                                                                                         <  1:417065/143‑1 (MQ=255)
                                                      ggTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                                                                        >  1:1052290/1‑143 (MQ=255)
                                                       gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                       >  1:981764/1‑143 (MQ=255)
                                                       gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                       >  1:698047/1‑143 (MQ=255)
                                                       cgcggctATGGAATGGATATGTTGCTGATATAGCGGATTAGCCATGTTTTATATCCAGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                       <  2:305578/137‑1 (MQ=255)
                                                        tCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttttt                                                                                      >  1:574635/1‑143 (MQ=255)
                                                        tCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCATGGGCACAAAAAAGGCCCCCGATTGAGGGGCTTCGAATCGGTGATCCACGCGAATAAAAACGACAGGCCACCGTCctttct                                                                                      >  2:11571/1‑141 (MQ=255)
                                                           tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAg                                                                                   >  1:231642/1‑143 (MQ=255)
                                                                  aactGATAGGTGTCTGATATAGCGGAATAGCCATCTTGTAGCACCAGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGTGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                                                            <  2:88793/140‑1 (MQ=255)
                                                                  aaagtatagGTTTCTGTATTAGCGAATTAGCCATCTTTTATCACCTGACGCCTGGGCAAAAAAAAGCCCCTCGATGGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                                                            <  2:697551/138‑1 (MQ=255)
                                                                  aaaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                                                            <  1:684866/143‑1 (MQ=255)
                                                                  aaaTGATAGGTGTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                                                            <  1:418029/143‑1 (MQ=255)
                                                                      gctatGTTTCTGATATAGCGGATTAGCCCTCTTTTATCTCCTGACGCCTGGGCAAAAAGACGCCCCTCGATTGAGTGGCTGGTAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGa                                                                        <  1:558739/141‑1 (MQ=255)
                                                                        tatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTACTGACGCCTGGGCAAAAAAGAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAACACGGCAGGCCAGCGTCTTTTTTAAGACGCGGTAAGACa                                                                      >  2:44414/1‑143 (MQ=255)
                                                                          tGTTTCTGATATCGCGGATCAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                                                                    <  1:439045/143‑1 (MQ=255)
                                                                          tGTTGCTGAGATAGCGGATGAGCCATCTTTTATCTCGTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                                                                    <  1:55930/143‑1 (MQ=255)
                                                                          gggtgcTGATATAGGGGATTAGCCATCTTTTAGCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                                                                    <  1:79664/138‑1 (MQ=255)
                                                                           ggttCTGACATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTTGGCAAACAAAAGCCCCTCGATTGAGGGGCGGGGAATGGGTGATAAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                                                   <  1:2392/141‑1 (MQ=255)
                                                                           gTTGCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                                                   <  1:929932/143‑1 (MQ=255)
                                                                            tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGCCGCGGTAAGACAAAAt                                                                  >  1:388613/1‑143 (MQ=255)
                                                                            tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAt                                                                  >  1:940454/1‑143 (MQ=255)
                                                                               cTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                               >  1:644505/1‑143 (MQ=255)
                                                                               cTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                               >  1:589553/1‑143 (MQ=255)
                                                                               cTGAGAAAGCGGATGAGCCAGCTTTCATCCCCTGCCGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                               <  1:278241/143‑1 (MQ=255)
                                                                                 gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                                                             <  1:727208/143‑1 (MQ=255)
                                                                                 gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGACCCTCGAATGAGTGGCGGGGGATGGGTGATCCACCGCAAGAAAAACGGCAGGCCAGAGCATATTTTCAGCCGCGGTAAGACCAACTGTCGa                                                             >  2:267125/1‑143 (MQ=255)
                                                                                  atatAGCGTCATAGCGATCGGTTAGCTCCAGGCCCCGGCGCAAACAACAGCCGCTCGATTGATGGGCTGTGATTGGGTGAAAAAAAAAAAAAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                            <  2:750062/143‑1 (MQ=255)
                                                                                  agataGCGGATTAGCCATCTTTCACCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                            <  1:370158/141‑1 (MQ=255)
                                                                                    ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGCTTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACACAATGTCaaaac                                                          >  2:735377/1‑141 (MQ=255)
                                                                                    ataGCGGATTAGCCATCGTTGACCTCCTGACGCCTGGACCCAACCAAGCCCCTCGATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGCCGCGGTAAGACAAAATGTCGAaca                                                          <  1:1005935/143‑1 (MQ=255)
                                                                                     taGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac                                                         >  2:832813/1‑143 (MQ=255)
                                                                                      aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                        >  2:385146/1‑143 (MQ=255)
                                                                                      aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                        <  1:875384/143‑1 (MQ=255)
                                                                                       gcgcATCAGCCATCTCTGATCTCCTGACGCCTGGGCAAAAAAAAGCTCCTCGATTGAGGGGCTGGGAATGGGATATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac                                                       <  1:73274/143‑1 (MQ=255)
                                                                                       gCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGCTTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGGCGAacacac                                                       >  1:178968/1‑143 (MQ=255)
                                                                                         ggATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg                                                     <  2:1052290/143‑1 (MQ=255)
                                                                                                ccATCTTTTATCTCCTGTCGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                              <  1:965583/143‑1 (MQ=255)
                                                                                                  aTCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCCGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACa                                            >  2:23583/1‑143 (MQ=255)
                                                                                                   tCTGTTATGTCCTGACGCCTGGGCAAAAAAAAGCCCGGCGATTGATGGGCTGGGAATGGGAGAACAACGGCAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAt                                           <  1:145466/143‑1 (MQ=255)
                                                                                                     ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGGGAGATAGTAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACGAAATGTGGAACACACTGAACGATACAtct                                         >  1:150214/1‑142 (MQ=255)
                                                                                                     ttttATCGCCCGGCGACGGAGCAAGAAAATGCGCCTCCGTGGGGGGGCTGGGGAGGGGTGATCAACAGAAAGAAAATCGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcc                                         <  2:1047219/143‑1 (MQ=255)
                                                                                                      tttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGCGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAGGTCGAACACACTGAACCATACAtcct                                        >  2:572140/1‑143 (MQ=255)
                                                                                                        tATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctcc                                      <  1:676622/143‑1 (MQ=255)
                                                                                                          tctcCTGACGCCTGTGCAAAAAGTCGCCGCTCGATTGAGGGGCTGGGAATTGGTGATCAACGGGAAGAAGAGCGGCAGGCCAGCGTCTTTTTTCAGGCGCGGTAAGACAAAATGTCGAACACAGTGAACCATACATCCTCCCg                                    <  2:781604/143‑1 (MQ=255)
                                                                                                           ctcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                   >  1:215711/1‑143 (MQ=255)
                                                                                                            tcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGc                                  <  1:274586/143‑1 (MQ=255)
                                                                                                            tcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGc                                  >  2:348399/1‑143 (MQ=255)
                                                                                                            tcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGc                                  <  1:309119/143‑1 (MQ=255)
                                                                                                               tGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCGAACGTCTTTTTTCAGACGCGGTAAGACAAAATGCCGAACACACTGAACCATACATCCTCCCGCCaac                               >  2:1026765/1‑142 (MQ=255)
                                                                                                               ggACGGCGGGGCAAAAAAAAGCCCCTGGATTGAGGGGCGGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCaaa                               <  2:133127/142‑1 (MQ=255)
                                                                                                                gACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAAt                              >  1:559984/1‑143 (MQ=255)
                                                                                                                   gCCTGGGCAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAAGACAGGCCAGCGTCCTTTTTCAGACGCGGTAAGACAAAAAGGCGAACACACTGAACCATACCTCCTTCCAGCAAATTg                            >  2:422987/1‑143 (MQ=255)
                                                                                                                   gCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCACCGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGt                           >  1:556927/1‑143 (MQ=255)
                                                                                                                    ccTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCTGCAAATTGTc                          >  1:846477/1‑143 (MQ=255)
                                                                                                                     cTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTTGGAAACCCAGAACCATACATCCTCCCGGCAACTTGTcc                         >  2:1029232/1‑143 (MQ=255)
                                                                                                                         gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGACAAACGGCATGCCAGCGTATCAATTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCa                     >  1:290161/1‑143 (MQ=255)
                                                                                                                          gaaaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCa                     <  1:138242/141‑1 (MQ=255)
                                                                                                                            aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCAtt                   <  1:665216/143‑1 (MQ=255)
                                                                                                                            aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGAACGCCAGCTTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAAAACACTGAACCATACATCCTCCCGGCAAATTGCACGGCAtt                   >  2:698988/1‑143 (MQ=255)
                                                                                                                            aaaaaaaaGCCCCTCGATGGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCAtt                   <  2:589553/143‑1 (MQ=255)
                                                                                                                             aagaacaGCCCCTAGATTGAGGGGCTGGGTATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTGTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATATATCCTCCCGGCAAATTGTCCGGCATta                  <  2:827826/137‑1 (MQ=255)
                                                                                                                              aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATTCTCCCGGCAAATTGTCCGGCATtat                 >  1:485297/1‑143 (MQ=255)
                                                                                                                              aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATta                  <  1:187346/142‑1 (MQ=255)
                                                                                                                              aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATta                  >  2:187346/1‑142 (MQ=255)
                                                                                                                               tcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGCACAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGCCACGATTGCCTTCACACTGGCCCAGCCATCCTCCCGGCGACTTGTCCGTCATtata                >  2:752235/5‑143 (MQ=255)
                                                                                                                                  caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc             >  1:884708/2‑143 (MQ=255)
                                                                                                                                  caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAAGTGTCCGGCATGATACTc             >  1:726135/2‑143 (MQ=255)
                                                                                                                                   aGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCa            >  1:637019/1‑143 (MQ=255)
                                                                                                                                     ccccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAATGGTAGTTCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATc          <  1:269518/143‑1 (MQ=255)
                                                                                                                                     ccccGCGATAGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATc          <  1:524938/143‑1 (MQ=255)
                                                                                                                                         tCGATTGAGGGGCGGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCa      <  2:684120/143‑1 (MQ=255)
                                                                                                                                             ttGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCCTTATACTCATCGTCAGAAg  >  2:672274/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAG  >  NC_000913/4472296‑4472578

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: