Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,772,236 A→T 6.4% pseudogene (947/1016 nt) yfjV ← CP4‑57 prophage; putative arsenite transporter, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,772,2360AT 6.4% 324.0 / 8.1 103pseudogene (947/1016 nt)yfjVCP4‑57 prophage; putative arsenite transporter, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base A (48/48);  new base T (2/5);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 104 total aligned reads displayed.

GCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTT  >  NC_000913/2772091‑2772385
                                                                                                                                                 |                                                                                                                                                     
gCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCTGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGc                                                                                                                                                   <  1:523148/150‑1 (MQ=255)
 cATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGcc                                                                                                                                                  <  1:1499344/150‑1 (MQ=255)
 cATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGcc                                                                                                                                                  <  1:910552/150‑1 (MQ=255)
    ctcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGt                                                                                                                                               <  1:1064054/150‑1 (MQ=255)
    ctcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGt                                                                                                                                               >  2:224510/1‑150 (MQ=255)
     tcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtg                                                                                                                                              <  1:1082569/150‑2 (MQ=255)
     tcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTa                                                                                                                                              <  2:470433/150‑1 (MQ=255)
      ccAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTAc                                                                                                                                             >  1:323511/1‑150 (MQ=255)
       cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgct                                                                                                                                            <  2:1065557/150‑4 (MQ=21)
         gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctgt                                                                                                                                          <  1:170275/150‑6 (MQ=21)
          cGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACATAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTTTCATGACTGTGCCGGtgctgtt                                                                                                                                         <  1:100529/150‑7 (MQ=21)
               tGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAat                                                                                                                                    >  1:546816/1‑150 (MQ=255)
               tGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGCCACAGTCATGACAATGCCGGTACGGAAGTAat                                                                                                                                    >  2:675145/1‑150 (MQ=255)
                gCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAata                                                                                                                                   >  2:304843/1‑150 (MQ=255)
                gCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAana                                                                                                                                   <  1:218952/150‑1 (MQ=255)
                 cGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAatat                                                                                                                                  <  1:1494997/150‑1 (MQ=255)
                     aGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcccc                                                                                                                              <  1:304843/150‑1 (MQ=255)
                         gCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGt                                                                                                                          >  2:884582/1‑150 (MQ=255)
                          ccACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTg                                                                                                                         <  1:1091934/150‑1 (MQ=255)
                           cacaGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGa                                                                                                                        <  1:998392/150‑1 (MQ=255)
                               ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGt                                                                                                                    >  1:1322034/1‑150 (MQ=255)
                                gCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTc                                                                                                                   >  2:334266/1‑150 (MQ=255)
                                  ggggTTGTGTTAGATGGTAATGTTGCTCATATCTGTATGTAATTAAATAGTGAAAGAGAGTGGCCACGCAAGGGCGGGCTGAGTGACAAACAGCTACGTCACAGTTATGACTATGACGGTACGGATGTAATTTACCGTAGTGTTTGTCat                                                                                                                 <  2:1504864/149‑1 (MQ=255)
                                  cGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCat                                                                                                                 >  2:138065/1‑150 (MQ=255)
                                       tgtgATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGCGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtttt                                                                                                            >  2:1313981/1‑150 (MQ=255)
                                        gtgATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttt                                                                                                           >  1:846566/1‑150 (MQ=255)
                                          gATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCt                                                                                                         >  2:1533270/1‑150 (MQ=255)
                                          gATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCt                                                                                                         >  1:309192/1‑150 (MQ=255)
                                          gATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTATCATATTTTTCt                                                                                                         >  2:68021/1‑150 (MQ=255)
                                           aTAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                                                                        <  2:833831/150‑1 (MQ=255)
                                            tAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTgg                                                                                                       <  2:1289087/150‑1 (MQ=255)
                                            tAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTgg                                                                                                       >  2:150788/1‑150 (MQ=255)
                                             aGATGGTAATGTTGCTCATATCAGTATCTACTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCCCCGGCACACGCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                                                                      <  2:1307391/150‑1 (MQ=255)
                                             aGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                                                                      >  1:577495/1‑150 (MQ=255)
                                                  gTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTTACAAACAGCACCGTCACAGTCATTACAATGCCGGTACTGAAGTAATATACACAAGTGATTGTCATATTTTTCTGGTCAAgt                                                                                                 <  2:658707/150‑2 (MQ=25)
                                                  gTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGc                                                                                                 <  1:657668/150‑1 (MQ=255)
                                                       gTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                                                                            >  2:305624/1‑150 (MQ=255)
                                                        ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGc                                                                                           >  2:251381/1‑150 (MQ=255)
                                                         tgatCATATCAGTATCTAATTACAAAGTGACATAGAGCCGCGACGACAGGGCGGCCATAGTGACAAACAGCACGGGCACAGTCATGACAATGGCGGTACGGAAGTAATATCCCCATGTGATTGTTATATTTTTCTTGGGAAGCACATggt                                                                                          <  2:323511/147‑3 (MQ=21)
                                                          gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAg                                                                                         <  2:1225284/150‑1 (MQ=255)
                                                               tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcaca                                                                                    <  1:1533270/150‑1 (MQ=255)
                                                               tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcaca                                                                                    <  1:1201786/150‑1 (MQ=255)
                                                                gtCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcag                                                                                   >  2:1040048/2‑150 (MQ=255)
                                                                 tCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagc                                                                                  <  2:11374/150‑1 (MQ=255)
                                                                 tCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagc                                                                                  >  1:1029845/1‑150 (MQ=255)
                                                                   aGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagcag                                                                                >  2:435128/1‑150 (MQ=255)
                                                                   aGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagcag                                                                                >  1:1381109/1‑150 (MQ=255)
                                                                    gTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCaga                                                                               >  1:1376572/1‑150 (MQ=255)
                                                                     tATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCagag                                                                              <  1:31545/150‑1 (MQ=255)
                                                                       tCTAATTAAAAAGTGACAGAGAGACTCCACGACAGCGCGGACAGAGTGACAAACAGCACCGGCACAGTCATGACAATACCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGtt                                                                            >  2:1173155/1‑150 (MQ=255)
                                                                        cTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTg                                                                           <  1:344049/150‑1 (MQ=255)
                                                                        cTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTg                                                                           >  2:631339/1‑150 (MQ=255)
                                                                         tAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGc                                                                          >  2:1228753/1‑150 (MQ=255)
                                                                          aaTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGcc                                                                         <  1:571089/150‑1 (MQ=255)
                                                                           aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCGTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                                                                        >  1:1348049/1‑150 (MQ=255)
                                                                            ttACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAg                                                                       <  1:483096/150‑1 (MQ=255)
                                                                              aCAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGAc                                                                     >  2:727377/1‑150 (MQ=255)
                                                                                    tGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCaa                                                               >  1:669/1‑150 (MQ=255)
                                                                                      aCAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATc                                                             >  2:129716/1‑150 (MQ=255)
                                                                                       cAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCg                                                            >  1:979483/1‑150 (MQ=255)
                                                                                             gccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGa                                                      >  1:223939/1‑150 (MQ=255)
                                                                                                     gccaGCGCGGCCAGAGTGACAAACAGCACCGGCATAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCggg                                              <  1:659388/150‑1 (MQ=255)
                                                                                                      ccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGc                                             <  2:871524/150‑1 (MQ=255)
                                                                                                      ccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGc                                             <  2:846566/150‑1 (MQ=255)
                                                                                                      ccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGc                                             >  1:650458/1‑150 (MQ=255)
                                                                                                      ccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGc                                             <  1:138065/150‑1 (MQ=255)
                                                                                                       caGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGcc                                            <  2:550831/150‑1 (MQ=255)
                                                                                                          cgcgGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTaa                                         <  1:178948/150‑1 (MQ=255)
                                                                                                              gCCAGAGTGACAAACAGCACCGGCACAGTCGTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATca                                     <  2:1348049/150‑1 (MQ=255)
                                                                                                              gCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATca                                     >  1:33679/1‑150 (MQ=255)
                                                                                                               ccAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGt                                                        >  1:1459970/1‑130 (MQ=255)
                                                                                                               ccAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGt                                                        <  2:1459970/130‑1 (MQ=255)
                                                                                                               ccAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATcac                                    <  1:1246158/150‑1 (MQ=255)
                                                                                                                cAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATcaca                                   <  2:1322034/150‑1 (MQ=255)
                                                                                                                 agagTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAg                                  <  2:1029845/150‑1 (MQ=255)
                                                                                                                        cAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                           >  1:1158118/1‑150 (MQ=255)
                                                                                                                        cAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                           >  2:1052011/1‑150 (MQ=255)
                                                                                                                            cAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc                                                                 <  1:308772/108‑1 (MQ=255)
                                                                                                                            cAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc                                                                 >  2:308772/1‑108 (MQ=255)
                                                                                                                             cgtaCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTg                      >  2:1446149/4‑150 (MQ=37)
                                                                                                                              gCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTgg                     >  1:1200657/1‑150 (MQ=255)
                                                                                                                              gCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTgg                     >  1:870960/1‑150 (MQ=255)
                                                                                                                                aCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCa                   >  1:49957/1‑150 (MQ=255)
                                                                                                                                aCCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCa                   >  2:365110/1‑150 (MQ=255)
                                                                                                                                 ccGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCAt                  <  2:358236/150‑1 (MQ=255)
                                                                                                                                  cGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATGTGATTGGctgt                 <  2:937051/150‑4 (MQ=255)
                                                                                                                                  cGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATa                 <  1:588192/150‑1 (MQ=255)
                                                                                                                                   ggCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTATCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATGTGATTGGctgtg                <  1:68021/150‑5 (MQ=255)
                                                                                                                                     cacaGTCATGACAATGCCGGTACGGAAGTAATATCCCCCAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTGGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATGTGATTGGctgtgat              <  2:661147/150‑7 (MQ=25)
                                                                                                                                     cacaGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATGTGATTGGctgtgat              <  1:221772/150‑7 (MQ=37)
                                                                                                                                     cacaGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATGTGATTGGctgtgat              <  2:1416167/150‑7 (MQ=37)
                                                                                                                                     cacaGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAAt              >  1:542900/1‑150 (MQ=255)
                                                                                                                                      acaGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATGTGATTGGctgtgatt             <  1:334266/150‑8 (MQ=37)
                                                                                                                                      acaGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATc             <  2:570093/150‑1 (MQ=255)
                                                                                                                                           cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGc        >  1:1428764/1‑150 (MQ=255)
                                                                                                                                           cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGc        >  2:272610/1‑150 (MQ=255)
                                                                                                                                             tGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCct      >  2:1280934/1‑150 (MQ=255)
                                                                                                                                              gACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCctc     <  1:272610/150‑1 (MQ=255)
                                                                                                                                                cAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCtt   >  1:1528762/1‑150 (MQ=255)
                                                                                                                                                 aaTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCttt  >  2:1224271/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTT  >  NC_000913/2772091‑2772385

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