Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,521 G→A 30.4% 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,5210GA30.4% 246.0 / 135.2 184P460S (CCA→TCA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (60/68);  new base A (27/29);  total (87/97)
Fisher's exact test for biased strand distribution p-value = 8.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

Only 100 of 191 total aligned reads displayed.

TAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTC  >  NC_000913/2890422‑2890620
                                                                                                   |                                                                                                   
 aaaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTggg                                                                                                    <  1:1070761/100‑1 (MQ=255)
 aaaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTggg                                                                                                    <  2:1422365/100‑1 (MQ=255)
 aaaGTGCGGATTACCAAAAAACAGCCAGTTTTTACCTGGCGTTTCGTCGGCGGCGCGTTGCTGCATCAAGGCGCGGAACTGTGCAATACCGGTGACTggg                                                                                                    <  2:1321166/100‑1 (MQ=255)
      gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCTTCGGCGGCGCTTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAAt                                                                                               <  1:1096390/100‑1 (MQ=255)
      gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAAt                                                                                               >  1:1237506/1‑100 (MQ=255)
      gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAAt                                                                                               >  2:574965/1‑100 (MQ=255)
      gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAAt                                                                                               <  1:2080253/100‑1 (MQ=255)
       cGGATTATCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCGTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATc                                                                                              <  1:3928335/100‑1 (MQ=255)
       cGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATc                                                                                              <  1:3287160/100‑1 (MQ=255)
        ggATTACCAAAAAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCa                                                                                             <  1:1806121/100‑1 (MQ=255)
         gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCAt                                                                                            <  1:3561691/100‑1 (MQ=255)
         gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCACGGAACGGCGCAATACCGGTGCCTGAGCCAATCAt                                                                                            >  1:243318/1‑100 (MQ=255)
           ttACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGTGCGTTGCTGCATAAAGGCGCGGAGCGGCGCAATACCGGTGCCTGGGCCAATCATaa                                                                                          <  1:1996374/100‑1 (MQ=255)
             aCCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAAt                                                                                         >  1:647287/1‑99 (MQ=255)
              ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCa                                                                                       >  2:1002252/1‑100 (MQ=255)
                aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                     >  2:1596610/1‑100 (MQ=255)
                   gAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCggg                                                                                  >  1:70731/1‑100 (MQ=255)
                    aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCgggg                                                                                 <  1:3497167/100‑1 (MQ=255)
                     aCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCgggg                                                                                 >  2:3162552/1‑99 (MQ=255)
                      cAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGt                                                                                >  1:577524/1‑99 (MQ=255)
                      cAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGtt                                                                               <  1:3841540/100‑1 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGttac                                                                             >  1:2850071/1‑98 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTc                                                                             >  2:1225123/1‑100 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTc                                                                             >  2:1677245/1‑100 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTc                                                                             >  1:2795414/1‑100 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTc                                                                             >  1:521750/1‑100 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTc                                                                             <  2:3883688/100‑1 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTc                                                                             <  2:3314797/100‑1 (MQ=255)
                          cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTg                                                                           >  1:3690940/1‑100 (MQ=255)
                          cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTg                                                                           >  2:1186769/1‑100 (MQ=255)
                            gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGa                                                                         <  1:1474583/100‑1 (MQ=255)
                                 tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggc                                                                    <  1:1150961/100‑1 (MQ=255)
                                 tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtgg                                                                     >  2:1283079/1‑99 (MQ=255)
                                 tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtg                                                                      <  2:2373738/98‑1 (MQ=255)
                                 tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAAGACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggc                                                                    <  2:1250123/100‑1 (MQ=255)
                                   ccTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctg                                                                  <  1:77244/100‑1 (MQ=255)
                                    cTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctgg                                                                 >  1:447669/1‑100 (MQ=255)
                                     tGGCGCTTCGTCGGCGGCGCTTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGAGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                                <  1:1394226/100‑1 (MQ=255)
                                     tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                                >  2:3020150/1‑100 (MQ=255)
                                      ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                               >  2:104720/1‑100 (MQ=255)
                                       ncgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAg                                                              >  1:3230424/2‑100 (MQ=255)
                                         gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCCGGc                                                            >  2:272397/1‑100 (MQ=255)
                                         gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGc                                                            <  2:1937924/100‑1 (MQ=255)
                                          cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCg                                                           >  2:2108611/1‑100 (MQ=255)
                                          cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCg                                                           <  1:1459885/100‑1 (MQ=255)
                                          cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCg                                                           <  1:3044890/100‑1 (MQ=255)
                                              gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTTGTTTGGCTGGCAGGCGGAAg                                                       <  2:1850445/100‑1 (MQ=255)
                                              gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                       <  2:2656198/100‑1 (MQ=255)
                                                cggcggCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGAGGaa                                                        >  1:1632772/1‑97 (MQ=255)
                                                 ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCAAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTa                                                    >  1:905198/1‑100 (MQ=255)
                                                 ggcggcCCGTCGCTGCATCAACGCGCGGAACGACGTAATGCCCGTGGGTGGGCCAATCATAATCACCGTGGTGTCTGGATTGCATGTCGGGCGGAAg                                                       <  2:3567864/97‑1 (MQ=255)
                                                  gcggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGATTCTGGATTGGCTGGCAGGCGGAAGTTAt                                                   >  1:1275711/1‑100 (MQ=255)
                                                   cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATc                                                  >  2:616931/1‑100 (MQ=255)
                                                   cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGATTGGCTGGCAGGCGGAAGTTATCg                                                 >  2:3371954/1‑100 (MQ=255)
                                                     gcgcGTTGCTGCATAAAGGCGCGGAACGGCGTAATATCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCg                                                 <  1:1535525/99‑1 (MQ=255)
                                                     gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                                <  1:2332190/100‑1 (MQ=255)
                                                      cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGtt                                               >  2:3682941/1‑100 (MQ=255)
                                                        cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTAt                                             >  2:3045818/1‑100 (MQ=255)
                                                        cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTATCTGGATTGGCTGGCAGGCGGAAGTTATCGTTAt                                             >  1:2039313/1‑100 (MQ=255)
                                                         gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGttn                                              <  2:2224587/98‑2 (MQ=255)
                                                            gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                         <  1:2785326/100‑1 (MQ=255)
                                                            gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                         <  1:3198641/100‑1 (MQ=255)
                                                            gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                         <  2:1469784/100‑1 (MQ=255)
                                                            gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTACCGTTATGTTc                                         >  2:2592080/1‑100 (MQ=255)
                                                            gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                         <  1:3403276/100‑1 (MQ=255)
                                                             ctgcATAAAGGCGCGGAACGGCGCAATTCCGTTGTTTTATTCTATCATTATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCg                                        <  1:1018071/100‑1 (MQ=255)
                                                              tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                         >  2:3277446/1‑98 (MQ=255)
                                                              tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                       >  1:2109131/1‑100 (MQ=255)
                                                               gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCTAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                      <  1:1104729/100‑1 (MQ=255)
                                                               gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                      >  2:2075611/1‑100 (MQ=255)
                                                               gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTATATGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                      >  1:1718004/1‑100 (MQ=255)
                                                                cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                     >  1:932368/1‑100 (MQ=255)
                                                                cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                     <  2:2955732/100‑1 (MQ=255)
                                                                cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                     <  2:3402097/100‑1 (MQ=255)
                                                                 aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                     <  1:540395/99‑1 (MQ=255)
                                                                  tAAAGGCGCGGAACGGCGCAATACCGGTTCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATaaa                                   >  2:1752048/1‑100 (MQ=255)
                                                                  tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATaaa                                   >  2:2488868/1‑100 (MQ=255)
                                                                   aaaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAAt                                  >  1:2140526/1‑100 (MQ=255)
                                                                    aaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                                 >  1:2694711/1‑100 (MQ=255)
                                                                        cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGAATGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                             >  1:1959273/1‑100 (MQ=255)
                                                                        cgcgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                             <  1:3046279/100‑1 (MQ=255)
                                                                         gcgGAACGGCGCAATACCGGTGGCTGGGCCAATCATAATCATCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgc                              <  2:2497555/98‑1 (MQ=255)
                                                                         gcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCgg                            <  1:857944/100‑1 (MQ=255)
                                                                         gcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCgg                            <  1:1971093/100‑1 (MQ=255)
                                                                              aCGGCGGAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTc                       <  2:885961/100‑1 (MQ=255)
                                                                                ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                     >  1:2473166/1‑100 (MQ=255)
                                                                                 gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGcc                    >  1:3530296/1‑100 (MQ=255)
                                                                                  cgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                   >  1:3836564/1‑100 (MQ=255)
                                                                                  cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                   <  1:1318425/100‑1 (MQ=255)
                                                                                   gcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCTATAAATACGCGGACTTCGCCct                  >  1:1619131/1‑100 (MQ=255)
                                                                                     aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGTGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctct                <  2:1756017/100‑1 (MQ=255)
                                                                                       tACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttc              <  2:662854/100‑1 (MQ=255)
                                                                                        aCCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttct             >  2:853551/1‑100 (MQ=255)
                                                                                        aCCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttc              <  2:900779/99‑1 (MQ=255)
                                                                                         ccGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTcttctt            <  1:1439421/100‑1 (MQ=255)
                                                                                         ccGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTcttctt            <  1:2259825/100‑1 (MQ=255)
                                                                                          cGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCttcttc           <  2:3336934/100‑1 (MQ=255)
                                                                                          cGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCttcttc           <  1:3351937/100‑1 (MQ=255)
                                                                                              gCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcg       <  1:2995612/100‑1 (MQ=255)
                                                                                                cTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgcg     >  1:3657513/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTC  >  NC_000913/2890422‑2890620

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