Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,584,362 A→G intergenic (+53/‑42) yrhA → / → yrhD putative uncharacterized protein YrhA/uncharacterized protein YrhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,584,3620AG99.7% 1514.6 / ‑5.0 341intergenic (+53/‑42)yrhA/yrhDputative uncharacterized protein YrhA/uncharacterized protein YrhD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (186/154);  minor base T (1/0);  total (187/154)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 341 total aligned reads displayed.

GATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGT  >  NC_000913/3584265‑3584459
                                                                                                 |                                                                                                 
gATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCca                                                                                                >  2:5631569/1‑101 (MQ=255)
gATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCca                                                                                                <  2:5630756/101‑1 (MQ=255)
 aTTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCcag                                                                                               <  1:5631266/101‑1 (MQ=255)
  tttCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCcagc                                                                                              >  1:5631575/1‑101 (MQ=255)
  tttCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCcagc                                                                                              <  1:5631122/101‑1 (MQ=255)
    tctcTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCcagcag                                                                                            <  1:5630655/101‑1 (MQ=255)
    tctcTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCcagcag                                                                                            >  2:5631577/1‑101 (MQ=255)
    tctcTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCcagcag                                                                                            <  1:5631065/101‑1 (MQ=255)
        ttCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGaaaa                                                                                        >  1:5631581/1‑101 (MQ=255)
        ttCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGaaaa                                                                                        <  1:5631076/101‑1 (MQ=255)
         tCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAAt                                                                                       >  2:5631583/1‑101 (MQ=255)
          cATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAAtt                                                                                      <  1:5630751/101‑1 (MQ=255)
          cATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAAtt                                                                                      <  1:5631205/101‑1 (MQ=255)
           aTTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTc                                                                                     >  2:5631587/1‑101 (MQ=255)
             tCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCt                                                                                   <  2:5630862/101‑1 (MQ=255)
             tCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCt                                                                                   <  1:5631093/101‑1 (MQ=255)
              cTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCtt                                                                                  <  2:5630945/101‑1 (MQ=255)
                 aaaCGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGtt                                                                               <  1:5630914/101‑1 (MQ=255)
                 aaaCGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGtt                                                                               >  2:5631594/1‑101 (MQ=255)
                  aaCGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTc                                                                              >  2:5631597/1‑101 (MQ=255)
                  aaCGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTc                                                                              >  1:5631598/1‑101 (MQ=255)
                    cGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGa                                                                            >  2:5631602/1‑101 (MQ=255)
                     gAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAGAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGaa                                                                           <  1:5631105/101‑1 (MQ=255)
                     gAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGaa                                                                           >  1:5631607/1‑101 (MQ=255)
                      aaaTCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGaaa                                                                          <  1:5631136/101‑1 (MQ=255)
                      aaaTCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGaaa                                                                          >  2:5631608/1‑101 (MQ=255)
                        aTCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAAtt                                                                        <  1:5631030/101‑1 (MQ=255)
                         tCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAAttt                                                                       <  2:5631018/101‑1 (MQ=255)
                         tCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAAttt                                                                       >  1:5631612/1‑101 (MQ=255)
                          cATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTa                                                                      >  2:5631614/1‑101 (MQ=255)
                               aCAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGt                                                                 <  1:5630942/101‑1 (MQ=255)
                                cAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTg                                                                <  1:5631047/101‑1 (MQ=255)
                                cAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTg                                                                >  2:5631619/1‑101 (MQ=255)
                                 aGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGc                                                               <  2:5630904/101‑1 (MQ=255)
                                 aGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGc                                                               >  1:5631622/1‑101 (MQ=255)
                                   ttgttCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCaa                                                             >  1:5631627/1‑101 (MQ=255)
                                   ttgttCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCaa                                                             <  1:5631151/101‑1 (MQ=255)
                                    tgttCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAg                                                            >  2:5631628/1‑101 (MQ=255)
                                      ttCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCt                                                          >  2:5631632/1‑101 (MQ=255)
                                        cATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTAt                                                        >  1:5631636/1‑101 (MQ=255)
                                          taataaTCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATgg                                                      >  2:5631641/1‑101 (MQ=255)
                                           aataatCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGa                                                     >  1:5631643/1‑101 (MQ=255)
                                            ataatCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGaa                                                    >  1:5631646/1‑101 (MQ=255)
                                             taatCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGaaa                                                   >  2:5631651/1‑101 (MQ=255)
                                              aatCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGaaaa                                                  >  1:5631655/1‑101 (MQ=255)
                                               atCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGaaaaa                                                 <  1:5631073/101‑1 (MQ=255)
                                               atCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGaaaaa                                                 >  2:5631659/1‑101 (MQ=255)
                                               atCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGaaaaa                                                 >  2:5631656/1‑101 (MQ=255)
                                                tCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAAc                                                >  2:5631664/1‑101 (MQ=255)
                                                 cAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACt                                               <  2:5631324/101‑1 (MQ=255)
                                                 cAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACt                                               <  1:5631010/101‑1 (MQ=255)
                                                  aaaaGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACtt                                              >  2:5631668/1‑101 (MQ=255)
                                                  aaaaGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACtt                                              <  2:5631034/101‑1 (MQ=255)
                                                   aaaGGAAACAAAAATCATTGTAAACCATTCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTa                                             >  1:5631669/1‑101 (MQ=255)
                                                        aaaCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACg                                        <  2:5630958/101‑1 (MQ=255)
                                                         aaCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACgg                                       >  2:5631676/1‑101 (MQ=255)
                                                          aCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGa                                      <  1:5631170/101‑1 (MQ=255)
                                                          aCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGa                                      >  1:5631677/1‑101 (MQ=255)
                                                          aCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGa                                      <  1:5631024/101‑1 (MQ=255)
                                                           cAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGaa                                     <  2:5631041/101‑1 (MQ=255)
                                                           cAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGaa                                     <  2:5631150/101‑1 (MQ=255)
                                                           cAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGaa                                     <  1:5630973/101‑1 (MQ=255)
                                                            aaaaaTCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGAAAc                                    >  2:5631681/1‑101 (MQ=255)
                                                             aaaaTCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACt                                   >  1:5631685/1‑101 (MQ=255)
                                                              aaaTCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTa                                  >  2:5631687/1‑101 (MQ=255)
                                                               aaTCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTac                                 >  1:5631689/1‑101 (MQ=255)
                                                                 tCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTacac                               <  2:5631060/101‑1 (MQ=255)
                                                                 tCATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTacac                               <  2:5631008/101‑1 (MQ=255)
                                                                  cATTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACt                              >  2:5631696/1‑101 (MQ=255)
                                                                   aTTGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTc                             >  1:5631698/1‑101 (MQ=255)
                                                                    ttGTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTcc                            >  2:5631700/1‑101 (MQ=255)
                                                                      gTAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCtt                          <  2:5631117/101‑1 (MQ=255)
                                                                       tAAACCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTa                         >  1:5631705/1‑101 (MQ=255)
                                                                         aaCCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAg                       >  2:5631710/1‑101 (MQ=255)
                                                                         aaCCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAg                       >  1:5631711/1‑101 (MQ=255)
                                                                          aCCATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGt                      >  2:5631714/1‑101 (MQ=255)
                                                                            cATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTga                    >  1:5631720/1‑101 (MQ=255)
                                                                            cATGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTga                    >  2:5631718/1‑101 (MQ=255)
                                                                              tGCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTgaga                  <  1:5630996/101‑1 (MQ=255)
                                                                               gCCATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGAt                 >  2:5631724/1‑101 (MQ=255)
                                                                                ccATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATg                <  2:5631172/101‑1 (MQ=255)
                                                                                ccATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATg                <  2:5631481/101‑1 (MQ=255)
                                                                                 cATCTTTTTAGAGATTGCCAGCAGAAATTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATgg               >  2:5631727/1‑101 (MQ=255)
                                                                                 cATCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATgg               <  2:5631083/101‑1 (MQ=255)
                                                                                  aTCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGa              <  2:5631104/101‑1 (MQ=255)
                                                                                  aTCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGa              >  1:4011121/1‑101 (MQ=255)
                                                                                   tCTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGAc             <  1:5631331/101‑1 (MQ=255)
                                                                                    cTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGAcc            >  1:5631730/1‑101 (MQ=255)
                                                                                    cTTTTTAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGAcc            <  1:5630845/101‑1 (MQ=255)
                                                                                     tttttAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCg           <  2:5631353/101‑1 (MQ=255)
                                                                                     tttttAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCg           >  1:5631736/1‑101 (MQ=255)
                                                                                        ttAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAg        >  2:5631738/1‑101 (MQ=255)
                                                                                         tAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTTCACTCCTTaa                        <  2:5630984/84‑1 (MQ=255)
                                                                                         tAGAGATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGa       <  2:5631111/101‑1 (MQ=255)
                                                                                          agagATTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGAc      >  1:5631742/1‑101 (MQ=255)
                                                                                           gagaTTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACa     >  2:5631743/1‑101 (MQ=255)
                                                                                            agaTTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACAt    >  1:5631748/1‑101 (MQ=255)
                                                                                            agaTTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACAt    <  2:5631141/101‑1 (MQ=255)
                                                                                             gaTTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATg   <  1:5631491/101‑1 (MQ=255)
                                                                                             gaTTGCCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATg   >  1:5631753/1‑101 (MQ=255)
                                                                                                 |                                                                                                 
GATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGT  >  NC_000913/3584265‑3584459

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