Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,890,521 G→A 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,5210GA98.0% 321.6 / ‑2.5 99P460S (CCA→TCA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (39/57);  total (40/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGC  >  NC_000913/2890424‑2890616
                                                                                                 |                                                                                               
aaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGc                                                                                               <  2:2895946/100‑1 (MQ=255)
 aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                              <  2:2891298/100‑1 (MQ=255)
 aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                              <  2:1386342/100‑1 (MQ=255)
 aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                              <  2:1138519/100‑1 (MQ=255)
 aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                              <  1:2246786/100‑1 (MQ=255)
  gTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                              >  2:2061106/1‑99 (MQ=255)
  gTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCa                                                                                             >  2:2200589/1‑100 (MQ=255)
   tGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCaa                                                                                            >  2:719153/1‑100 (MQ=255)
    gCGGATTACCTAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGTGCAGAACGGCGCAATACCGGTGCCTGAGCCAAt                                                                                           <  1:1901896/100‑1 (MQ=255)
      ggATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCa                                                                                         <  2:978008/100‑1 (MQ=255)
      ggATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCa                                                                                         <  2:787518/100‑1 (MQ=255)
      ggATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCa                                                                                         >  2:2309284/1‑100 (MQ=255)
       gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCAt                                                                                        <  1:1061051/100‑1 (MQ=255)
       gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCAt                                                                                        <  1:1061335/100‑1 (MQ=255)
        aTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATa                                                                                       <  1:2466180/100‑1 (MQ=255)
        aTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATa                                                                                       <  2:2083837/100‑1 (MQ=255)
          tACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAAt                                                                                     >  1:612191/1‑100 (MQ=255)
           aCCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATc                                                                                    >  2:1723032/1‑100 (MQ=255)
           aCCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATc                                                                                    >  2:1109399/1‑100 (MQ=255)
            ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCa                                                                                   >  2:2838739/1‑100 (MQ=255)
              aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                 <  2:2617107/100‑1 (MQ=255)
              aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGCCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                 >  1:529183/1‑100 (MQ=255)
               aaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCg                                                                                >  2:1802521/1‑100 (MQ=255)
               aaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCg                                                                                <  2:127829/100‑1 (MQ=255)
                aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCgg                                                                               <  2:2371904/100‑1 (MQ=255)
                aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCgg                                                                               <  2:1667048/100‑1 (MQ=255)
                    cAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGtt                                                                           <  2:1546900/100‑1 (MQ=255)
                      gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCGAATCATAATCACCGGGGATTc                                                                         >  1:338569/1‑100 (MQ=255)
                       ccAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCt                                                                        <  2:59306/100‑1 (MQ=255)
                        cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGAGGCTGCATAAAGGCGCGGAACGGCGCAATACGGGTGCGTGAGCCAATCATAATCACCGGGTTTTCTg                                                                       >  2:2756667/1‑100 (MQ=255)
                            ttttACCTGGCGCTTCGTCGGTGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                   <  2:827765/100‑1 (MQ=255)
                            ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                   >  2:1313391/1‑100 (MQ=255)
                            tgttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                   <  1:1755295/98‑1 (MQ=255)
                             tttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtg                                                                  <  2:665793/100‑1 (MQ=255)
                              ttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtgg                                                                 <  2:666635/100‑1 (MQ=255)
                               tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggc                                                                >  2:376585/1‑100 (MQ=255)
                                 ccTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctg                                                              >  1:2402150/1‑100 (MQ=255)
                                  cTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctgg                                                             >  2:2838017/1‑100 (MQ=255)
                                  cTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctgg                                                             <  2:1083117/100‑1 (MQ=255)
                                   tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                            <  2:2618807/100‑1 (MQ=255)
                                   tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                            <  1:1477172/100‑1 (MQ=255)
                                   tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                            >  2:1320538/1‑100 (MQ=255)
                                   tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                            <  1:220853/100‑1 (MQ=255)
                                   tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                            >  1:2149941/1‑100 (MQ=255)
                                   tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                            >  2:1049666/1‑100 (MQ=255)
                                    ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                           <  2:1233330/100‑1 (MQ=255)
                                    ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                           >  1:297835/1‑100 (MQ=255)
                                       gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGc                                                        <  2:173873/100‑1 (MQ=255)
                                        cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCg                                                       >  2:1217489/1‑100 (MQ=255)
                                         ttCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCgg                                                      <  1:592689/100‑1 (MQ=255)
                                          tcgtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGa                                                     <  2:2459283/100‑1 (MQ=255)
                                            gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                   <  2:2265116/100‑1 (MQ=255)
                                             tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                  <  2:2151412/100‑1 (MQ=255)
                                              cggcggCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGtt                                                 <  2:393249/100‑1 (MQ=255)
                                               ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTa                                                <  2:1541914/100‑1 (MQ=255)
                                                gcggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTAt                                               >  1:1375112/1‑100 (MQ=255)
                                                 cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATc                                              >  1:2439072/1‑100 (MQ=255)
                                                 cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATc                                              <  2:2620808/100‑1 (MQ=255)
                                                 cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATAACCGGGGTTTCTGGATGGGCTGGCAGGCGGAACTTATc                                              >  1:2413394/1‑100 (MQ=255)
                                                  ggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCg                                             <  2:1558196/100‑1 (MQ=255)
                                                  ggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCg                                             <  2:2767835/100‑1 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                            <  2:2214415/100‑1 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                            <  2:2581786/100‑1 (MQ=255)
                                                     gcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTa                                          <  2:1875769/100‑1 (MQ=255)
                                                       gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATg                                        >  2:805748/1‑100 (MQ=255)
                                                       gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATg                                        <  2:1712307/100‑1 (MQ=255)
                                                          gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                     >  1:1115847/1‑100 (MQ=255)
                                                          gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                     >  1:1203795/1‑100 (MQ=255)
                                                           ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCg                                    <  1:1072939/100‑1 (MQ=255)
                                                            tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                   <  1:719147/100‑1 (MQ=255)
                                                            tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGCTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                   >  1:2529580/1‑100 (MQ=255)
                                                             gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                  >  1:102403/1‑100 (MQ=255)
                                                             gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                  >  2:1262562/1‑100 (MQ=255)
                                                              cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                 <  2:511330/100‑1 (MQ=255)
                                                              cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                 <  1:443598/100‑1 (MQ=255)
                                                                tAAAGGCGAAGAACGGCGCAATAGCGGTAACTGCGCCGATCATAAGCACCGGCGTTTTTGGATTGGCAGGCAGGCGGAAGCTATCTTTATGTTCGATaaa                               <  1:2478330/100‑1 (MQ=255)
                                                                  aaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                             <  1:1723006/100‑1 (MQ=255)
                                                                  aaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                             <  1:2857679/100‑1 (MQ=255)
                                                                    ggCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcg                           >  2:1638896/1‑100 (MQ=255)
                                                                    ggCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcg                           <  1:2875341/100‑1 (MQ=255)
                                                                      cgcgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                         <  1:2838701/100‑1 (MQ=255)
                                                                      cgcgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                         <  2:597660/100‑1 (MQ=255)
                                                                         ggAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGAc                      >  2:1608471/1‑100 (MQ=255)
                                                                         ggAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGAc                      >  2:1155412/1‑100 (MQ=255)
                                                                             cGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCg                  <  1:740518/100‑1 (MQ=255)
                                                                              ngcgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                 >  1:2799512/2‑100 (MQ=255)
                                                                               gcgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGcc                >  1:1323091/1‑100 (MQ=255)
                                                                                cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc               >  1:1162929/1‑100 (MQ=255)
                                                                                cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc               <  2:655636/100‑1 (MQ=255)
                                                                                cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc               <  2:2149968/100‑1 (MQ=255)
                                                                                cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc               <  1:1802495/100‑1 (MQ=255)
                                                                                    aTACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtctt           <  2:212259/100‑1 (MQ=255)
                                                                                        cGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCttcttc       <  1:960470/100‑1 (MQ=255)
                                                                                        cGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCttcttc       <  2:612197/100‑1 (MQ=255)
                                                                                        cGGTGCCTGAGCCAATCATAAACACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCttcttc       >  1:1786213/1‑100 (MQ=255)
                                                                                         ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGACGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc      >  1:2629067/1‑100 (MQ=255)
                                                                                           tGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAc    >  2:2814309/1‑100 (MQ=255)
                                                                                            gCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcg   <  1:2200562/100‑1 (MQ=255)
                                                                                             ccTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgc  >  2:2352187/1‑100 (MQ=255)
                                                                                                 |                                                                                               
AAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGC  >  NC_000913/2890424‑2890616

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