Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 882 254.4 12957588 98.0% 12698436 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,233,673 C→T G63G (GGC→GGT yjbE → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,233,6730CT96.1% 772.7 / 8.2 229G63G (GGC→GGTyjbEhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/7);  new base T (86/134);  total (88/141)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

Only 100 of 243 total aligned reads displayed.

GCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTG  >  NC_000913/4233584‑4233762
                                                                                          |                                                                                         
gcaACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGtggtggg                                                                                        >  2:4270046/1‑92 (MQ=255)
  aaCGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGa                                                                                      >  1:4986433/1‑93 (MQ=255)
  aaCGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGa                                                                                      <  2:3751671/93‑1 (MQ=255)
   aCGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGAt                                                                                     <  1:3601337/93‑1 (MQ=255)
   aCGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGAt                                                                                     <  2:390427/93‑1 (MQ=255)
     gTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATgg                                                                                   <  2:2830288/93‑1 (MQ=255)
     gTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATgg                                                                                   <  2:5332029/93‑1 (MQ=255)
      tCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGt                                                                                  >  2:2902906/1‑93 (MQ=255)
      tCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGt                                                                                  >  1:5721404/1‑93 (MQ=255)
        ggTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTc                                                                                <  1:423781/93‑1 (MQ=255)
        ggTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTc                                                                                <  2:601167/93‑1 (MQ=255)
        ggTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTc                                                                                >  2:3430248/1‑93 (MQ=255)
         gTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCt                                                                               >  1:3245336/1‑93 (MQ=255)
          tttCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTa                                                                              >  1:3847852/1‑93 (MQ=255)
          tttCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGGGGGTGTCGCGCTTGCGGCAACCGGTGGTGGGGATGGTTCTa                                                                              >  1:5278720/1‑93 (MQ=255)
           ttCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTaa                                                                             <  2:5068178/93‑1 (MQ=255)
            tCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGGTCTAAt                                                                            >  2:3616640/1‑93 (MQ=255)
              ggcggGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATAc                                                                          <  1:2975285/93‑1 (MQ=255)
              ggcggGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATAc                                                                          <  2:2437373/93‑1 (MQ=255)
              ggcggGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATAc                                                                          >  2:3110149/1‑93 (MQ=255)
                cggGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCg                                                                        <  2:3323537/93‑1 (MQ=255)
                cggGGAGTGCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCg                                                                        <  1:2338619/93‑1 (MQ=255)
                 gggtAGTGCATCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCgg                                                                       <  1:2861619/93‑1 (MQ=255)
                 ggggAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCgg                                                                       >  2:2183905/1‑93 (MQ=255)
                 ggggAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCgg                                                                       <  2:3131568/93‑1 (MQ=255)
                  gggAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCggg                                                                      >  1:2222589/1‑93 (MQ=255)
                   ggAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGa                                                                     >  1:1892794/1‑93 (MQ=255)
                      gTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGacca                                                                  >  2:3686431/1‑93 (MQ=255)
                      gTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTATAATACCGGGacca                                                                  >  2:5302050/1‑93 (MQ=255)
                       ttCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGaccac                                                                 >  1:6014287/1‑93 (MQ=255)
                       ttCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGaccac                                                                 <  2:4753051/93‑1 (MQ=255)
                        tcctccGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGaccacc                                                                >  1:631212/1‑93 (MQ=255)
                         cctccGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGAccacca                                                               <  1:452708/93‑1 (MQ=255)
                          ctccGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACcaccac                                                              <  2:4980189/93‑1 (MQ=255)
                             cgcgACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAAc                                                           >  1:5094840/1‑93 (MQ=255)
                               cgACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCa                                                         >  2:4608430/1‑93 (MQ=255)
                               cgACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCa                                                         <  2:5648527/93‑1 (MQ=255)
                                gACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCAc                                                        <  1:3161714/93‑1 (MQ=255)
                                gACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCAc                                                        >  2:4324406/1‑93 (MQ=255)
                                gACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCAc                                                        >  2:1649243/1‑93 (MQ=255)
                                 aCCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACg                                                       >  2:6431613/1‑93 (MQ=255)
                                 aCCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACg                                                       <  1:3678686/93‑1 (MQ=255)
                                   cAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGac                                                     <  2:1999837/93‑1 (MQ=255)
                                    aGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGacc                                                    <  1:1598764/93‑1 (MQ=255)
                                     gCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGacca                                                   <  1:57769/93‑1 (MQ=255)
                                     gCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGacca                                                   >  2:5797872/1‑93 (MQ=255)
                                     gCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGacca                                                   >  2:4442038/1‑93 (MQ=255)
                                          aGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGt                                              >  1:2676177/1‑93 (MQ=255)
                                           gCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTa                                             <  2:5559443/93‑1 (MQ=255)
                                            cACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTAc                                            <  1:1425218/93‑1 (MQ=255)
                                            cACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTAc                                            >  2:1559196/1‑93 (MQ=255)
                                             aCCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTAcc                                           <  1:791799/93‑1 (MQ=255)
                                              ccGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTAccc                                          >  1:3667261/1‑93 (MQ=255)
                                               cGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATCCCGGGACCACCACAACCACGACCACCAGTACCCa                                         <  1:687066/93‑1 (MQ=255)
                                               cGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCa                                         <  1:6221018/93‑1 (MQ=255)
                                               cGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCa                                         <  1:6082346/93‑1 (MQ=255)
                                               cGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCa                                         <  1:4433713/93‑1 (MQ=255)
                                                gTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAg                                        <  2:2139686/93‑1 (MQ=255)
                                                  aaGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGta                                      <  1:2526808/93‑1 (MQ=255)
                                                   aGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGtaa                                     <  1:5310800/93‑1 (MQ=255)
                                                    gCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGtaat                                    <  2:1397835/93‑1 (MQ=255)
                                                     cTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGtaata                                   <  2:544393/93‑1 (MQ=255)
                                                     cTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGCTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGtaata                                   <  2:931853/93‑1 (MQ=255)
                                                      tCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGtaataa                                  >  2:3076987/1‑93 (MQ=255)
                                                       cggcggTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATaaa                                 <  2:6148029/93‑1 (MQ=255)
                                                        ggcggTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAg                                <  1:4537943/93‑1 (MQ=255)
                                                         gcggTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgt                               >  1:3377862/1‑93 (MQ=255)
                                                         gcggTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgt                               <  2:2331450/93‑1 (MQ=255)
                                                         gcggTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgt                               <  1:2713558/93‑1 (MQ=255)
                                                         gcggTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACAGGGACCACCACAACACCGACTACCAGTACCCAGTAATAAAgt                               >  1:5560157/1‑93 (MQ=255)
                                                          cggTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgta                              <  2:39853/93‑1 (MQ=255)
                                                           ggtgggtgGCGCGATTGGGCAAGCCGTTGTTGGTGATGGTTATAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtat                             <  1:877448/93‑1 (MQ=255)
                                                            gtgggtgTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatg                            <  1:970649/93‑1 (MQ=255)
                                                             tggttgTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgt                           <  2:4845141/93‑1 (MQ=255)
                                                             tgggtgTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgt                           <  2:5241680/93‑1 (MQ=255)
                                                               ggtgTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgtat                         <  1:3120619/93‑1 (MQ=255)
                                                               ggtgTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgtat                         >  1:291829/1‑93 (MQ=255)
                                                               ggtgTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgtat                         <  1:4198375/93‑1 (MQ=255)
                                                                gtgtCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATc                        <  2:3453815/93‑1 (MQ=255)
                                                                     gcgcTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCaa                   >  2:2628609/1‑93 (MQ=255)
                                                                      cgcTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCaaa                  >  2:492364/1‑93 (MQ=255)
                                                                       gcTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCaaaa                 <  2:6280388/93‑1 (MQ=255)
                                                                        cTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACCACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAaat                >  2:541589/1‑93 (MQ=255)
                                                                        cTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAaat                >  1:4095860/1‑93 (MQ=255)
                                                                          tGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAaataa              >  2:6058908/1‑93 (MQ=255)
                                                                            cGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAAtt            <  1:2265947/93‑1 (MQ=255)
                                                                            cGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAAtt            <  2:5918621/93‑1 (MQ=255)
                                                                             ggCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTc           <  2:4709027/93‑1 (MQ=255)
                                                                              gCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCg          >  2:1235049/1‑93 (MQ=255)
                                                                              gCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCg          >  2:2761130/1‑93 (MQ=255)
                                                                               cAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACACCCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGa         <  2:3988629/93‑1 (MQ=255)
                                                                                aaCCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAg        <  2:5391895/93‑1 (MQ=255)
                                                                                aaCCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAg        >  1:1711836/1‑93 (MQ=255)
                                                                                aaCCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAg        >  1:1045094/1‑93 (MQ=255)
                                                                                 aCCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGt       >  1:6426548/1‑93 (MQ=255)
                                                                                   ccGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTc      <  1:1775681/93‑1 (MQ=255)
                                                                                     ggtggtggtGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCAt    <  1:1256427/93‑1 (MQ=255)
                                                                                       tggtggTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTg  >  2:3296003/1‑93 (MQ=255)
                                                                                       tggtggTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTg  <  1:2182410/93‑1 (MQ=255)
                                                                                       tggtggTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAATCATTg  >  2:4236589/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTG  >  NC_000913/4233584‑4233762

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATC  >  NC_000913/4233580‑4233766
                                                                                             |                                                                                             
GGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGG                                                                                               >  SRR4302537.1340294/1‑94 (MQ=60)
 GCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTG                                                                                              <  SRR4302537.983387/94‑1 (MQ=60)
  CAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCCCCGTAAGCTCGGCGGTGGGGGTCGCGCTTGGGGCAACCGGTGGGGG                                                                                             >  SRR4302537.585892/1‑94 (MQ=60)
  CAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGG                                                                                             <  SRR4302537.3816257/94‑1 (MQ=60)
   AGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCCACCGGTGGCGGG                                                                                            >  SRR4302537.2979026/1‑94 (MQ=60)
     CAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGA                                                                                          >  SRR4302537.5020799/1‑94 (MQ=60)
       ACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATG                                                                                        <  SRR4302537.3626720/94‑1 (MQ=60)
         GTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGT                                                                                      >  SRR4302537.5760572/1‑94 (MQ=60)
         GTCGGTTTCGGCGGGGAGGTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGT                                                                                      <  SRR4302537.3915742/94‑1 (MQ=60)
            GGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCT                                                                                   <  SRR4302537.426722/94‑1 (MQ=60)
            GGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCT                                                                                   >  SRR4302537.3268206/1‑94 (MQ=60)
             GTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTA                                                                                  >  SRR4302537.5455584/1‑94 (MQ=60)
             GTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTA                                                                                  >  SRR4302537.3874793/1‑94 (MQ=60)
             GTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGGGGGTGTCGCGCTTGCGGCAACCGGTGGTGGGGATGGTTCTA                                                                                  >  SRR4302537.5315036/1‑94 (MQ=60)
              TTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAA                                                                                 >  SRR4302537.2244808/1‑94 (MQ=60)
                 CGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATAC                                                                              <  SRR4302537.4478434/94‑1 (MQ=60)
                 CGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATAC                                                                              <  SRR4302537.2179685/94‑1 (MQ=60)
                  GGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACC                                                                             <  SRR4302537.2996117/94‑1 (MQ=60)
                  GGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACC                                                                             <  SRR4302537.1931032/94‑1 (MQ=60)
                  GGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACC                                                                             >  SRR4302537.3627289/1‑94 (MQ=60)
                   GCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCG                                                                            >  SRR4302537.594519/1‑94 (MQ=60)
                    CGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGG                                                                           >  SRR4302537.2571696/1‑94 (MQ=60)
                    CGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGG                                                                           <  SRR4302537.532005/94‑1 (MQ=60)
                    CGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGG                                                                           >  SRR4302537.2472252/1‑94 (MQ=60)
                    CGGGGAGTGCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGG                                                                           <  SRR4302537.2355105/94‑1 (MQ=60)
                     gggtagtgcaTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGG                                                                          <  SRR4302537.2881684/84‑1 (MQ=60)
                     GGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGG                                                                          >  SRR4302537.867313/1‑94 (MQ=60)
                     GGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGG                                                                          <  SRR4302537.6141126/94‑1 (MQ=60)
                     GGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGG                                                                          >  SRR4302537.2238247/1‑94 (MQ=60)
                      GGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGA                                                                         >  SRR4302537.1905950/1‑94 (MQ=60)
                          GTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCAC                                                                     >  SRR4302537.6055607/1‑94 (MQ=60)
                           TTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACC                                                                    >  SRR4302537.635504/1‑94 (MQ=60)
                             CCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCAC                                                                  <  SRR4302537.455840/94‑1 (MQ=60)
                               TCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGAGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTCGTGGTTCTAATACCGGGcacacaacaa                                                                >  SRR4302537.2884781/1‑84 (MQ=60)
                                CCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAAC                                                               >  SRR4302537.5129925/1‑94 (MQ=60)
                                CCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAAC                                                               >  SRR4302537.3154723/1‑94 (MQ=60)
                                  GCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCA                                                             <  SRR4302537.325950/94‑1 (MQ=60)
                                    GACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACG                                                           <  SRR4302537.2319594/94‑1 (MQ=60)
                                    GACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACG                                                           <  SRR4302537.1303650/94‑1 (MQ=60)
                                    GACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACG                                                           <  SRR4302537.3183908/94‑1 (MQ=60)
                                     ACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGA                                                          <  SRR4302537.3704567/94‑1 (MQ=60)
                                       CAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACC                                                        <  SRR4302537.2923921/94‑1 (MQ=60)
                                        AGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCA                                                       >  SRR4302537.3092337/1‑94 (MQ=60)
                                        AGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCA                                                       <  SRR4302537.1609937/94‑1 (MQ=60)
                                         GCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCAC                                                      >  SRR4302537.2869592/1‑94 (MQ=60)
                                         GCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCAC                                                      <  SRR4302537.58195/94‑1 (MQ=60)
                                          CGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACC                                                     >  SRR4302537.517165/1‑94 (MQ=60)
                                           GTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCA                                                    >  SRR4302537.5781142/1‑94 (MQ=60)
                                            TCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAG                                                   <  SRR4302537.2714737/94‑1 (MQ=60)
                                            TCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAG                                                   <  SRR4302537.3535533/94‑1 (MQ=60)
                                             CAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGT                                                  >  SRR4302537.2694866/1‑94 (MQ=60)
                                             CAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGT                                                  <  SRR4302537.713240/94‑1 (MQ=60)
                                                CACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACC                                               <  SRR4302537.1435141/94‑1 (MQ=60)
                                                 ACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCC                                              >  SRR4302537.1720209/1‑94 (MQ=60)
                                                 ACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCC                                              <  SRR4302537.797233/94‑1 (MQ=60)
                                                 ACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCC                                              >  SRR4302537.3693069/1‑94 (MQ=60)
                                                   CGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATCCCGGGACCACCACAACCACGACCACCAGTACCCAG                                            <  SRR4302537.691776/94‑1 (MQ=60)
                                                   CGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAG                                            <  SRR4302537.4464330/94‑1 (MQ=60)
                                                   CGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAG                                            <  SRR4302537.1757458/94‑1 (MQ=60)
                                                   CGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAG                                            <  SRR4302537.5879502/94‑1 (MQ=60)
                                                   CGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAG                                            <  SRR4302537.6124134/94‑1 (MQ=60)
                                                   CGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAG                                            <  SRR4302537.6263845/94‑1 (MQ=60)
                                                    GTAAGCTCGGCGCTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGT                                           <  SRR4302537.875753/94‑1 (MQ=60)
                                                     TAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTA                                          <  SRR4302537.3078456/94‑1 (MQ=60)
                                                      AAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAA                                         <  SRR4302537.2544538/94‑1 (MQ=60)
                                                       AGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAAT                                        <  SRR4302537.5347283/94‑1 (MQ=60)
                                                       AGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAAT                                        >  SRR4302537.1873745/1‑94 (MQ=60)
                                                        GCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATA                                       <  SRR4302537.3282112/94‑1 (MQ=60)
                                                         CTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAA                                      <  SRR4302537.659070/94‑1 (MQ=60)
                                                            GGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGT                                   >  SRR4302537.3401801/1‑94 (MQ=60)
                                                            GGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGT                                   <  SRR4302537.4569360/94‑1 (MQ=60)
                                                            GGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGT                                   <  SRR4302537.5147619/94‑1 (MQ=60)
                                                            GGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACAGGGACCACCACAACACCGACTACCAGTACCCAGTAATAAAGT                                   >  SRR4302537.5598150/1‑94 (MQ=60)
                                                            GGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGT                                   <  SRR4302537.6125626/94‑1 (MQ=60)
                                                             GCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTA                                  <  SRR4302537.2732506/94‑1 (MQ=60)
                                                              CGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTAT                                 <  SRR4302537.369923/94‑1 (MQ=60)
                                                              CGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTAT                                 <  SRR4302537.5508381/94‑1 (MQ=60)
                                                              CGGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTAT                                 >  SRR4302537.3229999/1‑94 (MQ=60)
                                                               GGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATG                                >  SRR4302537.5054129/1‑94 (MQ=60)
                                                               GGTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATG                                >  SRR4302537.5628725/1‑94 (MQ=60)
                                                                GTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGT                               <  SRR4302537.977340/94‑1 (MQ=60)
                                                                GTGGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGT                               <  SRR4302537.6337033/94‑1 (MQ=60)
                                                                  GGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTAT                             >  SRR4302537.293860/1‑94 (MQ=60)
                                                                  GGGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTAT                             >  SRR4302537.1206774/1‑94 (MQ=60)
                                                                   GGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATC                            >  SRR4302537.4156790/1‑94 (MQ=60)
                                                                   GGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATC                            <  SRR4302537.4227688/94‑1 (MQ=60)
                                                                   GGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATC                            <  SRR4302537.3142440/94‑1 (MQ=60)
                                                                   GGTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATC                            <  SRR4302537.2218265/94‑1 (MQ=60)
                                                                    GTGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCC                           <  SRR4302537.1512473/94‑1 (MQ=60)
                                                                     TGTCGCGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCCCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCC                          <  SRR4302537.269303/94‑1 (MQ=60)
                                                                          CGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAA                     <  SRR4302537.3833818/94‑1 (MQ=60)
                                                                          CGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAA                     >  SRR4302537.1746355/1‑94 (MQ=60)
                                                                          CGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAA                     >  SRR4302537.2788035/1‑94 (MQ=60)
                                                                          CGCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAA                     <  SRR4302537.6435055/94‑1 (MQ=60)
                                                                           GCTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAAT                    >  SRR4302537.4124460/1‑94 (MQ=60)
                                                                            CTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATA                   >  SRR4302537.927232/1‑94 (MQ=60)
                                                                            CTTGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATA                   <  SRR4302537.4875303/94‑1 (MQ=60)
                                                                              TGCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCATTAATAAAGTATGTATCCCCAAAATAAT                 <  SRR4302537.4185505/94‑1 (MQ=60)
                                                                               GCGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATT                <  SRR4302537.4964880/94‑1 (MQ=60)
                                                                                CGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTC               <  SRR4302537.2281950/94‑1 (MQ=60)
                                                                                CGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTC               <  SRR4302537.4568702/94‑1 (MQ=60)
                                                                                CGGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTC               <  SRR4302537.6398632/94‑1 (MQ=60)
                                                                                 GGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGAGCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCG              >  SRR4302537.3725267/1‑94 (MQ=60)
                                                                                 GGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCG              <  SRR4302537.606323/94‑1 (MQ=60)
                                                                                 GGCAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCG              <  SRR4302537.146969/94‑1 (MQ=60)
                                                                                   CAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAG            >  SRR4302537.1052306/1‑94 (MQ=60)
                                                                                   CAACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAG            >  SRR4302537.1723791/1‑94 (MQ=60)
                                                                                    AACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGT           >  SRR4302537.6470218/1‑94 (MQ=60)
                                                                                    AACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGT           <  SRR4302537.4455914/94‑1 (MQ=60)
                                                                                     ACCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCACTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTC          <  SRR4302537.4339697/94‑1 (MQ=60)
                                                                                      CCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCA         <  SRR4302537.1788044/94‑1 (MQ=60)
                                                                                      CCGGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCA         >  SRR4302537.2310749/1‑94 (MQ=60)
                                                                                        GGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATT       <  SRR4302537.1265170/94‑1 (MQ=60)
                                                                                        GGTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATT       <  SRR4302537.4255132/94‑1 (MQ=60)
                                                                                         GTGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTG      <  SRR4302537.355548/94‑1 (MQ=60)
                                                                                          TGGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGC     <  SRR4302537.2197785/94‑1 (MQ=60)
                                                                                           GGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCA    >  SRR4302537.6411659/1‑94 (MQ=60)
                                                                                           GGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCA    >  SRR4302537.6440700/1‑94 (MQ=60)
                                                                                           GGTGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCA    <  SRR4302537.3189802/94‑1 (MQ=60)
                                                                                             TGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATC  <  SRR4302537.1784603/94‑1 (MQ=60)
                                                                                             TGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATC  <  SRR4302537.308780/94‑1 (MQ=60)
                                                                                             |                                                                                             
GGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATC  >  NC_000913/4233580‑4233766

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