Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,521 G→A 21.8% P460S (CCA→TCA)  cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,5210GA21.8% 219.3 / 51.9 110P460S (CCA→TCA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (42/43);  new base A (15/9);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 116 total aligned reads displayed.

AGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGT  >  NC_000913/2890425‑2890619
                                                                                                |                                                                                                  
aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGcc                                                                                                 >  1:664190/1‑100 (MQ=255)
 gTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                                 >  2:741342/1‑99 (MQ=255)
    cGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGCCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATc                                                                                             <  1:680445/100‑1 (MQ=255)
     gggttgccAAAGAACAGCCAGTTTTTCCCCGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCGATGCCGGTGCCTGGGCCAATCa                                                                                            <  1:190027/94‑1 (MQ=255)
      gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCAt                                                                                           <  1:2570626/100‑1 (MQ=255)
      gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCAt                                                                                           <  1:595412/100‑1 (MQ=255)
       aTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATa                                                                                          <  2:2427130/100‑1 (MQ=255)
         tACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAAt                                                                                        <  1:1044902/100‑1 (MQ=255)
           ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATaa                                                                                         >  2:2414723/1‑97 (MQ=255)
           ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCa                                                                                      >  1:1044924/1‑100 (MQ=255)
           ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCa                                                                                      >  1:1894705/1‑100 (MQ=255)
             aaaGAACAGCCAGTTTTTTTTTTTTTTTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCACTACCGGTGCCTGAGCCAATCATAATCAcc                                                                                    >  1:977313/1‑100 (MQ=255)
             aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCAcc                                                                                    >  2:1255746/1‑100 (MQ=255)
             aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCAcc                                                                                    >  1:1616126/1‑100 (MQ=255)
                gAACAGCCAGTTTTTACCTNGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCggg                                                                                 >  2:1537557/1‑100 (MQ=255)
                  aCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGt                                                                               >  2:521296/1‑100 (MQ=255)
                  aCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGt                                                                               >  2:67515/1‑100 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCGATACCGGTGCCTGTGCCAATCATAATCACCgggg                                                                                >  2:684305/1‑96 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCgggg                                                                                >  1:552223/1‑96 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTc                                                                            <  1:1575153/100‑1 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTc                                                                            >  2:2341263/1‑100 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGttt                                                                             >  1:2078552/1‑99 (MQ=255)
                       cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTg                                                                          <  1:811399/100‑1 (MQ=255)
                       cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCCTAATCACCGGGGTTTCTg                                                                          >  2:770337/1‑100 (MQ=255)
                        aGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTgg                                                                         <  1:1101453/100‑1 (MQ=255)
                         gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGa                                                                        <  1:2004683/100‑1 (MQ=255)
                           ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                      >  2:1491704/1‑100 (MQ=255)
                           ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                      <  2:1003549/100‑1 (MQ=255)
                           ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                      <  1:1491969/100‑1 (MQ=255)
                                  tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                               <  2:167769/100‑1 (MQ=255)
                                  tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                               <  2:2256349/100‑1 (MQ=255)
                                  tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                               <  1:536887/100‑1 (MQ=255)
                                  tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                               <  1:572656/100‑1 (MQ=255)
                                  tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                               >  2:860772/1‑100 (MQ=255)
                                  tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                               >  1:1366807/1‑100 (MQ=255)
                                     cgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCTAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAgg                                                            <  1:1378141/100‑1 (MQ=255)
                                     cgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAgg                                                            >  2:1675/1‑100 (MQ=255)
                                     cgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAgg                                                            >  2:501847/1‑100 (MQ=255)
                                      gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGCCTGGCAGGc                                                           >  1:2367205/1‑100 (MQ=255)
                                       cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAg                                                             >  1:1640754/1‑97 (MQ=255)
                                       cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCg                                                          >  1:1661491/1‑100 (MQ=255)
                                        ttCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCgg                                                         <  2:1021271/100‑1 (MQ=255)
                                           gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                      <  2:1636061/100‑1 (MQ=255)
                                           gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                      <  2:1981194/100‑1 (MQ=255)
                                           gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                      <  1:1025923/100‑1 (MQ=255)
                                           gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                      <  1:332856/100‑1 (MQ=255)
                                             cggcggCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGtt                                                    >  2:1523650/1‑100 (MQ=255)
                                              ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTa                                                   >  1:1021209/1‑100 (MQ=255)
                                               gcggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTAt                                                  >  2:2409012/1‑100 (MQ=255)
                                                cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATc                                                 <  1:1617483/100‑1 (MQ=255)
                                                 ggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCAGCGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCg                                                <  1:1713566/100‑1 (MQ=255)
                                                  gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                               >  1:2576431/1‑100 (MQ=255)
                                                  gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                               >  1:235654/1‑100 (MQ=255)
                                                     cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTAt                                            >  2:1359060/1‑100 (MQ=255)
                                                     cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTAt                                            <  2:1615094/100‑1 (MQ=255)
                                                     cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTAt                                            <  2:2503024/100‑1 (MQ=255)
                                                         gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                        <  2:885305/100‑1 (MQ=255)
                                                         gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                        >  2:1016870/1‑100 (MQ=255)
                                                           tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                      <  2:2098267/100‑1 (MQ=255)
                                                           tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                      >  2:1433480/1‑100 (MQ=255)
                                                            gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                     >  1:2123099/1‑100 (MQ=255)
                                                            gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                     >  1:1738278/1‑100 (MQ=255)
                                                            gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                     >  1:2559804/1‑100 (MQ=255)
                                                               tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                    >  1:1047104/1‑98 (MQ=255)
                                                                aaaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTACCGTTATGTTCGATaaa                                  >  2:129693/1‑99 (MQ=255)
                                                                aaaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACAGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAAt                                 <  1:1334985/100‑1 (MQ=255)
                                                                 aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                                >  1:1669960/1‑100 (MQ=255)
                                                                 aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                                >  2:1716477/1‑100 (MQ=255)
                                                                 aaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                                <  1:1231063/100‑1 (MQ=255)
                                                                   ggCGCGGAACGGCGCAATACCGGTTCCTGGGCCAATCATAATCACCGGGGTTTCTTGATTGGCTGGCAGGCGGAATTTATCGTTATGTTCGATAAATAcg                              <  1:756460/100‑1 (MQ=255)
                                                                     cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAtgc                             >  2:1025216/1‑99 (MQ=255)
                                                                     cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                            >  2:2512008/1‑100 (MQ=255)
                                                                     cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                            >  1:1394647/1‑100 (MQ=255)
                                                                     cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGCTATGTTCGATAAATAcgcg                            >  1:315857/1‑100 (MQ=255)
                                                                      gcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCgg                           <  1:731651/100‑1 (MQ=255)
                                                                             ggCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                    <  1:424159/100‑1 (MQ=255)
                                                                             ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                    <  2:1243085/100‑1 (MQ=255)
                                                                               cgcAATGTCGGTGCCTCGGGCAAGCACAATCACCGGGGTTTCTGGATTGGCTGGCACGCGCAAGTTATCGTTAGGTTCGATAAATACGCGGACTTCGcc                   <  1:2435640/99‑1 (MQ=255)
                                                                               cgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                  <  1:769875/100‑1 (MQ=255)
                                                                               cgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                  <  2:2045079/100‑1 (MQ=255)
                                                                               cgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                  >  1:1727849/1‑100 (MQ=255)
                                                                                gcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTGCTGGATTGGCTGGCAGGTGGAAGTTATAGTTATGGTCGATAAATACGCGGACTTCGCCct                 >  1:1818940/1‑100 (MQ=255)
                                                                                 cAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctc                <  2:1913147/100‑1 (MQ=255)
                                                                                 cAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctc                <  1:1537548/100‑1 (MQ=255)
                                                                                  aaTACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctct               <  1:807137/100‑1 (MQ=255)
                                                                                    tACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttc             <  1:2190962/100‑1 (MQ=255)
                                                                                      ccGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGAGGACTTCGCCCTcttctt           >  2:2530756/1‑100 (MQ=255)
                                                                                       cGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATCAATACGCGGACTTCGCCCTCttcttc          <  1:475742/100‑1 (MQ=255)
                                                                                       cGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCttcttc          <  2:1130467/100‑1 (MQ=255)
                                                                                        ggTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc         <  2:1770083/100‑1 (MQ=255)
                                                                                        ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc         >  1:1630178/1‑100 (MQ=255)
                                                                                        ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc         <  1:1456392/100‑1 (MQ=255)
                                                                                         gTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCa        <  1:2220066/100‑1 (MQ=255)
                                                                                            ccTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgc     >  2:2246466/1‑100 (MQ=255)
                                                                                             cTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgcg    >  1:2299999/1‑100 (MQ=255)
                                                                                             cTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTAATATTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgc     >  1:1814327/1‑99 (MQ=255)
                                                                                               gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGt  <  2:2262880/100‑1 (MQ=255)
                                                                                               gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGt  >  2:2030679/1‑100 (MQ=255)
                                                                                               gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGt  >  1:1717002/1‑100 (MQ=255)
                                                                                               gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGt  <  1:280495/100‑1 (MQ=255)
                                                                                                |                                                                                                  
AGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGT  >  NC_000913/2890425‑2890619

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