Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,946 A→G 7.5% intergenic (+19/+17) mcbR → / ← yncD DNA‑binding transcriptional dual regulator McbR/putative TonB‑dependent outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,9460AG 7.5% 236.7 / 7.4 80intergenic (+19/+17)mcbR/yncDDNA‑binding transcriptional dual regulator McbR/putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base G (4/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACAC  >  NC_000913/1520798‑1521093
                                                                                                                                                    |                                                                                                                                                   
ggAGCAAATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATAc                                                                                                                                                    >  2:405400/1‑150 (MQ=255)
   gCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCgg                                                                                                                                                 <  1:894671/150‑1 (MQ=255)
        taGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATc                                                                                                                                            <  1:1178083/150‑1 (MQ=255)
         aGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCa                                                                                                                                           >  2:382621/1‑150 (MQ=255)
         aGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCa                                                                                                                                           >  2:1204679/1‑150 (MQ=255)
          gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCgg                                                                                                                                                 <  1:279928/143‑1 (MQ=255)
          gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCgg                                                                                                                                                 >  2:279928/1‑143 (MQ=255)
          gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCaa                                                                                                                                          <  1:1428210/150‑1 (MQ=255)
            aaaCTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAtttact                                                                                                                                        <  1:183380/150‑6 (MQ=255)
            aaaCTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAgg                                                                                                                                        <  2:1043047/150‑1 (MQ=255)
                          tACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtc                                                                                                                          <  1:644878/150‑1 (MQ=255)
                          tACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAta                                                                                                                          >  1:491231/1‑149 (MQ=255)
                           aCTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtct                                                                                                                         >  1:1033986/1‑150 (MQ=255)
                              cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAtt                                                                                                                                            <  2:1181029/128‑2 (MQ=255)
                              cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAtt                                                                                                                                            >  1:1181029/1‑127 (MQ=255)
                                         aaGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCa                                                                                                           <  1:246138/150‑1 (MQ=255)
                                            cAAAATACAAAGATGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                        <  2:188265/150‑1 (MQ=255)
                                             aaaaGACAAAGAGGGATGCAGACATTGTATTGCTGAAATTATGCAACAAAATATTGCTATTTTATACGAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAATGTAATTACTCAAATCTCCACGCAATATTCATGcc                                                                                                       <  2:336732/150‑1 (MQ=255)
                                              aaaGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCg                                                                                                      >  1:413732/1‑150 (MQ=255)
                                              aaaGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCg                                                                                                      <  2:94444/150‑1 (MQ=255)
                                                aGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGac                                                                                                    <  2:208206/150‑1 (MQ=255)
                                                   cAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCCGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCACCACGCAATATTCATGCCGACAcc                                                                                                 >  2:1251681/1‑150 (MQ=255)
                                                   cAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACAcc                                                                                                 <  1:302616/150‑1 (MQ=255)
                                                     aaGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCAt                                                                                               >  1:98080/1‑150 (MQ=255)
                                                            aTGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCg                                                                                        <  1:586928/150‑1 (MQ=255)
                                                               cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                                     >  2:834541/1‑150 (MQ=255)
                                                               cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                                     >  2:954123/1‑150 (MQ=255)
                                                                         cTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGt                                                                           >  2:58364/1‑150 (MQ=255)
                                                                             cTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGt                                                                       >  2:740667/1‑150 (MQ=255)
                                                                                aaaTTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAat                                                                    >  1:1287059/1‑150 (MQ=255)
                                                                                 aaTTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAata                                                                   >  2:637245/1‑150 (MQ=255)
                                                                                     aTGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGc                                                               >  2:353844/1‑150 (MQ=255)
                                                                                       gCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGccc                                                             >  1:904395/1‑150 (MQ=255)
                                                                                        caacaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCg                                                            >  1:1035023/1‑150 (MQ=255)
                                                                                          acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                          <  2:586131/150‑1 (MQ=255)
                                                                                            aaaaTATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTg                                                        <  1:1204679/150‑1 (MQ=255)
                                                                                             aaaTATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGa                                                       >  2:1168547/1‑150 (MQ=255)
                                                                                              aaTATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGAc                                                      <  1:1165328/150‑1 (MQ=255)
                                                                                              aaTATTGCTATTTTATACCAGCAATAAAATCGTTAAAGTAAATCCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCAATCAATATTCATGCCGAGACCATTATTTCGTCCTGGCGAAGGTTCGTAATATCGCCCGTTTGAc                                                      <  2:296636/150‑1 (MQ=255)
                                                                                                   tGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCAtt                                                 <  1:523291/150‑1 (MQ=255)
                                                                                                    gCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTg                                                >  1:730303/1‑150 (MQ=255)
                                                                                                    gCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTg                                                <  2:1033986/150‑1 (MQ=255)
                                                                                                     cTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGa                                               <  1:1142188/150‑1 (MQ=255)
                                                                                                             taCCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACAGCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGAc                                        <  1:499714/149‑1 (MQ=255)
                                                                                                             taCCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACAGCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGAc                                        >  2:499714/1‑149 (MQ=255)
                                                                                                              aCCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAg                                      <  1:969484/150‑1 (MQ=255)
                                                                                                               ccAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGa                                     >  1:671413/1‑150 (MQ=255)
                                                                                                                cAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGaa                                    >  2:720985/1‑150 (MQ=255)
                                                                                                                 aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCTAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAc                                   <  2:141344/150‑1 (MQ=255)
                                                                                                                 aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAc                                   <  2:267140/150‑1 (MQ=255)
                                                                                                                    aaTACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGc                                                               <  2:1293314/119‑1 (MQ=255)
                                                                                                                    aaTACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGc                                                               >  1:1293314/1‑119 (MQ=255)
                                                                                                                     aTACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAAc                               >  2:1375050/1‑150 (MQ=255)
                                                                                                                        cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCAtt                                                 <  2:1357426/129‑1 (MQ=255)
                                                                                                                        cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCAtt                                                 >  1:1357426/1‑129 (MQ=255)
                                                                                                                         aaTCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTAt                           >  1:706977/1‑150 (MQ=255)
                                                                                                                           tCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTc                         <  2:1035023/150‑1 (MQ=255)
                                                                                                                            cGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATCTTCATGCCGCCCCCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTTACAATGACAGAACCAACGTATTCt                        >  2:652519/1‑150 (MQ=255)
                                                                                                                            cGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCt                        <  2:1067905/150‑1 (MQ=255)
                                                                                                                             gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACa                                       >  1:385928/1‑134 (MQ=255)
                                                                                                                             gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACa                                       <  2:385928/134‑1 (MQ=255)
                                                                                                                             gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGCACCAACGTATTCtt                       >  1:1272588/1‑150 (MQ=255)
                                                                                                                             gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCtt                       >  1:652348/1‑150 (MQ=255)
                                                                                                                             gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCtt                       >  1:320513/1‑150 (MQ=255)
                                                                                                                             gTTAAAGTAAATGCCCGGCAAATACCGGGAATCAAGGTAATTACTCAAATCTCCAGGCATTATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATTTCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCtt                       <  2:366016/150‑1 (MQ=255)
                                                                                                                              ttAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCttt                      <  2:1067329/150‑1 (MQ=255)
                                                                                                                               tAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTa                     <  2:491231/150‑1 (MQ=255)
                                                                                                                                aaaGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTAt                    <  1:1168547/150‑1 (MQ=255)
                                                                                                                                  aGTAAATGCCCGGCAAATACCGGGCATCAATGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCg                  <  1:715388/150‑1 (MQ=255)
                                                                                                                                  aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCg                  >  2:175744/1‑150 (MQ=255)
                                                                                                                                  aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCg                  <  1:721433/150‑1 (MQ=255)
                                                                                                                                   gTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGa                 <  2:904395/150‑1 (MQ=255)
                                                                                                                                   cttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGa                 >  1:888326/5‑150 (MQ=255)
                                                                                                                                      aaTGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATa              >  2:914487/1‑150 (MQ=255)
                                                                                                                                       aTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATaa             >  1:814552/1‑150 (MQ=255)
                                                                                                                                        tGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTAACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATaaa            >  2:250699/1‑150 (MQ=255)
                                                                                                                                        tGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATaaa            >  1:955184/1‑150 (MQ=255)
                                                                                                                                         gCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAAt           <  1:353844/150‑1 (MQ=255)
                                                                                                                                               cAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGa     >  1:158565/1‑150 (MQ=255)
                                                                                                                                                  aTACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGacac  <  1:1082473/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACAC  >  NC_000913/1520798‑1521093

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: