Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,467,145 A→G L649L (TTA→TTG mgtA → Mg2+ transport ATPase, P‑type 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,467,1450AG99.5% 1267.8 / ‑3.6 370L649L (TTA→TTGmgtAMg2+ transport ATPase, P‑type 1
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (168/199);  total (170/200)
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.01e-01

Only 100 of 400 total aligned reads displayed.

GATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCTGAT  >  NC_000913/4467055‑4467234
                                                                                           |                                                                                         
gATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGcg                                                                                         >  1:4494703/1‑93 (MQ=255)
 agccATAAAGAACGCATTTGTGCCTTACTGAAGCACGAACGGCATGTGGTTGGTTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTgcgc                                                                                        <  2:2087279/90‑1 (MQ=255)
 aTGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTgcgc                                                                                        <  1:2201784/93‑1 (MQ=255)
 aTGCATAAAGAACGCATTATAGCCTTACTGAAGCGCGAAGGACATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTgcgc                                                                                        <  2:3639801/93‑1 (MQ=255)
  tGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGcgcg                                                                                       >  2:8719390/1‑93 (MQ=255)
    cATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCCTGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCt                                                                                     <  1:8395008/93‑1 (MQ=255)
    cATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCt                                                                                     >  1:2542123/1‑93 (MQ=255)
    cATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCt                                                                                     >  1:5809019/1‑93 (MQ=255)
      tAAAGAACGCATTGTGACCTTACTGAAGTGCGAAGGGCTTGTGGTTGGCTTTATGGGCGATGGTATTAATGTTGTGCCCGCCTTTCGCGCTGc                                                                                   >  2:8040862/1‑93 (MQ=255)
      tAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGc                                                                                   >  1:659997/1‑93 (MQ=255)
       aaaGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCg                                                                                  <  2:5813675/93‑1 (MQ=255)
        aaGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCgg                                                                                 <  2:628611/93‑1 (MQ=255)
         aGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGa                                                                                >  2:680777/1‑93 (MQ=255)
         aGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGa                                                                                >  1:372155/1‑93 (MQ=255)
             cGCATTGTGACCTTACTGAAGCGCGAAGGGTATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATc                                                                            <  1:1657455/93‑1 (MQ=255)
              gCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCg                                                                           >  2:6622626/1‑93 (MQ=255)
              gCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCg                                                                           >  1:4843632/1‑93 (MQ=255)
              gCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCg                                                                           >  1:4293081/1‑93 (MQ=255)
              gCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCGCCGGatatag                                                                           >  1:8042439/1‑91 (MQ=255)
               cATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCgg                                                                          >  1:874477/1‑93 (MQ=255)
                aTTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGc                                                                         <  2:7760092/93‑1 (MQ=255)
                        ccTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCtg                                                                  >  1:8464065/1‑93 (MQ=255)
                        ccTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCtg                                                                  >  1:7521098/1‑93 (MQ=255)
                         cTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCtgt                                                                 <  2:1694244/93‑1 (MQ=255)
                         cTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCtgt                                                                 >  1:685896/1‑93 (MQ=255)
                          ttACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCtgtg                                                                <  2:346514/93‑1 (MQ=255)
                          ttACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCtgtg                                                                >  2:1243100/1‑93 (MQ=255)
                          ttACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCtgtg                                                                >  1:574643/1‑93 (MQ=255)
                           tACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTgg                                                               <  2:3935682/93‑1 (MQ=255)
                            aCTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGa                                                              >  1:5350666/1‑93 (MQ=255)
                            aCTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGa                                                              <  2:2614546/93‑1 (MQ=255)
                             cTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGAc                                                             <  1:591390/93‑1 (MQ=255)
                             cTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGAc                                                             <  1:779924/93‑1 (MQ=255)
                             cTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGAc                                                             <  1:5544926/93‑1 (MQ=255)
                              tGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACg                                                            >  1:4132257/1‑93 (MQ=255)
                                aaGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGgc                                                          >  2:1575522/1‑93 (MQ=255)
                                   cgcgAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGcgcg                                                       <  2:2398669/93‑1 (MQ=255)
                                   cgcgAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGcgcg                                                       >  2:5354804/1‑93 (MQ=255)
                                    gcgAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCgg                                                      <  1:2083147/93‑1 (MQ=255)
                                    gcgAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCgg                                                      >  1:2862129/1‑93 (MQ=255)
                                    gcgAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGACGGCGCggg                                                     >  1:8156188/1‑92 (MQ=255)
                                     cgAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGt                                                     >  1:5321835/1‑93 (MQ=255)
                                      gAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTa                                                    <  2:4519714/93‑1 (MQ=255)
                                       aaGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAg                                                   >  1:671181/1‑93 (MQ=255)
                                       aaGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAg                                                   >  1:7624208/1‑93 (MQ=255)
                                       aaGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAg                                                   >  2:7274592/1‑93 (MQ=255)
                                       aaGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTGGCGCGCTGCGGATATCGGCATTTCTGTGGACGCCGCGGTAg                                                   >  1:5152916/1‑93 (MQ=255)
                                        aGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGa                                                  >  2:7913243/1‑93 (MQ=255)
                                         gggCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGat                                                 >  2:6758095/1‑93 (MQ=255)
                                          ggCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGata                                                <  1:3095874/93‑1 (MQ=255)
                                          ggCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGata                                                >  2:8597299/1‑93 (MQ=255)
                                          ggCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGata                                                >  1:2716602/1‑93 (MQ=255)
                                          ggCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGata                                                >  1:3851243/1‑93 (MQ=255)
                                           gCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGatat                                               <  2:1899223/93‑1 (MQ=255)
                                           gCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGatat                                               <  2:2885284/93‑1 (MQ=255)
                                            gaTGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATAtt                                              <  1:6083670/92‑1 (MQ=255)
                                            cATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATAtt                                              <  1:4921193/93‑1 (MQ=255)
                                            cATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATAtt                                              <  2:3732905/93‑1 (MQ=255)
                                               gtgGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGcc                                           >  1:1285039/1‑93 (MQ=255)
                                                 ggTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCg                                         >  2:3084150/1‑93 (MQ=255)
                                                  gTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGt                                        >  1:966926/1‑93 (MQ=255)
                                                   ttGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTg                                       >  1:7430689/1‑93 (MQ=255)
                                                     ggCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGaa                                     <  2:4615564/93‑1 (MQ=255)
                                                      gCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAg                                    >  2:8427050/1‑93 (MQ=255)
                                                      gCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAg                                    >  2:2328341/1‑93 (MQ=255)
                                                       cTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGc                                   <  1:2032365/93‑1 (MQ=255)
                                                       cTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGc                                   <  1:1390843/93‑1 (MQ=255)
                                                       cTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGc                                   <  1:1027878/93‑1 (MQ=255)
                                                          tATGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGc                                <  2:7483195/93‑1 (MQ=255)
                                                           aTGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCt                               <  1:6464317/93‑1 (MQ=255)
                                                           aTGGGCGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCt                               <  2:4979606/93‑1 (MQ=255)
                                                                cGATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGatat                          >  2:5896098/1‑93 (MQ=255)
                                                                 gATGGTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATatc                         <  1:4680380/93‑1 (MQ=255)
                                                                  aTGGTATTAATAATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATatca                        <  1:1925648/93‑1 (MQ=255)
                                                                     gTATTAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATcc                     <  1:7127548/93‑1 (MQ=255)
                                                                     gTATTAATGATGAGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATatcattc                     >  2:4134204/1‑91 (MQ=255)
                                                                        ttAATGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCctgc                  <  2:7663987/93‑1 (MQ=255)
                                                                          aaTGATGCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCctgctg                >  2:5144997/1‑93 (MQ=255)
                                                                           atgatgCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCCTCCTGCTgg               >  2:1800048/1‑93 (MQ=255)
                                                                             gatgCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGaa             <  2:4673689/93‑1 (MQ=255)
                                                                             gatgCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGaa             >  1:7931746/1‑93 (MQ=255)
                                                                               tgCGCCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGaaaa           >  2:1748071/1‑93 (MQ=255)
                                                                                gcgcCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGaaaaa          <  2:6971188/93‑1 (MQ=255)
                                                                                gcgcCCGCCTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCATGCTGGgaaaa          >  2:6912866/1‑93 (MQ=255)
                                                                                  gcccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAg        <  2:3006482/93‑1 (MQ=255)
                                                                                   cccgccTTGCGCGCTGCTGATATCGGAATTACTGTGGTCGGCGCGGTAGTTATTGCCCGTGAAGCGGCTGATATCATCCTTCTTGAAAAAAGc       <  2:3314614/93‑1 (MQ=255)
                                                                                   cccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGc       <  1:3019808/93‑1 (MQ=255)
                                                                                    ccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGcc      <  2:5550076/93‑1 (MQ=255)
                                                                                    ccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGcc      <  1:3181908/93‑1 (MQ=255)
                                                                                    ccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGcc      <  2:1443198/93‑1 (MQ=255)
                                                                                    ccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGcc      <  1:239355/93‑1 (MQ=255)
                                                                                    ccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGcc      <  1:4417554/93‑1 (MQ=255)
                                                                                    ccgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGGAGATATTGCCCGTGAAGCGGCTGATATCATCGTGCTGGAAAAGAGcc      >  2:6433707/1‑93 (MQ=255)
                                                                                     cgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCt     >  1:5943712/1‑93 (MQ=255)
                                                                                     cgccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATAATCCTGCTGGAAAAAAGCCt     >  2:5661595/1‑93 (MQ=255)
                                                                                      gccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCTg    <  1:7041807/93‑1 (MQ=255)
                                                                                      gccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCTg    >  2:8224732/1‑93 (MQ=255)
                                                                                       ccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCTGa   <  1:1496479/93‑1 (MQ=255)
                                                                                       ccTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCTGa   >  1:346366/1‑93 (MQ=255)
                                                                                        cTTGCGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCTGAt  <  2:5562278/93‑1 (MQ=255)
                                                                                           |                                                                                         
GATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGGCGCGGTAGATATTGCCCGTGAAGCGGCTGATATCATCCTGCTGGAAAAAAGCCTGAT  >  NC_000913/4467055‑4467234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: