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 593 264.2 13311096 97.8% 13018251 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,807,191 G→A P23S (CCT→TCT)  b1726 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,807,1910GA90.0% 639.3 / 41.9 201P23S (CCT→TCT) b1726hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/7);  new base A (97/84);  total (110/91)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 203 total aligned reads displayed.

AGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCG  >  NC_000913/1807101‑1807283
                                                                                          |                                                                                            
aGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAta                                                                                            >  2:164057/1‑93 (MQ=255)
aGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAta                                                                                            >  2:1673780/1‑93 (MQ=255)
aGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAta                                                                                            <  2:3340980/93‑1 (MQ=255)
 gTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGtat                                                                                           >  2:4065962/1‑93 (MQ=255)
  tCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAtata                                                                                          <  1:730631/93‑1 (MQ=255)
   cGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaa                                                                                         <  2:6333615/93‑1 (MQ=255)
   cGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaa                                                                                         <  2:5860831/93‑1 (MQ=255)
   cGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaa                                                                                         >  1:49793/1‑93 (MQ=255)
    gAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaa                                                                                        <  2:4244126/93‑1 (MQ=255)
    gAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaa                                                                                        >  1:2922825/1‑93 (MQ=255)
     aGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaaa                                                                                       >  2:2198185/1‑93 (MQ=255)
     aGCATGACCGCAATAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaaa                                                                                       <  2:4743737/93‑1 (MQ=255)
      gCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaaaa                                                                                      >  1:6518913/1‑93 (MQ=255)
       cATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAAt                                                                                     <  2:4103065/93‑1 (MQ=255)
         tGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCa                                                                                   >  1:2251476/1‑93 (MQ=255)
           aCCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacc                                                                                 >  2:1793653/1‑93 (MQ=255)
            ccGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCaccc                                                                                <  2:3115578/93‑1 (MQ=255)
              gCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaa                                                                              >  1:5839373/1‑93 (MQ=255)
              gCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaa                                                                              <  1:3111156/93‑1 (MQ=255)
               cATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaac                                                                             >  1:2501861/1‑93 (MQ=255)
                 ttAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaaccc                                                                           >  1:1480006/1‑93 (MQ=255)
                     tGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCaaaa                                                                       >  1:2017505/1‑93 (MQ=255)
                        cttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATcc                                                                    <  1:1427802/93‑1 (MQ=255)
                        cttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATct                                                                    >  2:1877622/1‑92 (MQ=255)
                        cttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATcc                                                                    <  1:542609/93‑1 (MQ=255)
                         ttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCt                                                                   >  1:1160626/1‑93 (MQ=255)
                         ttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCt                                                                   >  2:268120/1‑93 (MQ=255)
                          tccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCtt                                                                  >  1:3129303/1‑93 (MQ=255)
                          tccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCtt                                                                  <  2:4324042/93‑1 (MQ=255)
                          tccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCtt                                                                  >  2:2877056/1‑93 (MQ=255)
                            cttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCtttt                                                                <  2:3875999/93‑1 (MQ=255)
                             ttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTc                                                               <  2:2072987/93‑1 (MQ=255)
                                 tttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCaaaa                                                           <  2:4224908/93‑1 (MQ=255)
                                 tttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCaaaa                                                           <  1:5465005/93‑1 (MQ=255)
                                  ttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAAc                                                          >  2:5398640/1‑93 (MQ=255)
                                  ttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAAc                                                          <  2:5436880/93‑1 (MQ=255)
                                     cTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAACGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGc                                                       >  2:2500843/1‑93 (MQ=255)
                                       gCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCaa                                                     >  1:5551192/1‑93 (MQ=255)
                                        ccGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAAt                                                    >  1:5768010/1‑93 (MQ=255)
                                        ccGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAAt                                                    >  1:2602715/1‑93 (MQ=255)
                                         cGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATg                                                   <  1:4877177/93‑1 (MQ=255)
                                         cGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATg                                                   <  2:5320974/93‑1 (MQ=255)
                                         cGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATg                                                   <  1:2663463/93‑1 (MQ=255)
                                          gggTATTCATACATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGc                                                  >  1:2115287/1‑93 (MQ=255)
                                          gggTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGc                                                  <  2:708007/93‑1 (MQ=255)
                                          gggTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGc                                                  >  1:5496771/1‑93 (MQ=255)
                                              aTTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTcc                                              <  1:4077465/93‑1 (MQ=255)
                                              aTTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTcc                                              <  1:9206/93‑1 (MQ=255)
                                               ttCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTGCa                                             >  1:3845627/1‑93 (MQ=255)
                                                  aTAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCt                                          <  2:2468969/93‑1 (MQ=255)
                                                  aTAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCt                                          >  2:5225621/1‑93 (MQ=255)
                                                    aaaTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttg                                        >  2:2104503/1‑93 (MQ=255)
                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  1:4259269/1‑93 (MQ=255)
                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      <  2:3763764/93‑1 (MQ=255)
                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  1:5251072/1‑93 (MQ=255)
                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  1:5372869/1‑93 (MQ=255)
                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  2:4476937/1‑93 (MQ=255)
                                                       tttCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgttg                                     >  2:1319609/1‑93 (MQ=255)
                                                         tCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGAt                                   >  2:4900791/1‑93 (MQ=255)
                                                         tCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGAt                                   <  1:3236236/93‑1 (MQ=255)
                                                          cAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATa                                  >  1:1897055/1‑93 (MQ=255)
                                                          cAGCAAAGAAATCAATGTAGAGATAAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATa                                  >  1:2474102/1‑93 (MQ=255)
                                                           aGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATAc                                 <  2:4577132/93‑1 (MQ=255)
                                                            gCAAAGAAGTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACg                                >  1:1894019/1‑93 (MQ=255)
                                                              aaaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTc                              >  1:4744578/1‑93 (MQ=255)
                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAGTCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             <  2:4666285/93‑1 (MQ=255)
                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             >  1:4120840/1‑93 (MQ=255)
                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             <  1:4026708/93‑1 (MQ=255)
                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             <  1:5093756/93‑1 (MQ=255)
                                                                aGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCAt                            <  2:4088311/93‑1 (MQ=255)
                                                                 gAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATg                           >  1:5418381/1‑93 (MQ=255)
                                                                 gAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATg                           >  2:1213233/1‑93 (MQ=255)
                                                                  aaaTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGa                          <  2:3514787/93‑1 (MQ=255)
                                                                  aaaTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGa                          >  1:5698544/1‑93 (MQ=255)
                                                                    aTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAACTCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGt                        <  2:3421650/93‑1 (MQ=255)
                                                                       aaTGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTtct                     >  1:5889404/1‑93 (MQ=255)
                                                                       aaTGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTtct                     >  1:1833189/1‑93 (MQ=255)
                                                                       aaTGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTtct                     >  1:1660304/1‑93 (MQ=255)
                                                                         tGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTcc                   >  1:3668021/1‑93 (MQ=255)
                                                                         tGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTcc                   <  2:5761855/93‑1 (MQ=255)
                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  1:100436/93‑1 (MQ=255)
                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  1:5272157/93‑1 (MQ=255)
                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  2:6083448/93‑1 (MQ=255)
                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  1:3332637/93‑1 (MQ=255)
                                                                            agagATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTg                >  2:6056845/1‑93 (MQ=255)
                                                                             gagaTCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGa               >  1:1116922/1‑93 (MQ=255)
                                                                              agaTCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGAt              >  2:6282135/1‑93 (MQ=255)
                                                                               gaTCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGAtt             >  2:3136088/1‑93 (MQ=255)
                                                                                 tCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgttg                                     <  1:124567/67‑1 (MQ=255)
                                                                                 tCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTaa           >  1:642441/1‑93 (MQ=255)
                                                                                  cAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAg          >  1:5717145/1‑93 (MQ=255)
                                                                                  cAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAg          <  2:1590408/93‑1 (MQ=255)
                                                                                   aaTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAATTTGCAATTCGTTCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGa         >  2:6640544/1‑93 (MQ=255)
                                                                                   aaTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGa         <  1:356013/93‑1 (MQ=255)
                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        <  1:5995118/93‑1 (MQ=255)
                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        <  2:6133244/93‑1 (MQ=255)
                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        >  1:4354396/1‑93 (MQ=255)
                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        <  2:642401/93‑1 (MQ=255)
                                                                                      gCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAcgc      <  2:5291083/93‑1 (MQ=255)
                                                                                          gTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGcgtcg  <  2:3973358/93‑1 (MQ=255)
                                                                                          |                                                                                            
AGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCG  >  NC_000913/1807101‑1807283

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAAGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCGT  >  NC_000913/1807099‑1807284
                                                                                            |                                                                                             
AAAGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAT                                                                                              <  SRR4302547.4111819/94‑1 (MQ=60)
 AAGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATA                                                                                             <  SRR4302547.1213789/94‑1 (MQ=60)
    TCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAA                                                                                          <  SRR4302547.744268/94‑1 (MQ=60)
    TCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAA                                                                                          >  SRR4302547.50722/1‑94 (MQ=60)
     CGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAA                                                                                         >  SRR4302547.2981100/1‑94 (MQ=60)
      TAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAA                                                                                        >  SRR4302547.506821/1‑94 (MQ=60)
       AGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAA                                                                                       <  SRR4302547.1121322/94‑1 (MQ=60)
       AGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAA                                                                                       >  SRR4302547.6654637/1‑94 (MQ=60)
         CATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATC                                                                                     >  SRR4302547.5241290/1‑94 (MQ=60)
          ATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCA                                                                                    >  SRR4302547.5531792/1‑94 (MQ=60)
          ATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCA                                                                                    >  SRR4302547.2296930/1‑94 (MQ=60)
           TGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCAC                                                                                   >  SRR4302547.266601/1‑94 (MQ=60)
              CCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCA                                                                                >  SRR4302547.2599325/1‑94 (MQ=60)
               CGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAA                                                                               <  SRR4302547.1816759/94‑1 (MQ=60)
               CGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAA                                                                               >  SRR4302547.4914541/1‑94 (MQ=60)
               CGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAA                                                                               >  SRR4302547.5960198/1‑94 (MQ=60)
                GCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAAC                                                                              >  SRR4302547.882981/1‑94 (MQ=60)
                GCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAAC                                                                              >  SRR4302547.2552405/1‑94 (MQ=60)
                GCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAAC                                                                              <  SRR4302547.3172811/94‑1 (MQ=60)
                  ATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCC                                                                            >  SRR4302547.1509558/1‑94 (MQ=60)
                    TAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAA                                                                          >  SRR4302547.4331286/1‑94 (MQ=60)
                      ATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAA                                                                        >  SRR4302547.2156357/1‑94 (MQ=60)
                      ATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAA                                                                        >  SRR4302547.2057641/1‑94 (MQ=60)
                         CCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCC                                                                     <  SRR4302547.6098739/94‑1 (MQ=60)
                         CCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCC                                                                     >  SRR4302547.4515164/1‑94 (MQ=60)
                          CTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCT                                                                    >  SRR4302547.1184102/1‑94 (MQ=60)
                          CTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCT                                                                    <  SRR4302547.1456364/94‑1 (MQ=60)
                          CTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCT                                                                    <  SRR4302547.552402/94‑1 (MQ=60)
                           TTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTT                                                                   >  SRR4302547.3191399/1‑94 (MQ=60)
                           TTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTT                                                                   <  SRR4302547.5940922/94‑1 (MQ=60)
                              CTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTC                                                                >  SRR4302547.5493124/1‑94 (MQ=60)
                              CTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTC                                                                >  SRR4302547.6392106/1‑94 (MQ=60)
                               TTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCA                                                               <  SRR4302547.3424650/94‑1 (MQ=60)
                               TTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCA                                                               >  SRR4302547.4196692/1‑94 (MQ=60)
                                   TTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAAC                                                           <  SRR4302547.5577901/94‑1 (MQ=60)
                                     TCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAG                                                         <  SRR4302547.5982579/94‑1 (MQ=60)
                                     TCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAG                                                         >  SRR4302547.1237926/1‑94 (MQ=60)
                                       CTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCA                                                       <  SRR4302547.1569643/94‑1 (MQ=60)
                                       CTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCA                                                       >  SRR4302547.3820845/1‑94 (MQ=60)
                                        TGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAA                                                      >  SRR4302547.5665722/1‑94 (MQ=60)
                                         GCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAAT                                                     >  SRR4302547.5129505/1‑94 (MQ=60)
                                         GCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAAT                                                     >  SRR4302547.2655173/1‑94 (MQ=60)
                                         GCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAAT                                                     >  SRR4302547.5887299/1‑94 (MQ=60)
                                          CCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATG                                                    <  SRR4302547.2730969/94‑1 (MQ=60)
                                           TGGGTATTCATACATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGC                                                   >  SRR4302547.2157436/1‑94 (MQ=60)
                                           CGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGC                                                   >  SRR4302547.5610261/1‑94 (MQ=60)
                                           CGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGC                                                   <  SRR4302547.4978222/94‑1 (MQ=60)
                                           CGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGC                                                   <  SRR4302547.2716982/94‑1 (MQ=60)
                                           CGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGC                                                   >  SRR4302547.2543751/1‑94 (MQ=60)
                                            GGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCG                                                  <  SRR4302547.3372522/94‑1 (MQ=60)
                                            GGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCG                                                  <  SRR4302547.5802436/94‑1 (MQ=60)
                                             GGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGT                                                 <  SRR4302547.322393/94‑1 (MQ=60)
                                              GTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTC                                                <  SRR4302547.1472537/94‑1 (MQ=60)
                                                ATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTGCA                                              >  SRR4302547.3923926/1‑94 (MQ=60)
                                                ATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCA                                              <  SRR4302547.4160315/94‑1 (MQ=60)
                                                ATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCA                                              <  SRR4302547.9377/94‑1 (MQ=60)
                                                  TCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGC                                            >  SRR4302547.1931834/1‑94 (MQ=60)
                                                    ATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTT                                          <  SRR4302547.2284222/94‑1 (MQ=60)
                                                     TAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTG                                         <  SRR4302547.297886/94‑1 (MQ=60)
                                                       AATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTT                                       >  SRR4302547.5359194/1‑94 (MQ=60)
                                                       AATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTT                                       >  SRR4302547.5483749/1‑94 (MQ=60)
                                                       AATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTT                                       >  SRR4302547.4346430/1‑94 (MQ=60)
                                                         acgacgctcttccgatcttgacattATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGA                                     <  SRR4302547.126833/69‑1 (MQ=60)
                                                           TCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATA                                   <  SRR4302547.4143058/94‑1 (MQ=60)
                                                           TCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATA                                   >  SRR4302547.1934576/1‑94 (MQ=60)
                                                           TCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATA                                   <  SRR4302547.3300827/94‑1 (MQ=60)
                                                           TCAGCAAAGAAATCAATGTAGAGATAAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATA                                   >  SRR4302547.2524126/1‑94 (MQ=60)
                                                            CAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATAC                                  <  SRR4302547.5610223/94‑1 (MQ=60)
                                                             AGCAAAGAAGTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACG                                 >  SRR4302547.1931451/1‑94 (MQ=60)
                                                             AGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACG                                 >  SRR4302547.5320955/1‑94 (MQ=60)
                                                              GCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGT                                >  SRR4302547.3173068/1‑94 (MQ=60)
                                                               CAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTACAGCTTGTTGATACGTA                               >  SRR4302547.6743318/1‑94 (MQ=60)
                                                               CAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTC                               >  SRR4302547.4842896/1‑94 (MQ=60)
                                                                AAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCA                              >  SRR4302547.4204580/1‑94 (MQ=60)
                                                                 AAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCAT                             <  SRR4302547.4108519/94‑1 (MQ=60)
                                                                 AAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCAT                             <  SRR4302547.5198806/94‑1 (MQ=60)
                                                                  AGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATG                            >  SRR4302547.5530466/1‑94 (MQ=60)
                                                                   GAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGA                           >  SRR4302547.5816243/1‑94 (MQ=60)
                                                                    AAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAG                          <  SRR4302547.5831596/94‑1 (MQ=60)
                                                                       TCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTC                       >  SRR4302547.1164464/1‑94 (MQ=60)
                                                                       TCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTC                       >  SRR4302547.5413905/1‑94 (MQ=60)
                                                                        CAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCT                      >  SRR4302547.1693044/1‑94 (MQ=60)
                                                                        CAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCT                      >  SRR4302547.6011332/1‑94 (MQ=60)
                                                                        CAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCT                      >  SRR4302547.1869431/1‑94 (MQ=60)
                                                                        CAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCT                      >  SRR4302547.4412219/1‑94 (MQ=60)
                                                                        CAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCT                      <  SRR4302547.5642429/94‑1 (MQ=60)
                                                                          ATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCC                    >  SRR4302547.3718157/1‑94 (MQ=60)
                                                                          ATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCC                    >  SRR4302547.3742830/1‑94 (MQ=60)
                                                                            GTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTT                  <  SRR4302547.102287/94‑1 (MQ=60)
                                                                            GTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTT                  <  SRR4302547.4059311/94‑1 (MQ=60)
                                                                            GTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTT                  <  SRR4302547.5380800/94‑1 (MQ=60)
                                                                            GTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTT                  <  SRR4302547.3400035/94‑1 (MQ=60)
                                                                             TAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTG                 >  SRR4302547.5876507/1‑94 (MQ=60)
                                                                              AGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGA                >  SRR4302547.1139417/1‑94 (MQ=60)
                                                                              AGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGA                <  SRR4302547.2459384/94‑1 (MQ=60)
                                                                               GAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGAT               <  SRR4302547.552798/94‑1 (MQ=60)
                                                                                  ATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAA            >  SRR4302547.654289/1‑94 (MQ=60)
                                                                                  ATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAA            >  SRR4302547.1414657/1‑94 (MQ=60)
                                                                                  ATCAATGCAGATATAAAAATAACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAA            >  SRR4302547.6599318/1‑94 (MQ=60)
                                                                                   TCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAG           >  SRR4302547.5835233/1‑94 (MQ=60)
                                                                                    CAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGA          >  SRR4302547.2024262/1‑94 (MQ=60)
                                                                                    CAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGA          <  SRR4302547.5363233/94‑1 (MQ=60)
                                                                                    CAATGCAGATATAAAAATCACACAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTcagt          >  SRR4302547.1911203/1‑90 (MQ=60)
                                                                                     AATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAC         <  SRR4302547.362616/94‑1 (MQ=60)
                                                                                     AATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAC         >  SRR4302547.4443849/1‑94 (MQ=60)
                                                                                      ATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACG        <  SRR4302547.6119524/94‑1 (MQ=60)
                                                                                          AGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTC    >  SRR4302547.4920944/1‑94 (MQ=60)
                                                                                           GATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCG   >  SRR4302547.3885088/1‑94 (MQ=60)
                                                                                           GATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCG   <  SRR4302547.5703977/94‑1 (MQ=60)
                                                                                           GATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCG   >  SRR4302547.3122880/1‑94 (MQ=60)
                                                                                            GTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCGT  >  SRR4302547.2541747/1‑94 (MQ=60)
                                                                                            ATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCGT  >  SRR4302547.6251025/1‑94 (MQ=60)
                                                                                            |                                                                                             
AAAGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCGT  >  NC_000913/1807099‑1807284

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: