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 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,959,281 G→A E330E (GAG→GAA rep → rep helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,2810GA99.5% 799.3 / ‑5.0 218E330E (GAG→GAAreprep helicase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (112/104);  total (113/105)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

Only 100 of 221 total aligned reads displayed.

GCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAA  >  NC_000913/3959191‑3959370
                                                                                          |                                                                                         
  gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGa                                                                                       >  1:477603/1‑93 (MQ=255)
   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      >  1:2544758/1‑93 (MQ=255)
   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      <  2:3540831/93‑1 (MQ=255)
   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      >  1:3514747/1‑93 (MQ=255)
   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      >  1:2528838/1‑93 (MQ=255)
    cTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATc                                                                                     >  2:4459604/1‑93 (MQ=255)
      gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGc                                                                                   >  1:947240/1‑93 (MQ=255)
      gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGc                                                                                   >  1:4885507/1‑93 (MQ=255)
       ttCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGcc                                                                                  >  1:2928046/1‑93 (MQ=255)
         ctcCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCa                                                                                >  2:1089190/1‑93 (MQ=255)
           ccGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATc                                                                              >  2:3431735/1‑93 (MQ=255)
             gAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCAc                                                                            >  1:3078695/1‑93 (MQ=255)
               aCTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACtt                                                                          >  1:5105157/1‑93 (MQ=255)
                cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTc                                                                         >  2:2513834/1‑93 (MQ=255)
                cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTc                                                                         >  2:566319/1‑93 (MQ=255)
                  gggTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGt                                                                       >  1:3043888/1‑93 (MQ=255)
                    gTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCa                                                                     >  1:5999390/1‑93 (MQ=255)
                     ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCaa                                                                    <  2:5867940/93‑1 (MQ=255)
                     ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCaa                                                                    <  1:463899/93‑1 (MQ=255)
                     ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCaa                                                                    <  1:4877585/93‑1 (MQ=255)
                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  <  1:5563972/93‑1 (MQ=255)
                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  >  1:4467734/1‑93 (MQ=255)
                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  <  2:5900178/93‑1 (MQ=255)
                          gcgcGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATaaaa                                                               >  1:1037537/1‑93 (MQ=255)
                          gcgcGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATaaaa                                                               <  2:624090/93‑1 (MQ=255)
                            gcgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACg                                                             >  2:319306/1‑93 (MQ=255)
                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            <  2:1764367/93‑1 (MQ=255)
                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            <  2:3624879/93‑1 (MQ=255)
                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            >  1:3102720/1‑93 (MQ=255)
                              ggAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCa                                                           <  1:3837008/93‑1 (MQ=255)
                              ggAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCa                                                           >  1:1964978/1‑93 (MQ=255)
                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  2:2131001/93‑1 (MQ=255)
                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  2:3978443/93‑1 (MQ=255)
                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          >  1:4751568/1‑93 (MQ=255)
                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  2:5043216/93‑1 (MQ=255)
                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  1:482350/93‑1 (MQ=255)
                                  cTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAc                                                       >  1:773437/1‑93 (MQ=255)
                                  cTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAc                                                       <  2:3444080/93‑1 (MQ=255)
                                  cTAAAAGAATTAAGCGCGAATAACGAAGAACATGAGGCTGAGGGCGTTAATGGCGAACTGATGCCCCATCACTTCGTCAATAAAACGCAgcaa                                                       >  2:800717/1‑90 (MQ=255)
                                     aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaaa                                                    >  1:4390450/1‑93 (MQ=255)
                                     aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaaa                                                    >  2:4042052/1‑93 (MQ=255)
                                     aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaaa                                                    <  1:5577308/93‑1 (MQ=255)
                                      aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAg                                                   >  1:6102937/1‑93 (MQ=255)
                                      aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATAACTTCGTCAATAAAACGCAGTACaaaa                                                   >  2:1529083/1‑92 (MQ=255)
                                       atcaTTGAGCGCGGATAGCGCAGTGCATGCGTCTGTGCGCGTTACTGGCAAACTGATCGCCCATCACTTCCTCAATAAAACGCAGTACAAAgc                                                  <  2:597851/90‑2 (MQ=255)
                                       aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGa                                                  <  2:305991/93‑1 (MQ=255)
                                       aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGa                                                  <  1:4509227/93‑1 (MQ=255)
                                       aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGa                                                  >  2:5285495/1‑93 (MQ=255)
                                        gTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAt                                                 >  1:1059270/1‑93 (MQ=255)
                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACCAAGAtt                                                >  2:1346169/1‑93 (MQ=255)
                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                >  1:677663/1‑93 (MQ=255)
                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                >  2:208659/1‑93 (MQ=255)
                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                <  1:5064109/93‑1 (MQ=255)
                                          aTTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTa                                               <  2:931703/93‑1 (MQ=255)
                                          aTTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTa                                               >  1:5234186/1‑93 (MQ=255)
                                           ttAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAc                                              <  2:614465/93‑1 (MQ=255)
                                           ttAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAc                                              <  1:1930576/93‑1 (MQ=255)
                                              aGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGcc                                           >  1:5790405/1‑93 (MQ=255)
                                                  cgAATAACGACGAACATGAGGCTGACCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                       <  1:2570518/93‑1 (MQ=255)
                                                  cgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCAttt                                       >  2:4258213/1‑92 (MQ=255)
                                                  cgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                       <  1:4613198/93‑1 (MQ=255)
                                                    aaTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCtt                                     <  1:4132438/93‑1 (MQ=255)
                                                     aTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCttt                                    <  1:5914247/93‑1 (MQ=255)
                                                      tAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTa                                   <  1:5397594/93‑1 (MQ=255)
                                                       aaCGAAGAACATGAGGCTGAGCGCGTTACTGGCCAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAt                                  <  1:5935652/93‑1 (MQ=255)
                                                          gaagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgc                               >  2:2656834/1‑93 (MQ=255)
                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              >  2:282422/1‑93 (MQ=255)
                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              <  2:2177350/93‑1 (MQ=255)
                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              <  1:2666532/93‑1 (MQ=255)
                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              <  2:2459970/93‑1 (MQ=255)
                                                            agaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCgg                             >  2:4817365/1‑93 (MQ=255)
                                                            agaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAACAAGATTACGCCATTCTTTATCGCgg                             >  2:3733857/1‑93 (MQ=255)
                                                             gaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGt                            <  2:1451426/93‑1 (MQ=255)
                                                               aCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTaa                          >  2:2102552/1‑93 (MQ=255)
                                                                cATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAc                         >  2:2382915/1‑93 (MQ=255)
                                                                cATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAc                         >  1:4093706/1‑93 (MQ=255)
                                                                 aTGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAcc                        <  1:442862/93‑1 (MQ=255)
                                                                 aTGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAcc                        <  2:3598109/93‑1 (MQ=255)
                                                                  tGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCa                       >  2:4107612/1‑93 (MQ=255)
                                                                   gAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCAt                      <  2:4214893/93‑1 (MQ=255)
                                                                      gCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAg                   <  1:1818019/93‑1 (MQ=255)
                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  <  1:4704098/93‑1 (MQ=255)
                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  >  2:5754171/1‑93 (MQ=255)
                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  <  1:3456434/93‑1 (MQ=255)
                                                                        tGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTc                 >  2:3124887/1‑93 (MQ=255)
                                                                         gAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcg                >  1:5055877/1‑93 (MQ=255)
                                                                          aGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgc               >  1:1231497/1‑93 (MQ=255)
                                                                           gcgcGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTCACCAACAGTcgcg              >  1:1707953/1‑93 (MQ=255)
                                                                           gcgcGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgcg              >  2:5534956/1‑93 (MQ=255)
                                                                           gcgcGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgcg              >  1:1518182/1‑93 (MQ=255)
                                                                            cgcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCgg             <  1:1337902/93‑1 (MQ=255)
                                                                             gcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCggg            <  2:3774748/93‑1 (MQ=255)
                                                                               gTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtg          >  2:1006375/1‑93 (MQ=255)
                                                                                ttACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtgt         >  2:595764/1‑93 (MQ=255)
                                                                                ttACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtgt         >  2:2245058/1‑93 (MQ=255)
                                                                                  aCTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGttt       <  2:4723633/93‑1 (MQ=255)
                                                                                  aCTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGttt       >  2:6049859/1‑93 (MQ=255)
                                                                                    tGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGa     >  2:4646822/1‑93 (MQ=255)
                                                                                     ggCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGaa    >  1:2288396/1‑93 (MQ=255)
                                                                                       cGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGaaaa  <  1:2086246/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAA  >  NC_000913/3959191‑3959370

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ 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
AAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTC  >  NC_000913/3959188‑3959374
                                                                                             |                                                                                             
AAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAA                                                                                               >  SRR4302532.2595731/1‑94 (MQ=60)
AAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAA                                                                                               >  SRR4302532.6170600/1‑94 (MQ=60)
 AAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAAC                                                                                              <  SRR4302532.3368644/94‑1 (MQ=60)
  AGCGTCTGTCCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACT                                                                                             <  SRR4302532.5437106/94‑1 (MQ=60)
   GCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTG                                                                                            <  SRR4302532.4770647/94‑1 (MQ=60)
    CGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGA                                                                                           >  SRR4302532.494257/1‑94 (MQ=60)
     GTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAT                                                                                          >  SRR4302532.3632426/1‑94 (MQ=60)
     GTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAT                                                                                          >  SRR4302532.388267/1‑94 (MQ=60)
     GTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAT                                                                                          >  SRR4302532.2629632/1‑94 (MQ=60)
     GTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAT                                                                                          >  SRR4302532.2613222/1‑94 (MQ=60)
       CTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCG                                                                                        <  SRR4302532.3625238/94‑1 (MQ=60)
        TGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGC                                                                                       >  SRR4302532.979647/1‑94 (MQ=60)
        TGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGC                                                                                       >  SRR4302532.4412109/1‑94 (MQ=60)
        TGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGC                                                                                       >  SRR4302532.5053446/1‑94 (MQ=60)
         GTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCC                                                                                      >  SRR4302532.3025447/1‑94 (MQ=60)
          TTCTCCGGACTGGGATATGGCGCGGAGCTAAAAGTATTAAGAGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTAGCGAACTGATCACCC                                                                                     >  SRR4302532.4468214/1‑94 (MQ=60)
               CGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCAC                                                                                >  SRR4302532.3180440/1‑94 (MQ=60)
               CGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCAC                                                                                >  SRR4302532.4246671/1‑94 (MQ=60)
               CGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCAC                                                                                >  SRR4302532.1375488/1‑94 (MQ=60)
                GAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACT                                                                               <  SRR4302532.62958/94‑1 (MQ=60)
                 AACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTT                                                                              >  SRR4302532.5280663/1‑94 (MQ=60)
                  ACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTC                                                                             >  SRR4302532.2145917/1‑94 (MQ=60)
                    TGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGT                                                                           >  SRR4302532.3144338/1‑94 (MQ=60)
                      GGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCA                                                                         >  SRR4302532.6206595/1‑94 (MQ=60)
                        TTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAAT                                                                       <  SRR4302532.4563281/94‑1 (MQ=60)
                        TTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAAT                                                                       <  SRR4302532.480054/94‑1 (MQ=60)
                        TTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAAT                                                                       <  SRR4302532.5045282/94‑1 (MQ=60)
                         TATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCTCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATA                                                                      <  SRR4302532.299061/94‑1 (MQ=60)
                         TATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATA                                                                      >  SRR4302532.4620054/1‑94 (MQ=60)
                          ATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAA                                                                     <  SRR4302532.5292137/94‑1 (MQ=60)
                          ATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAA                                                                     <  SRR4302532.5757166/94‑1 (MQ=60)
                            GGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAA                                                                   >  SRR4302532.1072538/1‑94 (MQ=60)
                             GCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAAC                                                                  <  SRR4302532.1776559/94‑1 (MQ=60)
                               GCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGC                                                                >  SRR4302532.3205619/1‑94 (MQ=60)
                                CGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCA                                                               >  SRR4302532.2030244/1‑94 (MQ=60)
                                CGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCA                                                               >  SRR4302532.5058759/1‑94 (MQ=60)
                                 GGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAG                                                              >  SRR4302532.1630754/1‑94 (MQ=60)
                                 GGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAG                                                              >  SRR4302532.4914028/1‑94 (MQ=60)
                                 GGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAG                                                              <  SRR4302532.3966553/94‑1 (MQ=60)
                                 GGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAG                                                              <  SRR4302532.4792693/94‑1 (MQ=60)
                                  GAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGT                                                             <  SRR4302532.499180/94‑1 (MQ=60)
                                  GAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGT                                                             >  SRR4302532.1886272/1‑94 (MQ=60)
                                   AGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGACCATCACTTCGTCAATAAAACGCAGTA                                                            >  SRR4302532.1492037/1‑94 (MQ=60)
                                    GCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAC                                                           >  SRR4302532.2515579/1‑94 (MQ=60)
                                    GCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAC                                                           >  SRR4302532.800101/1‑94 (MQ=60)
                                    GCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAC                                                           >  SRR4302532.800103/1‑94 (MQ=60)
                                     CTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACA                                                          >  SRR4302532.5965553/1‑94 (MQ=60)
                                       AAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAA                                                        >  SRR4302532.4539718/1‑94 (MQ=60)
                                        AAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAG                                                       <  SRR4302532.5770944/94‑1 (MQ=60)
                                        AAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAG                                                       >  SRR4302532.6313544/1‑94 (MQ=60)
                                          AGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAT                                                     <  SRR4302532.4663187/94‑1 (MQ=60)
                                          AGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAT                                                     <  SRR4302532.6075667/94‑1 (MQ=60)
                                          AGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAT                                                     >  SRR4302532.1094961/1‑94 (MQ=60)
                                           GTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATT                                                    >  SRR4302532.700970/1‑94 (MQ=60)
                                           GTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATT                                                    >  SRR4302532.5244696/1‑94 (MQ=60)
                                           GTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATT                                                    <  SRR4302532.5047273/94‑1 (MQ=60)
                                            TATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTA                                                   <  SRR4302532.5987934/94‑1 (MQ=60)
                                            TATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTA                                                   <  SRR4302532.5238005/94‑1 (MQ=60)
                                            TATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTA                                                   >  SRR4302532.5414682/1‑94 (MQ=60)
                                             ATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCGATAAAACGCAGTACAAAGATTAC                                                  >  SRR4302532.4640626/1‑94 (MQ=60)
                                             ATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAC                                                  >  SRR4302532.5292609/1‑94 (MQ=60)
                                             ATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAC                                                  <  SRR4302532.523958/94‑1 (MQ=60)
                                             ATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAC                                                  <  SRR4302532.3583579/94‑1 (MQ=60)
                                              TTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACG                                                 <  SRR4302532.2771238/94‑1 (MQ=60)
                                              TTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACG                                                 <  SRR4302532.1994833/94‑1 (MQ=60)
                                                AAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCC                                               >  SRR4302532.5991730/1‑94 (MQ=60)
                                                   CGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAACGATTACCCCATT                                            >  SRR4302532.3348816/1‑94 (MQ=60)
                                                   CGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATT                                            >  SRR4302532.1744035/1‑94 (MQ=60)
                                                   CGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATT                                            <  SRR4302532.6340383/94‑1 (MQ=60)
                                                     CGAATAACGACGAACATGAGGCTGACCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCT                                          <  SRR4302532.2656341/94‑1 (MQ=60)
                                                     CGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCT                                          <  SRR4302532.4770500/94‑1 (MQ=60)
                                                     CAAATAACGAAGAGAGTGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCT                                          <  SRR4302532.2656340/94‑1 (MQ=60)
                                                      GAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTT                                         >  SRR4302532.4997965/1‑94 (MQ=60)
                                                       AATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTT                                        <  SRR4302532.4271072/94‑1 (MQ=60)
                                                        ATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTA                                       <  SRR4302532.6119046/94‑1 (MQ=60)
                                                         TAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAT                                      <  SRR4302532.2322533/94‑1 (MQ=60)
                                                         TAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAT                                      <  SRR4302532.5584615/94‑1 (MQ=60)
                                                         TAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAT                                      <  SRR4302532.4823465/94‑1 (MQ=60)
                                                         TAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAT                                      <  SRR4302532.6134134/94‑1 (MQ=60)
                                                         TAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAT                                      <  SRR4302532.2260169/94‑1 (MQ=60)
                                                          AACGAAGAACATGAGGCTGAGCGCGTTACTGGCCAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATC                                     <  SRR4302532.6141069/94‑1 (MQ=60)
                                                             GAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCG                                  <  SRR4302532.557866/94‑1 (MQ=60)
                                                              AAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGG                                 <  SRR4302532.2756034/94‑1 (MQ=60)
                                                              AAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGG                                 <  SRR4302532.2667764/94‑1 (MQ=60)
                                                              AAGAACATGAGGCCGAGCGCGTTACTGGCGAACTGATCGCCCAGCACTTCGTCAACAAAACGCAGTACAAAGATTACGCCATCCACTATCGCGG                                 >  SRR4302532.5679448/1‑94 (MQ=60)
                                                               AGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGT                                >  SRR4302532.4724527/1‑94 (MQ=60)
                                                                  ACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAC                             >  SRR4302532.4231135/1‑94 (MQ=60)
                                                                  ACATGAGGCTGAGCGCGTTACTGGCGAACTGATAGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTCAC                             >  SRR4302532.49054/1‑94 (MQ=60)
                                                                   CATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACC                            >  SRR4302532.932670/1‑94 (MQ=60)
                                                                    ATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCA                           <  SRR4302532.458179/94‑1 (MQ=60)
                                                                    ATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCA                           <  SRR4302532.1237227/94‑1 (MQ=60)
                                                                       AGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCA                        <  SRR4302532.2014315/94‑1 (MQ=60)
                                                                        GGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAG                       <  SRR4302532.1090615/94‑1 (MQ=60)
                                                                        GGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCCTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAG                       <  SRR4302532.4507610/94‑1 (MQ=60)
                                                                         GCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGT                      <  SRR4302532.1878622/94‑1 (MQ=60)
                                                                          CTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTGAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTC                     <  SRR4302532.3181193/94‑1 (MQ=60)
                                                                          CTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTC                     <  SRR4302532.5024185/94‑1 (MQ=60)
                                                                          CTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTC                     <  SRR4302532.3572184/94‑1 (MQ=60)
                                                                          CTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTC                     <  SRR4302532.4864587/94‑1 (MQ=60)
                                                                          CTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTC                     >  SRR4302532.3030481/1‑94 (MQ=60)
                                                                          CTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAACATCAGTC                     >  SRR4302532.6164205/1‑94 (MQ=60)
                                                                           TGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCG                    >  SRR4302532.5229403/1‑94 (MQ=60)
                                                                            GAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGC                   >  SRR4302532.1272700/1‑94 (MQ=60)
                                                                            GAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGC                   >  SRR4302532.3157629/1‑94 (MQ=60)
                                                                             AGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTCACCAACAGTCGCG                  >  SRR4302532.1765532/1‑94 (MQ=60)
                                                                             AGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCG                  >  SRR4302532.1569458/1‑94 (MQ=60)
                                                                               CGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGG                <  SRR4302532.1382642/94‑1 (MQ=60)
                                                                               CGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGG                >  SRR4302532.2623427/1‑94 (MQ=60)
                                                                                GCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGT               <  SRR4302532.2650669/94‑1 (MQ=60)
                                                                                 CGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTG              <  SRR4302532.992304/94‑1 (MQ=60)
                                                                                 CGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTG              >  SRR4302532.3393430/1‑94 (MQ=60)
                                                                                 CGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTG              <  SRR4302532.695436/94‑1 (MQ=60)
                                                                                  GTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGGGT             >  SRR4302532.417464/1‑94 (MQ=60)
                                                                                   TTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTT            >  SRR4302532.4984905/1‑94 (MQ=60)
                                                                                     ACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAGCGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTG          <  SRR4302532.5765249/94‑1 (MQ=60)
                                                                                      CTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGA         <  SRR4302532.4592557/94‑1 (MQ=60)
                                                                                       TGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAA        >  SRR4302532.2364344/1‑94 (MQ=60)
                                                                                        GGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAA       <  SRR4302532.4930678/94‑1 (MQ=60)
                                                                                          CGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGTAAAA     >  SRR4302532.4397187/1‑94 (MQ=60)
                                                                                          CGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAA     <  SRR4302532.2155809/94‑1 (MQ=60)
                                                                                           GAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAAT    <  SRR4302532.825972/94‑1 (MQ=60)
                                                                                           GAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAAT    <  SRR4302532.1031691/94‑1 (MQ=60)
                                                                                             ACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTC  >  SRR4302532.5167408/1‑94 (MQ=60)
                                                                                             ACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTC  >  SRR4302532.3049246/1‑94 (MQ=60)
                                                                                             |                                                                                             
AAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTC  >  NC_000913/3959188‑3959374

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