Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,499,2680TC25.0% 235.6 / 88.2 132intergenic (+67/‑180)cysG/yhfLfused siroheme synthase 1,3‑dimethyluroporphyriongen III dehydrogenase and siroheme ferrochelatase/uroporphyrinogen methyltransferase/small lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (46/53);  new base C (17/16);  total (63/69)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

Only 100 of 132 total aligned reads displayed.

CGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCC  >  NC_000913/3499171‑3499366
                                                                                                 |                                                                                                  
cgATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAg                                                                                                  >  1:915578/1‑100 (MQ=255)
 gATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGc                                                                                                 <  1:412055/100‑1 (MQ=255)
      actgaactgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCtt                                                                                            >  1:146437/1‑100 (MQ=255)
       ctgaactgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTta                                                                                           <  1:1908506/100‑1 (MQ=255)
       ctgaactgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTta                                                                                           >  2:2021894/1‑100 (MQ=255)
         gaactgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTtaat                                                                                         >  1:1493998/1‑100 (MQ=255)
           actgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTtaataa                                                                                       <  2:1987414/100‑1 (MQ=255)
            ctgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTaataaa                                                                                      >  2:361911/1‑100 (MQ=255)
            ctgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTaataaa                                                                                      >  2:1145110/1‑100 (MQ=255)
             tgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTaataaat                                                                                     <  2:2325153/100‑1 (MQ=255)
             tgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTaataaat                                                                                     >  2:2185857/1‑100 (MQ=255)
              ggTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTaataaata                                                                                    <  2:824979/100‑1 (MQ=255)
              ggTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTaataaata                                                                                    <  2:227887/100‑1 (MQ=255)
               gTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAatat                                                                                   <  2:869349/100‑1 (MQ=255)
               gTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAatat                                                                                   <  1:486893/100‑1 (MQ=255)
                 tctcCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAAtatta                                                                                 <  2:146435/100‑1 (MQ=255)
                   tcCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAAtattat                                                                                >  1:573011/1‑99 (MQ=255)
                    ccAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCt                                                                              >  2:1529978/1‑100 (MQ=255)
                        ccATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTa                                                                          <  1:2074399/100‑1 (MQ=255)
                           ttAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATg                                                                       <  1:401713/100‑1 (MQ=255)
                           ttAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATg                                                                       <  1:192007/100‑1 (MQ=255)
                             aaTTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATggg                                                                     <  2:1922180/100‑1 (MQ=255)
                              aTTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGt                                                                    <  2:1493994/100‑1 (MQ=255)
                                 tAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTtat                                                                 <  2:89755/100‑1 (MQ=255)
                                   aCCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTatat                                                               >  1:1666757/1‑100 (MQ=255)
                                   aCCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTatat                                                               <  2:42567/100‑1 (MQ=255)
                                    cccGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATAtt                                                              >  1:1953726/1‑100 (MQ=255)
                                     ccGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTc                                                             >  2:1566032/1‑100 (MQ=255)
                                     ccGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTc                                                             <  2:1391177/100‑1 (MQ=255)
                                       gggCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCAc                                                           <  2:1563066/100‑1 (MQ=255)
                                        ggCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCAcc                                                          >  1:230438/1‑100 (MQ=255)
                                        ggCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCAcc                                                          >  1:1126423/1‑100 (MQ=255)
                                         nccAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCg                                                         >  2:1169470/2‑100 (MQ=255)
                                          ccAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGt                                                        <  1:90275/100‑1 (MQ=255)
                                           cAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGtt                                                       <  1:924092/100‑1 (MQ=255)
                                            agagAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTa                                                      >  2:977115/1‑100 (MQ=255)
                                              agaATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTAc                                                     <  2:2154619/99‑1 (MQ=255)
                                              agaATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACa                                                    <  2:1701227/100‑1 (MQ=255)
                                                aaTTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGc                                                  <  2:395357/100‑1 (MQ=255)
                                                    cTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGAc                                               >  2:64769/1‑99 (MQ=255)
                                                    cTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACt                                              >  2:1341372/1‑100 (MQ=255)
                                                    cTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACt                                              <  2:727636/100‑1 (MQ=255)
                                                     tGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCAGCTAATGGGTTATATTCACCGTTACAGCGACtt                                             >  1:197637/1‑100 (MQ=255)
                                                      ggCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTa                                            >  1:956689/1‑100 (MQ=255)
                                                       gCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTAc                                           <  2:2214735/100‑1 (MQ=255)
                                                          ttcttAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACaaa                                        <  1:1313300/100‑1 (MQ=255)
                                                           tcttAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACaaaa                                       >  1:1427325/1‑100 (MQ=255)
                                                           tcttAATGACGCTTTTTATATATTCAACAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACaaa                                        <  2:95668/99‑1 (MQ=255)
                                                            cttAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAAt                                      <  1:245461/100‑1 (MQ=255)
                                                            cttAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAAt                                      >  2:183669/1‑100 (MQ=255)
                                                            cttAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAAt                                      >  1:25357/1‑100 (MQ=255)
                                                            cttAATGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACTGTTACAGCGACTTACAAAAt                                      <  2:1425875/100‑1 (MQ=255)
                                                             ttAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAAtt                                     <  1:2059751/100‑1 (MQ=255)
                                                              tAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTa                                    >  1:1944221/1‑100 (MQ=255)
                                                               aaTGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTaa                                   >  1:446599/1‑100 (MQ=255)
                                                                aTGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAg                                  >  1:2002292/1‑100 (MQ=255)
                                                                 tGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGt                                 <  1:1194985/100‑1 (MQ=255)
                                                                 tGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGt                                 <  2:969465/100‑1 (MQ=255)
                                                                 tGACTCTTTTTATATATTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGt                                 >  1:2010391/1‑100 (MQ=255)
                                                                  gACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTc                                <  1:1685805/100‑1 (MQ=255)
                                                                   aCTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCa                               <  1:137568/100‑1 (MQ=255)
                                                                   aCTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCa                               >  2:2094429/1‑100 (MQ=255)
                                                                    ctctTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCAt                              >  1:1297657/1‑100 (MQ=255)
                                                                     tctTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCAtt                             >  1:1752335/1‑100 (MQ=255)
                                                                     tctTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCAtt                             <  2:573010/100‑1 (MQ=255)
                                                                      ctTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTc                            >  2:1723747/1‑100 (MQ=255)
                                                                       tttttATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCg                           <  2:197635/100‑1 (MQ=255)
                                                                         tttATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCt                         >  1:825444/1‑100 (MQ=255)
                                                                          ttATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTc                        <  2:625728/100‑1 (MQ=255)
                                                                            atatatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCAc                      <  2:1354653/100‑1 (MQ=255)
                                                                            atatatTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCAc                      <  1:632470/100‑1 (MQ=255)
                                                                             tatatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACCTACAAAATTAAGTCATTCGCTCACt                     <  1:64771/100‑1 (MQ=255)
                                                                             tatatTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACt                     <  1:1529982/100‑1 (MQ=255)
                                                                              ntatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTg                    >  2:406164/2‑100 (MQ=255)
                                                                              atatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTg                    >  2:684560/1‑100 (MQ=255)
                                                                              atatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTg                    >  2:267987/1‑100 (MQ=255)
                                                                              atatTCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCtagatc                    >  1:2263837/1‑95 (MQ=255)
                                                                               tatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGa                   >  1:1540445/1‑100 (MQ=255)
                                                                               tatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGa                   >  1:479676/1‑100 (MQ=255)
                                                                               tatTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGa                   >  1:1777910/1‑100 (MQ=255)
                                                                                ttccgatctaTAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGAc                  <  2:1708678/91‑1 (MQ=255)
                                                                                atTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGAc                  >  2:11146/1‑100 (MQ=255)
                                                                                 ttCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACg                 <  1:1655121/100‑1 (MQ=255)
                                                                                 ttCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACg                 <  2:4436/100‑1 (MQ=255)
                                                                                 ttCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACg                 >  1:1922548/1‑100 (MQ=255)
                                                                                 ttCAGCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACg                 <  1:1424471/100‑1 (MQ=255)
                                                                                  tCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGc                <  1:1013014/100‑1 (MQ=255)
                                                                                    aGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCa              >  1:388701/1‑100 (MQ=255)
                                                                                    aGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCa              >  1:1729803/1‑100 (MQ=255)
                                                                                     gCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTAGATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCaa             <  2:230436/100‑1 (MQ=255)
                                                                                     gCAAATAAAACACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCaa             >  2:1886682/1‑100 (MQ=255)
                                                                                      cAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAAt            >  1:826989/1‑100 (MQ=255)
                                                                                       aaaTAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGATCACTGACGCCAAtt           >  1:2297448/1‑100 (MQ=255)
                                                                                         aTAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGAcagatcggaa         >  1:1708683/1‑91 (MQ=255)
                                                                                          tAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGATGCCAATTaaa        >  1:919462/1‑100 (MQ=255)
                                                                                           aaaaCATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAat       <  1:267989/100‑1 (MQ=255)
                                                                                             aaCACAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAatat     <  1:1145112/100‑1 (MQ=255)
                                                                                              aCATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAatat     <  1:183671/99‑1 (MQ=255)
                                                                                               cATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTc   >  2:434629/1‑100 (MQ=255)
                                                                                               cATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTc   >  1:1318156/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCC  >  NC_000913/3499171‑3499366

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: ‑