Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,132 T→G 8.0% A75A (GCA→GCC) ‡ cutA ← copper binding protein CutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,1320TG 8.0% 289.7 / 14.7 100A75A (GCA→GCC) ‡cutAcopper binding protein CutA
Reads supporting (aligned to +/- strand):  ref base T (55/37);  new base G (5/3);  total (60/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTA  >  NC_000913/4364984‑4365277
                                                                                                                                                    |                                                                                                                                                 
ccGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTg                                                                                                                                                  <  1:27749/150‑1 (MQ=255)
   gCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTGCAGGCATTCCAGCAGTGCct                                                                                                                                               >  1:684849/1‑150 (MQ=255)
     aaaaaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctgc                                                                                                                                             <  2:140428/150‑1 (MQ=255)
       aaaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctgctg                                                                                                                                           <  1:217096/150‑1 (MQ=255)
          cGGCAGGGCTGCAAAGTCGCAGGCTCCGCGCAACGATGCGGTGAGCCATGAGAGCTAAGCTGGGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGCTGAGACTTCAGGCATTCCAGCAGGCACTGCTggaa                                                                                                                                        <  2:1393183/150‑3 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                        >  1:878104/1‑150 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                        <  2:777421/150‑1 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGAGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                        >  2:1435294/1‑150 (MQ=255)
          cGGAAGTGCTGCAAAGTACCATGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                        <  2:1309488/150‑1 (MQ=255)
           ggAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTga                                                                                                                                       <  2:689742/150‑1 (MQ=255)
            gAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGGCACTGCTggaatg                                                                                                                                      <  1:1097019/150‑5 (MQ=255)
            gAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGGCACTGCTggaatg                                                                                                                                      <  2:842344/150‑5 (MQ=255)
             aaGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTgaga                                                                                                                                     <  2:890674/150‑1 (MQ=255)
               gTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATa                                                                                                                                   <  2:427011/150‑1 (MQ=255)
                     cAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAg                                                                                                                             >  1:381597/1‑150 (MQ=255)
                     cAAAGTAGCAGGATCAGCGTAAAGACGCGTTGAGCCATGAGAGGTAATCTGTGTCTCAGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAg                                                                                                                             >  2:726972/1‑150 (MQ=255)
                           aGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGTTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGATATACGGTAGTTTTTa                                                                                                                       >  2:605817/1‑150 (MQ=255)
                            gCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaa                                                                                                                      >  1:1443493/1‑150 (MQ=255)
                             cAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaaa                                                                                                                     >  2:381588/1‑150 (MQ=255)
                                  tCAGCGTGAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAatcat                                                                                                                >  2:789620/1‑150 (MQ=255)
                                    tgCGTATATATGCGTGGAGCCTTGAGAGGTCATGTGGGTCTTCGTGTGTAACAGCTAAATGCAGAAGTTCGGGGGTTTGATTTTGTTCATGACACTTCAGGCATTCCAGCAGTGCCTGCTGGTGATTTACGGTAGTTTTTATAATCATCt                                                                                                              <  2:1441102/149‑1 (MQ=255)
                                     gCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTg                                                                                                             >  2:1571972/1‑150 (MQ=255)
                                        tAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCAc                                                                                                          >  1:472535/1‑150 (MQ=255)
                                            gATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCg                                                                                                      >  1:948375/1‑150 (MQ=255)
                                             aTGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGt                                                                                                     >  2:741587/1‑150 (MQ=255)
                                                 gTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTc                                                                                                 >  1:964986/1‑150 (MQ=255)
                                                  ttGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCt                                                                                                <  2:374733/150‑1 (MQ=255)
                                                   tGAGCCATGAGAGTTAATCTGTGTCTCCTTGTGTAACAGGTAAAACCAGAAGTTCATGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCtt                                                                                               >  2:634962/1‑150 (MQ=255)
                                                   tGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCtt                                                                                               <  1:1436633/150‑1 (MQ=255)
                                                     aGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGCCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGc                                                                                             >  1:300278/1‑150 (MQ=255)
                                                      gCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCt                                                                                            >  2:1295913/1‑150 (MQ=255)
                                                      gCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCt                                                                                            >  2:1383716/1‑150 (MQ=255)
                                                      gCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCt                                                                                            <  2:775058/150‑1 (MQ=255)
                                                        cATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTcc                                                                                          <  1:634480/150‑1 (MQ=255)
                                                          tGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAatcatc                                                                                                               >  2:1457054/1‑127 (MQ=255)
                                                          tGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAatcatc                                                                                                               <  1:1457054/127‑1 (MQ=255)
                                                          tGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAg                                                                                        <  2:102412/150‑1 (MQ=255)
                                                             gagGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCtt                                                                                     >  1:1624842/1‑150 (MQ=255)
                                                             gagGTAATCTGTGTCTCAGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCtt                                                                                     <  1:726972/150‑1 (MQ=255)
                                                              agGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTa                                                                                    <  1:1479379/150‑1 (MQ=255)
                                                                 tAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCt                                                                                 >  1:1408982/1‑150 (MQ=255)
                                                                   aTCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTc                                                                               >  2:935908/1‑150 (MQ=255)
                                                                    tCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTcc                                                                              >  1:880898/1‑150 (MQ=255)
                                                                     cTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTccc                                                                             <  1:381588/150‑1 (MQ=255)
                                                                      tgtgTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCa                                                                            <  2:1209905/150‑1 (MQ=255)
                                                                      tgtgTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCa                                                                            <  2:183072/150‑1 (MQ=255)
                                                                       gtgtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAg                                                                           <  2:964986/150‑1 (MQ=255)
                                                                       gtgtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAg                                                                           >  2:573926/1‑150 (MQ=255)
                                                                          tctcCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTaa                                                                        >  2:328959/1‑150 (MQ=255)
                                                                             ccGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATag                                                                     >  1:483449/1‑150 (MQ=255)
                                                                             ccGTGTGTAACAGGTAAAACCACAAGTTACGGGGTTTGATATGCATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATag                                                                     >  2:533092/1‑150 (MQ=255)
                                                                                tgtgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagaga                                                                  >  2:347181/1‑150 (MQ=255)
                                                                                tgtgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATAGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAACagaga                                                                  >  2:1473532/1‑150 (MQ=255)
                                                                                 gtgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagagag                                                                 >  2:804049/1‑150 (MQ=255)
                                                                                   gtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagagagag                                                               >  2:1540889/1‑150 (MQ=255)
                                                                                    tAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAgg                                                              >  2:671781/1‑150 (MQ=255)
                                                                                    tAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCATGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAgg                                                              >  1:844289/1‑150 (MQ=255)
                                                                                      aCAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTa                                                            >  1:1150318/1‑150 (MQ=255)
                                                                                         ggTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                         >  1:1382551/1‑150 (MQ=255)
                                                                                             aaaCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCgg                                                     >  1:682076/1‑150 (MQ=255)
                                                                                              aaCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCggg                                                    >  2:1137738/1‑150 (MQ=255)
                                                                                               aCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCgggg                                                   >  2:1256195/1‑150 (MQ=255)
                                                                                                    aaGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCaa                                              <  1:279288/150‑1 (MQ=255)
                                                                                                     aGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAg                                             <  1:1571972/150‑1 (MQ=255)
                                                                                                         ccGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcg                                         >  1:246057/1‑150 (MQ=255)
                                                                                                           ggggTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGTCGGGGATCAAGGTcgcg                                       <  2:878104/150‑1 (MQ=255)
                                                                                                            gggTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCgcgc                                      <  2:300278/150‑1 (MQ=255)
                                                                                                             ggTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCa                                     <  1:671781/150‑1 (MQ=255)
                                                                                                              gTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAg                                    >  1:1053903/1‑150 (MQ=255)
                                                                                                              gTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAg                                    <  2:1443493/150‑1 (MQ=255)
                                                                                                                 tGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgcc                                 <  1:1435294/150‑1 (MQ=255)
                                                                                                                   atatGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccgc                               >  1:1454584/1‑150 (MQ=255)
                                                                                                                    tatGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccgcc                              >  1:970257/1‑150 (MQ=255)
                                                                                                                      tGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAg                            <  1:957474/150‑1 (MQ=255)
                                                                                                                      tGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAg                            >  1:1162710/1‑150 (MQ=255)
                                                                                                                      tGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAg                            <  2:569192/150‑1 (MQ=255)
                                                                                                                       ggATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGt                           >  2:319073/1‑150 (MQ=255)
                                                                                                                         atgatgAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGttt                         <  1:533092/150‑1 (MQ=255)
                                                                                                                         atgatgAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGGTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGttt                         >  1:251641/1‑150 (MQ=255)
                                                                                                                          tggtgaGGCTTCGGGCATTCCAGCGGTGCCTGCTGGTGGGATACGGTAGGTTTTAGAGTCGTGTGCACTTGGTATTCTTGCGGCAGCTTGCCTGCGGAGGAATAGAGGGAGGGAGGGCGGGGGGTCAGGGTCGCGCAGGCCGCCAGtttt                        >  2:1582880/4‑150 (MQ=255)
                                                                                                                          tgatgaGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGtttt                        <  2:966095/150‑1 (MQ=255)
                                                                                                                            atgaGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTc                      <  1:1272926/150‑1 (MQ=255)
                                                                                                                               agaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGc                   >  2:200857/1‑150 (MQ=255)
                                                                                                                               agaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGc                   >  2:1588519/1‑150 (MQ=255)
                                                                                                                                gaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                  <  1:573926/150‑1 (MQ=255)
                                                                                                                                gaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                  >  1:240886/1‑150 (MQ=255)
                                                                                                                                gaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                  <  1:130473/150‑1 (MQ=255)
                                                                                                                                  cTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAg                >  1:1509752/1‑150 (MQ=255)
                                                                                                                                  cTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACATTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAg                >  1:202961/1‑150 (MQ=255)
                                                                                                                                     cAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCAc             >  1:558867/1‑150 (MQ=255)
                                                                                                                                     cAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCAc             >  2:955687/1‑150 (MQ=255)
                                                                                                                                        gCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACttt          <  2:1287733/150‑1 (MQ=255)
                                                                                                                                         ctcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCCCTTCGTATTCTTGCTCCAGCTTACCTTCCCAGGAATAGAGCGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTGCCGCCCGCACGTTg         >  2:1632412/5‑150 (MQ=255)
                                                                                                                                         ctcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTg         >  2:66812/5‑150 (MQ=255)
                                                                                                                                         cATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAACACTTTg         >  2:49060/1‑150 (MQ=255)
                                                                                                                                          tcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg        >  2:99390/4‑150 (MQ=255)
                                                                                                                                          tcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCAAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg        >  2:670183/4‑150 (MQ=255)
                                                                                                                                          aTTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg        >  1:1470867/1‑150 (MQ=255)
                                                                                                                                            accAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcg      >  2:911328/2‑150 (MQ=255)
                                                                                                                                                gcagTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGGCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTa  <  2:1463279/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTA  >  NC_000913/4364984‑4365277

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