Sample Resequencing Stats

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

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,965,062 A→T H344L (CAC→CTC)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,0620AT96.5% 866.1 / 17.4 255H344L (CAC→CTC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (3/6);  new base T (139/107);  total (142/113)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 260 total aligned reads displayed.

ACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTC  >  NC_000913/3964971‑3965151
                                                                                           |                                                                                         
aCCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCAc                                                                                          <  2:1853049/93‑1 (MQ=255)
  cATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCct                                                                                        >  1:2058455/1‑93 (MQ=255)
   aTTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGctcctc                                                                                       >  2:3439712/1‑93 (MQ=255)
     tATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCctctc                                                                                     >  2:2104657/1‑93 (MQ=255)
       tCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCtctctc                                                                                   <  1:4060318/93‑1 (MQ=255)
         gcgACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGt                                                                                 >  1:3151932/1‑93 (MQ=255)
         gcgACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGt                                                                                 <  2:3186099/93‑1 (MQ=255)
          cgACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTa                                                                                >  2:6455826/1‑93 (MQ=255)
          cgACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTa                                                                                <  1:287797/93‑1 (MQ=255)
          cgACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTa                                                                                <  1:2654972/93‑1 (MQ=255)
           gACGGCGCTTATCGATACCGGTTCTAAAATGGACGGAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTaa                                                                               >  2:4683587/1‑93 (MQ=255)
           gACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTCCAGGCAACATGGAACTGCTCCTCTCTCGTaa                                                                               <  2:5661098/93‑1 (MQ=255)
            aCGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAg                                                                              >  2:1520347/1‑93 (MQ=255)
            aCGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAg                                                                              >  1:3187565/1‑93 (MQ=255)
               gcgcTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATc                                                                           >  1:6191/1‑93 (MQ=255)
               gcgcTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATc                                                                           >  2:3924884/1‑93 (MQ=255)
                   ttATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTg                                                                       >  1:1097463/1‑93 (MQ=255)
                    tATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGa                                                                      >  2:3295003/1‑93 (MQ=255)
                    tATCCATACCGGTGCTAAAATGGACGAAGTTATCTACGGAGATTTGAAAGGGACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGa                                                                      <  2:58045/93‑1 (MQ=255)
                     aTCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGaa                                                                     <  1:3600317/93‑1 (MQ=255)
                      tCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATAGCTGaaa                                                                    >  2:1876019/1‑93 (MQ=255)
                       cGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGaaaa                                                                   <  1:3439615/93‑1 (MQ=255)
                        gATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGaaaaa                                                                  >  1:5006513/1‑93 (MQ=255)
                         aTACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAAc                                                                 <  1:972744/93‑1 (MQ=255)
                          tACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAAcg                                                                >  2:2991844/1‑93 (MQ=255)
                          tACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAAcg                                                                >  1:1441281/1‑93 (MQ=255)
                          tACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAAcg                                                                >  2:1013159/1‑93 (MQ=255)
                            ccGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAAcgcg                                                              <  1:4467767/93‑1 (MQ=255)
                             cGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGt                                                             <  2:1406252/93‑1 (MQ=255)
                             cGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGt                                                             >  1:610778/1‑93 (MQ=255)
                              ggTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTc                                                            >  1:3956971/1‑93 (MQ=255)
                              ggTGCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTc                                                            <  1:3489117/93‑1 (MQ=255)
                               gTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCt                                                           <  1:305562/93‑1 (MQ=255)
                               gTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCt                                                           >  1:4357454/1‑93 (MQ=255)
                               gTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCt                                                           <  2:715419/93‑1 (MQ=255)
                               gTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCt                                                           >  1:5889656/1‑93 (MQ=255)
                                ttCTAAAATGGACGAAGTTATCTACGAAGAGTTTACAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCtt                                                          >  1:694215/1‑93 (MQ=255)
                                ttCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCtt                                                          <  2:3256225/93‑1 (MQ=255)
                                ttCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCtt                                                          >  2:957038/1‑93 (MQ=255)
                                ttCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCtt                                                          <  1:426413/93‑1 (MQ=255)
                                ttCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCtt                                                          <  1:5969834/93‑1 (MQ=255)
                                 tCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTc                                                         >  1:562993/1‑93 (MQ=255)
                                 tCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTc                                                         <  1:1755342/93‑1 (MQ=255)
                                 tCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTc                                                         >  2:4326013/1‑93 (MQ=255)
                                 tCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTc                                                         >  1:6353047/1‑93 (MQ=255)
                                  cTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTcc                                                        >  1:5347321/1‑93 (MQ=255)
                                   tAAAATGTACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGTAACTGCTCCTCTCTCGTAATATCGCTGAAAAACGCGTCTTccc                                                       <  2:612748/93‑1 (MQ=255)
                                   tAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTccc                                                       <  1:2262336/93‑1 (MQ=255)
                                   tAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTccc                                                       <  2:5345467/93‑1 (MQ=255)
                                   tAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTccc                                                       >  1:1913703/1‑93 (MQ=255)
                                   tAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTccc                                                       <  1:5469234/93‑1 (MQ=255)
                                    aaaaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCg                                                      <  2:5305000/93‑1 (MQ=255)
                                     aaaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCgg                                                     >  1:1598540/1‑93 (MQ=255)
                                     aaaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCgc                                                     >  1:3745884/1‑92 (MQ=255)
                                      aaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGc                                                    >  2:234218/1‑93 (MQ=255)
                                       aTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCt                                                   >  2:5525648/1‑93 (MQ=255)
                                       aTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCt                                                   >  2:3516505/1‑93 (MQ=255)
                                               aGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTa                                           >  1:4900611/1‑93 (MQ=255)
                                                  tATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACaa                                        <  1:1032618/93‑1 (MQ=255)
                                                  tATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACaa                                        <  1:1566091/93‑1 (MQ=255)
                                                  tATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACaa                                        >  2:1578521/1‑93 (MQ=255)
                                                    tCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                      <  2:5725518/93‑1 (MQ=255)
                                                    tCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                      >  2:5777170/1‑93 (MQ=255)
                                                    tCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                      >  1:220834/1‑93 (MQ=255)
                                                     cTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCg                                     <  1:1855994/93‑1 (MQ=255)
                                                      tACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGt                                    >  1:1241684/1‑93 (MQ=255)
                                                       aCGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGtt                                   >  2:3986897/1‑93 (MQ=255)
                                                       aCGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGtt                                   <  2:4610863/93‑1 (MQ=255)
                                                          aaGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                >  2:1227940/1‑93 (MQ=255)
                                                          aaGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                >  2:4397437/1‑93 (MQ=255)
                                                          aaGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                >  1:4767140/1‑93 (MQ=255)
                                                          aaGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                >  2:2549348/1‑93 (MQ=255)
                                                            gagTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGt                              >  1:1179673/1‑93 (MQ=255)
                                                             agTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTa                             <  1:3654113/93‑1 (MQ=255)
                                                               tttAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAcc                           <  1:242699/93‑1 (MQ=255)
                                                                ttAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAccc                          <  1:4289244/93‑1 (MQ=255)
                                                                  aaaGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGt                        >  2:3852126/1‑93 (MQ=255)
                                                                     ggTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTaaa                     >  2:2976663/1‑93 (MQ=255)
                                                                      gTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAaag                    >  2:224580/1‑93 (MQ=255)
                                                                        aCAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAaagaa                  <  2:2260779/93‑1 (MQ=255)
                                                                         cAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAaagaag                 >  2:2749715/1‑93 (MQ=255)
                                                                           ggCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAagag               >  1:6460139/1‑93 (MQ=255)
                                                                           ggCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAagag               >  2:4261310/1‑93 (MQ=255)
                                                                            gCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgc              <  2:3652054/93‑1 (MQ=255)
                                                                             cAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgct             >  2:5916769/1‑93 (MQ=255)
                                                                             cAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgct             <  1:559040/93‑1 (MQ=255)
                                                                              aaCATGGAACTGCTCCTCTCTNGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctg            <  1:2567323/93‑1 (MQ=255)
                                                                               aCATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgc           >  1:2778513/1‑93 (MQ=255)
                                                                                 aTGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTc         >  1:3366732/1‑93 (MQ=255)
                                                                                  tGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCa        <  1:4618779/93‑1 (MQ=255)
                                                                                  tGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCa        <  2:6353047/93‑1 (MQ=255)
                                                                                    gAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACg      >  1:1804872/1‑93 (MQ=255)
                                                                                     aaCTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa     >  2:128714/1‑93 (MQ=255)
                                                                                      aCTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGAc    >  2:4772629/1‑93 (MQ=255)
                                                                                      aCTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCCTCTTCCCGCCTATCGACTACAACCGTCCTGGTACCCGTAAAGAAGAGCTGCTCACGAc    <  1:5764284/93‑1 (MQ=255)
                                                                                       cTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACACCCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACt   <  2:3342188/93‑1 (MQ=255)
                                                                                       cTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACt   >  1:4417962/1‑93 (MQ=255)
                                                                                       cTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCAGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACt   >  1:508948/1‑93 (MQ=255)
                                                                                        tGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTc  <  2:4248806/93‑1 (MQ=255)
                                                                                        tGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTc  <  2:3366715/93‑1 (MQ=255)
                                                                                           |                                                                                         
ACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTC  >  NC_000913/3964971‑3965151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ 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
GACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGA  >  NC_000913/3964970‑3965155
                                                                                            |                                                                                             
GACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTC                                                                                              >  SRR4302537.1392487/1‑94 (MQ=60)
  CCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCT                                                                                            >  SRR4302537.2072907/1‑94 (MQ=60)
  CCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAATTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCT                                                                                            <  SRR4302537.4222613/94‑1 (MQ=60)
     TTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTC                                                                                         >  SRR4302537.196008/1‑94 (MQ=60)
     TTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCCACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTC                                                                                         <  SRR4302537.3482362/94‑1 (MQ=60)
      TATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCT                                                                                        <  SRR4302537.487540/94‑1 (MQ=60)
      TATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCT                                                                                        <  SRR4302537.4741163/94‑1 (MQ=60)
       ATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTC                                                                                       >  SRR4302537.6070551/1‑94 (MQ=60)
       ATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTC                                                                                       >  SRR4302537.2092682/1‑94 (MQ=60)
       ATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTC                                                                                       <  SRR4302537.4460784/94‑1 (MQ=60)
        TCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCG                                                                                      <  SRR4302537.4088694/94‑1 (MQ=60)
        TCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCG                                                                                      <  SRR4302537.5457725/94‑1 (MQ=60)
         CGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGT                                                                                     >  SRR4302537.3174065/1‑94 (MQ=60)
         CGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGT                                                                                     <  SRR4302537.6421645/94‑1 (MQ=60)
          GCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTA                                                                                    >  SRR4302537.128999/1‑94 (MQ=60)
          GCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTA                                                                                    <  SRR4302537.5825179/94‑1 (MQ=60)
           CGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAA                                                                                   >  SRR4302537.1469383/1‑94 (MQ=60)
           CGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAA                                                                                   <  SRR4302537.3541915/94‑1 (MQ=60)
           CGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAA                                                                                   <  SRR4302537.289797/94‑1 (MQ=60)
           CGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAA                                                                                   <  SRR4302537.2673531/94‑1 (MQ=60)
            GACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAG                                                                                  >  SRR4302537.3209987/1‑94 (MQ=60)
               GGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATC                                                                               >  SRR4302537.3673792/1‑94 (MQ=60)
               GGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATC                                                                               >  SRR4302537.6237/1‑94 (MQ=60)
                GCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCG                                                                              >  SRR4302537.1601134/1‑94 (MQ=60)
                GCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCG                                                                              <  SRR4302537.1870923/94‑1 (MQ=60)
                 GGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGC                                                                             <  SRR4302537.1216452/94‑1 (MQ=60)
                 CGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGC                                                                             >  SRR4302537.5748033/1‑94 (MQ=60)
                 CGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGC                                                                             <  SRR4302537.4035223/94‑1 (MQ=60)
                 CGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTATCTCGTAAGATCGC                                                                             >  SRR4302537.2787034/1‑94 (MQ=60)
                 CGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGGTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGC                                                                             >  SRR4302537.5830800/1‑94 (MQ=60)
                   CTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTG                                                                           >  SRR4302537.514671/1‑94 (MQ=60)
                   CTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTG                                                                           >  SRR4302537.1105046/1‑94 (MQ=60)
                   CTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTG                                                                           >  SRR4302537.454182/1‑94 (MQ=60)
                    TTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAATGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGA                                                                          <  SRR4302537.3481116/94‑1 (MQ=60)
                    TTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGA                                                                          >  SRR4302537.1674911/1‑94 (MQ=60)
                    TTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGA                                                                          >  SRR4302537.3190854/1‑94 (MQ=60)
                     TATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAA                                                                         >  SRR4302537.1831548/1‑94 (MQ=60)
                     TATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAA                                                                         >  SRR4302537.4205472/1‑94 (MQ=60)
                      ATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAA                                                                        <  SRR4302537.1795193/94‑1 (MQ=60)
                      ATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAA                                                                        <  SRR4302537.3625692/94‑1 (MQ=60)
                      ATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAA                                                                        <  SRR4302537.4326184/94‑1 (MQ=60)
                       TCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAA                                                                       <  SRR4302537.233684/94‑1 (MQ=60)
                        CGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAA                                                                      >  SRR4302537.3241500/1‑94 (MQ=60)
                        CGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAA                                                                      <  SRR4302537.3463958/94‑1 (MQ=60)
                        CGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAA                                                                      >  SRR4302537.5040994/1‑94 (MQ=60)
                          ATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACG                                                                    >  SRR4302537.2916395/1‑94 (MQ=60)
                          ATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACG                                                                    <  SRR4302537.979449/94‑1 (MQ=60)
                          ATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACG                                                                    <  SRR4302537.6360582/94‑1 (MQ=60)
                          ATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACG                                                                    >  SRR4302537.1451327/1‑94 (MQ=60)
                          ATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACG                                                                    <  SRR4302537.5355329/94‑1 (MQ=60)
                            ACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCG                                                                  >  SRR4302537.4307977/1‑94 (MQ=60)
                             CCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGT                                                                 <  SRR4302537.2799284/94‑1 (MQ=60)
                             CCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGT                                                                 <  SRR4302537.4498665/94‑1 (MQ=60)
                             CCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGT                                                                 >  SRR4302537.614950/1‑94 (MQ=60)
                              CGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAACAACGCGTC                                                                >  SRR4302537.3038864/1‑94 (MQ=60)
                              CGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTC                                                                >  SRR4302537.3984657/1‑94 (MQ=60)
                               GGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCT                                                               >  SRR4302537.4387616/1‑94 (MQ=60)
                               GGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCT                                                               >  SRR4302537.5930065/1‑94 (MQ=60)
                               GGTGCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCT                                                               <  SRR4302537.3513797/94‑1 (MQ=60)
                                GTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTACAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTT                                                              >  SRR4302537.698974/1‑94 (MQ=60)
                                GTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTT                                                              <  SRR4302537.307681/94‑1 (MQ=60)
                                 TTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTC                                                             >  SRR4302537.566867/1‑94 (MQ=60)
                                 TTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTC                                                             <  SRR4302537.2453568/94‑1 (MQ=60)
                                 TTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTC                                                             >  SRR4302537.6396418/1‑94 (MQ=60)
                                 TTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTC                                                             <  SRR4302537.429374/94‑1 (MQ=60)
                                 TTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTC                                                             <  SRR4302537.6010810/94‑1 (MQ=60)
                                  TCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCC                                                            >  SRR4302537.5383959/1‑94 (MQ=60)
                                  TCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCC                                                            >  SRR4302537.1056634/1‑94 (MQ=60)
                                  TCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCC                                                            <  SRR4302537.266283/94‑1 (MQ=60)
                                  TCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCC                                                            <  SRR4302537.1767576/94‑1 (MQ=60)
                                  TCTAAAATGGACGAAGTTATCTACGAAGAGGTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCC                                                            >  SRR4302537.6001220/1‑94 (MQ=60)
                                   CTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCC                                                           >  SRR4302537.2435878/1‑94 (MQ=60)
                                   CTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCC                                                           >  SRR4302537.1926985/1‑94 (MQ=60)
                                   CTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCC                                                           <  SRR4302537.5599900/94‑1 (MQ=60)
                                    TAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCG                                                          <  SRR4302537.2278301/94‑1 (MQ=60)
                                    TAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCG                                                          <  SRR4302537.5506541/94‑1 (MQ=60)
                                     AAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGG                                                         >  SRR4302537.1609710/1‑94 (MQ=60)
                                     AAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGC                                                         >  SRR4302537.3772190/1‑94 (MQ=60)
                                        ATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTA                                                      >  SRR4302537.4476988/1‑94 (MQ=60)
                                           GACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCG                                                   <  SRR4302537.4942930/94‑1 (MQ=60)
                                            ACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGA                                                  <  SRR4302537.4582009/94‑1 (MQ=60)
                                             CGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGAC                                                 <  SRR4302537.2154298/94‑1 (MQ=60)
                                              GAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACT                                                <  SRR4302537.6470526/94‑1 (MQ=60)
                                              GAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACA                                                >  SRR4302537.4057170/1‑94 (MQ=60)
                                               AAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTA                                               >  SRR4302537.5529908/1‑94 (MQ=60)
                                               AAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTA                                               >  SRR4302537.4934441/1‑94 (MQ=60)
                                                  TTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAA                                            <  SRR4302537.5098702/94‑1 (MQ=60)
                                                   TATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAC                                           <  SRR4302537.1577018/94‑1 (MQ=60)
                                                   TATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAC                                           >  SRR4302537.1159161/1‑94 (MQ=60)
                                                   TATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAC                                           <  SRR4302537.936673/94‑1 (MQ=60)
                                                   TATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAC                                           <  SRR4302537.1039732/94‑1 (MQ=60)
                                                    ATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACC                                          <  SRR4302537.3431207/94‑1 (MQ=60)
                                                    ATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACC                                          >  SRR4302537.222409/1‑94 (MQ=60)
                                                    ATCGACGAAGAGTTTAAAGGTACAGGCAGCATGGATCTGCTCCTCTCTCGTAACATCGCTCAAAACCGCGTCTGCCCGGCTATCGACTACAACC                                          <  SRR4302537.3358428/94‑1 (MQ=60)
                                                      CTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGT                                        >  SRR4302537.1250334/1‑94 (MQ=60)
                                                      CTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGT                                        <  SRR4302537.1868923/94‑1 (MQ=60)
                                                        ACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTC                                      >  SRR4302537.2531507/1‑94 (MQ=60)
                                                         CGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCT                                     <  SRR4302537.1881662/94‑1 (MQ=60)
                                                         CGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCT                                     <  SRR4302537.2790169/94‑1 (MQ=60)
                                                         CGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCT                                     >  SRR4302537.4924496/1‑94 (MQ=60)
                                                          GAAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTG                                    >  SRR4302537.4800060/1‑94 (MQ=60)
                                                           AAGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGG                                   <  SRR4302537.1868905/94‑1 (MQ=60)
                                                            AGAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGT                                  >  SRR4302537.1187883/1‑94 (MQ=60)
                                                             GAGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTA                                 <  SRR4302537.3921804/94‑1 (MQ=60)
                                                              AGTTTAAAGGTACAGTCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAC                                <  SRR4302537.332971/94‑1 (MQ=60)
                                                              AGTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAC                                <  SRR4302537.3679830/94‑1 (MQ=60)
                                                               GTTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACC                               >  SRR4302537.1462856/1‑94 (MQ=60)
                                                                TTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCC                              <  SRR4302537.244418/94‑1 (MQ=60)
                                                                TTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCC                              >  SRR4302537.3171239/1‑94 (MQ=60)
                                                                TTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCC                              >  SRR4302537.3171238/1‑94 (MQ=60)
                                                                TTTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCC                              >  SRR4302537.704031/1‑94 (MQ=60)
                                                                 TTAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCG                             <  SRR4302537.4319059/94‑1 (MQ=60)
                                                                  TAAAGGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGT                            <  SRR4302537.3720951/94‑1 (MQ=60)
                                                                      GGTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAG                        >  SRR4302537.661159/1‑94 (MQ=60)
                                                                       GTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGA                       <  SRR4302537.143476/94‑1 (MQ=60)
                                                                       GTACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGA                       <  SRR4302537.5114210/94‑1 (MQ=60)
                                                                        TACAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAA                      <  SRR4302537.3570964/94‑1 (MQ=60)
                                                                          CAGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCAACTACAACCGTTCTGGTACCCGTAAAGAAGA                    >  SRR4302537.5633780/1‑94 (MQ=60)
                                                                          CAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGA                    <  SRR4302537.3757828/94‑1 (MQ=60)
                                                                           AGGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAG                   >  SRR4302537.6503917/1‑94 (MQ=60)
                                                                            GGCAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGC                  >  SRR4302537.5686034/1‑94 (MQ=60)
                                                                              CAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTG                <  SRR4302537.562889/94‑1 (MQ=60)
                                                                              CAACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTACCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTG                <  SRR4302537.1121917/94‑1 (MQ=60)
                                                                               AACATGGAACTGCTCCTCTCTNGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGC               <  SRR4302537.2585310/94‑1 (MQ=60)
                                                                               AACATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGC               >  SRR4302537.2797942/1‑94 (MQ=60)
                                                                                 CATGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTC             >  SRR4302537.3390580/1‑94 (MQ=60)
                                                                                 CATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTC             >  SRR4302537.6107928/1‑94 (MQ=60)
                                                                                   TGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCAC           <  SRR4302537.4650766/94‑1 (MQ=60)
                                                                                   TGGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCAC           >  SRR4302537.2247108/1‑94 (MQ=60)
                                                                                    GGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACG          >  SRR4302537.1817467/1‑94 (MQ=60)
                                                                                    GGAACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACG          >  SRR4302537.1676729/1‑94 (MQ=60)
                                                                                      AACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGAC        >  SRR4302537.6167294/1‑94 (MQ=60)
                                                                                       ACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCAagaat       >  SRR4302537.4364103/1‑89 (MQ=60)
                                                                                       ACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACT       >  SRR4302537.4448463/1‑94 (MQ=60)
                                                                                       ACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACT       <  SRR4302537.4396562/94‑1 (MQ=60)
                                                                                       ACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACT       >  SRR4302537.6254700/1‑94 (MQ=60)
                                                                                       ACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCAGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACT       >  SRR4302537.512444/1‑94 (MQ=60)
                                                                                       ACTGCTCCTCTCTCGTAAGATCGCTGAAAAACGCCTCTTCCCGCCTATCGACTACAACCGTCCTGGTACCCGTAAAGAAGAGCTGCTCACGACT       <  SRR4302537.5803767/94‑1 (MQ=60)
                                                                                         TGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCA     >  SRR4302537.916630/1‑94 (MQ=60)
                                                                                         TGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCA     >  SRR4302537.3842213/1‑94 (MQ=60)
                                                                                         TGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCA     >  SRR4302537.1798/1‑94 (MQ=60)
                                                                                         TGCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCA     >  SRR4302537.4631579/1‑94 (MQ=60)
                                                                                          GCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAG    >  SRR4302537.274929/1‑94 (MQ=60)
                                                                                          GCTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAG    >  SRR4302537.4083402/1‑94 (MQ=60)
                                                                                           CTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGG   >  SRR4302537.5287181/1‑94 (MQ=60)
                                                                                           CTCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGG   >  SRR4302537.2775435/1‑94 (MQ=60)
                                                                                            TCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGA  <  SRR4302537.5588302/94‑1 (MQ=60)
                                                                                            |                                                                                             
GACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGA  >  NC_000913/3964970‑3965155

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