Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,258,242 (A)8→9 intergenic (‑232/+43) cyuP ← / ← tdcG putative D/L‑serine transporter/L‑serine deaminase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,258,2341.A96.9% 226.6 / ‑0.5 65intergenic (‑224/+51)cyuP/tdcGputative D/L‑serine transporter/L‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (39/24);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

AGTTATTAGAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCAC‑‑AAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATAC  >  NC_000913/3258095‑3258366
                                                                                                                                            |                                                                                                                                     
aGTTATTAGAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATcac                                                                                                                                        <  1:546592/140‑1 (MQ=255)
     ttaGAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCAC‑‑aaaaa                                                                                                                                 <  1:314832/140‑1 (MQ=255)
      taGAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCAC‑‑aaaaaa                                                                                                                                <  1:303650/140‑1 (MQ=255)
        gAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCAC‑‑aaaaaaaa                                                                                                                              <  2:1813101/140‑1 (MQ=255)
           tAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCAC‑‑aaaaaaaaagg                                                                                                                           >  1:481771/1‑137 (MQ=255)
            aaaGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCAC‑‑aaaaaaaaaggt                                                                                                                          >  2:1528366/1‑136 (MQ=255)
              aGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGc                                                                                                                         >  1:290157/1‑140 (MQ=255)
              aGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGc                                                                                                                         >  2:899625/1‑140 (MQ=255)
                  aTGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACAt                                                                                                                     >  2:357840/1‑140 (MQ=255)
                        cgcgcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGc                                                                                                               <  2:465822/140‑1 (MQ=255)
                        cgcgcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGc                                                                                                               <  1:899625/140‑1 (MQ=255)
                          cgcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCAc                                                                                                             >  2:376376/1‑140 (MQ=255)
                           gcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCAcc                                                                                                            <  1:1304619/140‑1 (MQ=255)
                             gAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCa                                                                                                          <  2:214961/140‑1 (MQ=255)
                                gTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTtgtac                                                                                                               <  1:226837/132‑3 (MQ=255)
                                gTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTtgtac                                                                                                               >  2:226837/1‑130 (MQ=255)
                                gTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAgg                                                                                                       >  1:1495889/1‑140 (MQ=255)
                                    gagaAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGa                                                                                                   <  2:481771/140‑1 (MQ=255)
                                          aTCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTg                                                                                             >  1:185426/1‑140 (MQ=255)
                                           tCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGa                                                                                            >  1:141050/1‑140 (MQ=255)
                                                acaATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCa                                                                                       >  1:1147520/1‑140 (MQ=255)
                                                             tGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAg                                                                          >  2:72112/1‑140 (MQ=255)
                                                              gAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGa                                                                         >  2:1235376/1‑140 (MQ=255)
                                                              gAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGa                                                                         >  1:1498566/1‑140 (MQ=255)
                                                                  tAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCAc                                                                     >  1:1630226/1‑140 (MQ=255)
                                                                      ggAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTa                                                                 <  2:1498566/140‑1 (MQ=255)
                                                                      ggAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTa                                                                 <  2:185426/140‑1 (MQ=255)
                                                                      ggAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTa                                                                 <  2:1147520/140‑1 (MQ=255)
                                                                               tAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGt                                                        >  2:410611/1‑140 (MQ=255)
                                                                                  gAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGcc                                                           >  2:1778402/1‑134 (MQ=255)
                                                                                  gAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGcc                                                           <  1:1778402/134‑1 (MQ=255)
                                                                                  gAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcct                                                     <  1:468849/140‑1 (MQ=255)
                                                                                    aTTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcctcc                                                   >  1:1309322/1‑140 (MQ=255)
                                                                                    aTTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcctcc                                                   <  1:357840/140‑1 (MQ=255)
                                                                                       ttCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCcgcg                                                >  1:1800743/1‑140 (MQ=255)
                                                                                       ttCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCcgcg                                                >  1:1261335/1‑140 (MQ=255)
                                                                                       ttCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCcgcg                                                >  1:340882/1‑140 (MQ=255)
                                                                                         ctactaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTg                                              <  2:1495889/140‑1 (MQ=255)
                                                                                              aaaTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGtt                                         <  2:1261335/140‑1 (MQ=255)
                                                                                                 ttAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTcg                                      <  1:1235376/140‑1 (MQ=255)
                                                                                                  tAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTcgc                                     <  1:376376/140‑1 (MQ=255)
                                                                                                    aTCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCgg                                   >  1:1440443/1‑140 (MQ=255)
                                                                                                         aGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAttt                              >  2:1625244/1‑140 (MQ=255)
                                                                                                         aGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAttt                              >  2:70671/1‑140 (MQ=255)
                                                                                                          ggAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTa                             <  1:37838/140‑1 (MQ=255)
                                                                                                          ggAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTa                             <  1:410611/140‑1 (MQ=255)
                                                                                                           gagaATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAGGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTAt                            >  2:1047916/1‑140 (MQ=255)
                                                                                                                tatTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTcgcg                                    <  2:597589/127‑1 (MQ=255)
                                                                                                                tatTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTcgcg                                    >  1:597589/1‑127 (MQ=255)
                                                                                                                     attcTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATc                  >  2:1812858/1‑140 (MQ=255)
                                                                                                                        ctctATTTGCGGCAGATCACAAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCtt                >  2:1410075/1‑140 (MQ=255)
                                                                                                                        ctctATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCttt               >  1:358658/1‑140 (MQ=255)
                                                                                                                        ctctATTTGCGGCAGATCAC‑‑AAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTg              >  1:1005548/1‑140 (MQ=255)
                                                                                                                        ctctATTTGCGGCAGATCAC‑‑AAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTg              >  2:970837/1‑140 (MQ=255)
                                                                                                                          ctATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGc             <  2:358658/140‑1 (MQ=255)
                                                                                                                          ctATTTGCGGCAGATCACAAAAAAAAAGGTGCACAATTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGGTGTTTCGCGGTATTTATCGTTCATATCTTTGc             >  2:1661070/1‑140 (MQ=255)
                                                                                                                           tATTTGCGTCAGATCAC‑‑AAAAAAAAGGTGCACATTCGGGCACCCAAGGAAGAAACCTGAGAGAAATGTCCGTCGGAATCCACTTTAAAGGCCAGTCCTCCGCGTGACGTTTCGCGGTATTTATCGTTCACAGCTTTGcag           <  2:1185179/140‑3 (MQ=255)
                                                                                                                           tATTTGCGGCAGATCACAAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGc             <  1:676123/140‑1 (MQ=255)
                                                                                                                           tATTTGCGGCAGATCACAAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGc             >  2:676123/1‑140 (MQ=255)
                                                                                                                           tATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGcc            >  2:1382091/1‑140 (MQ=255)
                                                                                                                           tATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCATTCATATCTTTGcc            >  1:171270/1‑140 (MQ=255)
                                                                                                                            aTTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCg           <  2:1309322/140‑1 (MQ=255)
                                                                                                                            aTTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCg           >  1:80940/1‑140 (MQ=255)
                                                                                                                             tttGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCgg          >  2:1373368/1‑140 (MQ=255)
                                                                                                                               tGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGtt        >  1:284699/1‑140 (MQ=255)
                                                                                                                                gCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGttt       <  2:284699/140‑1 (MQ=255)
                                                                                                                                 cGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTc      >  2:729442/1‑140 (MQ=255)
                                                                                                                                   gCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCAt    >  2:564065/1‑140 (MQ=255)
                                                                                                                                   gCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCAt    >  2:1318899/1‑140 (MQ=255)
                                                                                                                                     aGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATAc  >  2:347142/1‑140 (MQ=255)
                                                                                                                                       aTCACAAAAAAAAAGCTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCg           <  2:487839/129‑1 (MQ=255)
                                                                                                                                       aTCACAAAAAAAAAGCTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCg           >  1:487839/1‑129 (MQ=255)
                                                                                                                                            |                                                                                                                                     
AGTTATTAGAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCAC‑‑AAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATAC  >  NC_000913/3258095‑3258366

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