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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,997,095 C→T A137T (GCC→ACC)  b1919 ← putative 1‑aminocyclopropane‑1‑carboxylate deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,997,0950CT94.7% 464.1 / 6.8 153A137T (GCC→ACC) b1919putative 1‑aminocyclopropane‑1‑carboxylate deaminase
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base T (69/75);  total (72/81)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

Only 100 of 158 total aligned reads displayed.

GCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTG  >  NC_000913/1997005‑1997184
                                                                                          |                                                                                         
 cgcgcCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCAtt                                                                                        <  2:3555939/93‑1 (MQ=255)
  gcgcCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTg                                                                                       <  1:133695/93‑1 (MQ=255)
    gcCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTggg                                                                                     >  2:1356501/1‑93 (MQ=255)
     ccAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGa                                                                                    >  2:2392485/1‑93 (MQ=255)
      cAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGAt                                                                                   <  1:3863159/93‑1 (MQ=255)
      cAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGAt                                                                                   >  1:912147/1‑93 (MQ=255)
      cAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGAt                                                                                   >  2:1282185/1‑93 (MQ=255)
      cAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGAt                                                                                   >  2:2232872/1‑93 (MQ=255)
       aaGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATc                                                                                  >  1:4472478/1‑93 (MQ=255)
       aaGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATc                                                                                  >  2:2327656/1‑93 (MQ=255)
       aaGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATc                                                                                  >  2:4294972/1‑93 (MQ=255)
        agagCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCg                                                                                 <  2:267543/93‑1 (MQ=255)
           gCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTc                                                                              <  2:1637176/93‑1 (MQ=255)
           gCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTc                                                                              <  1:881011/93‑1 (MQ=255)
            cATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCa                                                                             >  1:2458814/1‑93 (MQ=255)
               tAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTg                                                                          >  1:1934200/1‑93 (MQ=255)
                aGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGc                                                                         <  1:778783/93‑1 (MQ=255)
                aGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGc                                                                         <  1:2556082/93‑1 (MQ=255)
                   aCCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCCGTGCGTc                                                                      >  2:3651501/1‑93 (MQ=255)
                    ccGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGCGCGTCg                                                                     >  2:221074/1‑93 (MQ=255)
                     cGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCCGTGCGgcgc                                                                    >  2:2752319/1‑93 (MQ=255)
                     cGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGc                                                                    >  2:4096149/1‑93 (MQ=255)
                       ccAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGcac                                                                  >  1:3070156/1‑93 (MQ=255)
                        cAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGcaca                                                                 >  1:4056234/1‑93 (MQ=255)
                         aaCCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACAt                                                                >  1:1118092/1‑93 (MQ=255)
                         aaCCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACAt                                                                <  1:1152859/93‑1 (MQ=255)
                           ccGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACAttt                                                              >  1:3212513/1‑93 (MQ=255)
                           ccGGAATGACATACGGGCGAAAGCCTTGTGCTTCCACTCGCGTCGCCAGGTCTTCCAGTTGGGTATTGGGATCAGTCCGTGCGTCGCACAttt                                                              >  2:3492837/1‑93 (MQ=255)
                            cGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTc                                                             <  1:2293935/93‑1 (MQ=255)
                            cGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTc                                                             <  1:1276679/93‑1 (MQ=255)
                             ggAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCa                                                            <  1:598780/93‑1 (MQ=255)
                              gAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCaa                                                           <  2:2968991/93‑1 (MQ=255)
                                aTGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATc                                                         >  2:3485648/1‑93 (MQ=255)
                                   aCATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTgg                                                      >  1:3747523/1‑93 (MQ=255)
                                    cATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTggg                                                     >  1:4351288/1‑93 (MQ=255)
                                     aTACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGt                                                    >  2:3020316/1‑93 (MQ=255)
                                      tACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGtt                                                   >  2:3263803/1‑92 (MQ=255)
                                      tACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTa                                                   >  1:3847574/1‑93 (MQ=255)
                                      tACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTa                                                   >  1:3422141/1‑93 (MQ=255)
                                      tACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTa                                                   >  2:3428977/1‑93 (MQ=255)
                                       aCGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTAt                                                  <  2:4535927/93‑1 (MQ=255)
                                        cGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTAtt                                                 <  1:2604518/93‑1 (MQ=255)
                                         gggCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTg                                                >  1:3534364/1‑93 (MQ=255)
                                            ggAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAAc                                             <  2:3846712/92‑1 (MQ=255)
                                            cGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAAc                                             <  2:286931/93‑1 (MQ=255)
                                              aaaGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAg                                           >  1:2563161/1‑93 (MQ=255)
                                              aaaGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAg                                           >  1:1188328/1‑93 (MQ=255)
                                               aaGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGa                                          >  2:3508511/1‑93 (MQ=255)
                                               aaGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGa                                          <  2:3240335/93‑1 (MQ=255)
                                               aaGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGa                                          <  2:2601638/93‑1 (MQ=255)
                                                aGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGAt                                         <  2:4273963/93‑1 (MQ=255)
                                                aGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGAt                                         <  2:3406822/93‑1 (MQ=255)
                                                 gCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGGCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATc                                        <  2:1303069/93‑1 (MQ=255)
                                                   cTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCa                                      >  2:2508746/1‑93 (MQ=255)
                                                   cTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCa                                      >  1:1405316/1‑93 (MQ=255)
                                                    ttGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAg                                     >  1:3207643/1‑93 (MQ=255)
                                                     tgtgCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGc                                    <  1:3334985/93‑1 (MQ=255)
                                                        gCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGcaac                                 <  2:3836570/93‑1 (MQ=255)
                                                          ttCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGcaacaa                               >  2:3160954/1‑93 (MQ=255)
                                                           tCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAg                              <  1:2918977/93‑1 (MQ=255)
                                                           tCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAg                              >  1:1769495/1‑93 (MQ=255)
                                                           tCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAg                              <  1:1580256/93‑1 (MQ=255)
                                                            cGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGc                             <  1:4623497/93‑1 (MQ=255)
                                                            cGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGc                             <  1:2390183/93‑1 (MQ=255)
                                                              aCTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGa                           >  2:3885007/1‑93 (MQ=255)
                                                              aCTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGa                           >  2:4638947/1‑93 (MQ=255)
                                                              aCTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGa                           >  2:2992103/1‑93 (MQ=255)
                                                              aCTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGAATTGAACAGATCCAGCAACAAGCGa                           >  1:1128980/1‑93 (MQ=255)
                                                               caagcgTCTCCAGCTCTTCCAGTTGGGAATGGGGATCGGTCATTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGAt                          <  1:1279222/90‑1 (MQ=255)
                                                               cTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGAt                          <  2:1506249/93‑1 (MQ=255)
                                                               cTCGCGTCGCCAGCTCTCCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGAt                          <  1:2807664/93‑1 (MQ=255)
                                                                tCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGAtt                         <  2:3675982/93‑1 (MQ=255)
                                                                tCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGAtt                         <  2:790285/93‑1 (MQ=255)
                                                                tCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGAtt                         <  1:2198329/93‑1 (MQ=255)
                                                                tCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGAtt                         <  2:3528589/93‑1 (MQ=255)
                                                                    gTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCg                     >  1:956723/1‑93 (MQ=255)
                                                                      cGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGtt                   <  1:1680381/93‑1 (MQ=255)
                                                                      cGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGtt                   >  2:2188159/1‑93 (MQ=255)
                                                                        ccAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTgg                 >  1:3827502/1‑93 (MQ=255)
                                                                        ccAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTgg                 >  1:2289178/1‑93 (MQ=255)
                                                                        ccAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTgg                 >  2:4000197/1‑93 (MQ=255)
                                                                         cAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGt                >  1:4489207/1‑93 (MQ=255)
                                                                         cAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGt                <  1:2201574/93‑1 (MQ=255)
                                                                           gCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTa              >  1:3985761/1‑93 (MQ=255)
                                                                           gCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACAGTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTa              >  1:3295216/1‑93 (MQ=255)
                                                                            ctctTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTaa             >  1:2839639/1‑93 (MQ=255)
                                                                              ctTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAAt           <  1:739529/93‑1 (MQ=255)
                                                                               ttCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATa          <  2:709717/93‑1 (MQ=255)
                                                                               ttCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATa          <  2:1029950/93‑1 (MQ=255)
                                                                               ttCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATa          >  2:3968594/1‑93 (MQ=255)
                                                                               ttCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATa          <  2:4191025/93‑1 (MQ=255)
                                                                                   aGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGttt      <  2:1145032/93‑1 (MQ=255)
                                                                                   aGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGttt      <  1:3036489/93‑1 (MQ=255)
                                                                                   aGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGttt      >  2:661767/1‑93 (MQ=255)
                                                                                   aGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGttt      <  2:3340133/93‑1 (MQ=255)
                                                                                    gTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGtttt     >  2:4355050/1‑93 (MQ=255)
                                                                                     ttGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTc    <  1:4587283/93‑1 (MQ=255)
                                                                                     ttGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTc    <  2:1429337/93‑1 (MQ=255)
                                                                                      tGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCt   <  2:1667952/93‑1 (MQ=255)
                                                                                       gggTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTg  <  1:4380024/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTG  >  NC_000913/1997005‑1997184

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTAGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGT  >  NC_000913/1997002‑1997188
                                                                                             |                                                                                             
CTAGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGT                                                                                               >  SRR4302535.4655401/1‑94 (MQ=60)
 TAGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTA                                                                                              <  SRR4302535.4506759/94‑1 (MQ=60)
  AGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTAT                                                                                             <  SRR4302535.4365425/94‑1 (MQ=60)
  AGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTAT                                                                                             <  SRR4302535.3285799/94‑1 (MQ=60)
   GCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATT                                                                                            <  SRR4302535.1313956/94‑1 (MQ=60)
    CGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTG                                                                                           <  SRR4302535.2214161/94‑1 (MQ=60)
    CGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTG                                                                                           <  SRR4302535.1211162/94‑1 (MQ=60)
     GCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGG                                                                                          <  SRR4302535.342259/94‑1 (MQ=60)
     GCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGG                                                                                          <  SRR4302535.134630/94‑1 (MQ=60)
        CCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGAT                                                                                       <  SRR4302535.3784029/94‑1 (MQ=60)
        CCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGAT                                                                                       >  SRR4302535.918634/1‑94 (MQ=60)
         CAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATC                                                                                      >  SRR4302535.4503173/1‑94 (MQ=60)
         CAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATC                                                                                      <  SRR4302535.3889735/94‑1 (MQ=60)
              GCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCA                                                                                 >  SRR4302535.3175634/1‑94 (MQ=60)
              GCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCA                                                                                 <  SRR4302535.887248/94‑1 (MQ=60)
              GCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCA                                                                                 >  SRR4302535.2476546/1‑94 (MQ=60)
               CATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAG                                                                                >  SRR4302535.2508221/1‑94 (MQ=60)
                ATTAGAACCGGCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGT                                                                               <  SRR4302535.1159100/94‑1 (MQ=60)
                 TTAGAACCGCCACCCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTG                                                                              <  SRR4302535.1903104/94‑1 (MQ=60)
                 TTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTG                                                                              >  SRR4302535.1948130/1‑94 (MQ=60)
                 TTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTG                                                                              >  SRR4302535.3656450/1‑94 (MQ=60)
                   AGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCG                                                                            <  SRR4302535.784254/94‑1 (MQ=60)
                   AGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCG                                                                            <  SRR4302535.2574489/94‑1 (MQ=60)
                      ACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCG                                                                         >  SRR4302535.4309392/1‑94 (MQ=60)
                      ACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCG                                                                         <  SRR4302535.1370487/94‑1 (MQ=60)
                         GCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCAC                                                                      >  SRR4302535.3091598/1‑94 (MQ=60)
                          CCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACA                                                                     >  SRR4302535.4084004/1‑94 (MQ=60)
                           CAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACAT                                                                    >  SRR4302535.1126068/1‑94 (MQ=60)
                           CAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACAT                                                                    <  SRR4302535.875550/94‑1 (MQ=60)
                            AACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATT                                                                   <  SRR4302535.1161098/94‑1 (MQ=60)
                             ACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTT                                                                  >  SRR4302535.3234761/1‑94 (MQ=60)
                             ACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTT                                                                  >  SRR4302535.2920888/1‑94 (MQ=60)
                               CGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCA                                                                <  SRR4302535.2310342/94‑1 (MQ=60)
                               CGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCA                                                                <  SRR4302535.1285794/94‑1 (MQ=60)
                                GGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAA                                                               <  SRR4302535.602966/94‑1 (MQ=60)
                                  AATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATC                                                             >  SRR4302535.540443/1‑94 (MQ=60)
                                  AATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGGATTGGGATCGGTCAGTGCGTCGCACATTTCAATC                                                             >  SRR4302535.2233027/1‑94 (MQ=60)
                                   ATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCT                                                            <  SRR4302535.3884348/94‑1 (MQ=60)
                                     GACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGG                                                          >  SRR4302535.3773449/1‑94 (MQ=60)
                                      ACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGG                                                         >  SRR4302535.4381061/1‑94 (MQ=60)
                                        ATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTA                                                       >  SRR4302535.3874059/1‑94 (MQ=60)
                                        ATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTA                                                       >  SRR4302535.3445945/1‑94 (MQ=60)
                                        ATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTA                                                       >  SRR4302535.2711090/1‑94 (MQ=60)
                                           CGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTG                                                    >  SRR4302535.3558902/1‑94 (MQ=60)
                                           CGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTG                                                    >  SRR4302535.72569/1‑94 (MQ=60)
                                           CGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTG                                                    <  SRR4302535.2623214/94‑1 (MQ=60)
                                               CGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACA                                                <  SRR4302535.979114/94‑1 (MQ=60)
                                                GAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAG                                               >  SRR4302535.1196838/1‑94 (MQ=60)
                                                GAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAG                                               >  SRR4302535.2581608/1‑94 (MQ=60)
                                                   AGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATC                                            >  SRR4302535.3082975/1‑94 (MQ=60)
                                                    GCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCC                                           <  SRR4302535.3308622/94‑1 (MQ=60)
                                                     CCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCA                                          <  SRR4302535.953549/94‑1 (MQ=60)
                                                     CCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCA                                          >  SRR4302535.1415361/1‑94 (MQ=60)
                                                      CTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAG                                         >  SRR4302535.3229861/1‑94 (MQ=60)
                                                        TGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCA                                       <  SRR4302535.3358211/94‑1 (MQ=60)
                                                             TTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAG                                  >  SRR4302535.1782323/1‑94 (MQ=60)
                                                              TCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGC                                 <  SRR4302535.2939576/94‑1 (MQ=60)
                                                              TCGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGC                                 <  SRR4302535.1591663/94‑1 (MQ=60)
                                                               CGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCG                                <  SRR4302535.4654898/94‑1 (MQ=60)
                                                               CGACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCG                                <  SRR4302535.2407369/94‑1 (MQ=60)
                                                                GACTCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGAATTGAACAGATCCAGCAACAAGCGA                               >  SRR4302535.1137041/1‑94 (MQ=60)
                                                                  caaGCGTCTCCAGCTCTTCCAGTTGGGAATGGGGATCGGTCATTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATT                             <  SRR4302535.1288356/91‑1 (MQ=60)
                                                                  CTCGCGTCGCCAGCTCTCCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATT                             <  SRR4302535.2827621/94‑1 (MQ=60)
                                                                   TCGCGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTG                            <  SRR4302535.2213992/94‑1 (MQ=60)
                                                                   TCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTG                            <  SRR4302535.2326869/94‑1 (MQ=60)
                                                                      CGTCGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCG                         >  SRR4302535.963513/1‑94 (MQ=60)
                                                                         CGCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTG                      <  SRR4302535.1692552/94‑1 (MQ=60)
                                                                          GGCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGG                     <  SRR4302535.257771/94‑1 (MQ=60)
                                                                          GCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGG                     >  SRR4302535.2305529/1‑94 (MQ=60)
                                                                          GCCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGG                     >  SRR4302535.3853905/1‑94 (MQ=60)
                                                                           CCAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGT                    >  SRR4302535.4520024/1‑94 (MQ=60)
                                                                            CAGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTT                   <  SRR4302535.2217265/94‑1 (MQ=60)
                                                                             AGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTA                  >  SRR4302535.4013000/1‑94 (MQ=60)
                                                                             AGCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACAGTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTA                  >  SRR4302535.3318126/1‑94 (MQ=60)
                                                                              GCTCTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAA                 >  SRR4302535.2859782/1‑94 (MQ=60)
                                                                                 CTTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATA              <  SRR4302535.744727/94‑1 (MQ=60)
                                                                                  TTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAG             <  SRR4302535.1651042/94‑1 (MQ=60)
                                                                                  TTCCAGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAG             <  SRR4302535.4072967/94‑1 (MQ=60)
                                                                                   TCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGT            <  SRR4302535.4424903/94‑1 (MQ=60)
                                                                                      AGTTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTT         <  SRR4302535.3057757/94‑1 (MQ=60)
                                                                                        TTGGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCT       <  SRR4302535.4618520/94‑1 (MQ=60)
                                                                                        GTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCT       <  SRR4302535.4570867/94‑1 (MQ=60)
                                                                                          GGGTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGC     <  SRR4302535.4409991/94‑1 (MQ=60)
                                                                                            GTATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGG   >  SRR4302535.2556930/1‑94 (MQ=60)
                                                                                             TATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGT  <  SRR4302535.1653415/94‑1 (MQ=60)
                                                                                             |                                                                                             
CTAGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGT  >  NC_000913/1997002‑1997188

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