Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,082,500 T→A 7.2% F86L (TTT→TTA) ‡ efeO → inactive ferrous ion transporter EfeUOB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,082,5000TA 7.2% 535.7 / 22.8 151F86L (TTT→TTA) ‡efeOinactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base T (70/70);  new base A (6/5);  total (76/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 157 total aligned reads displayed.

CGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTAAGTCTTGGTGGTGC  >  NC_000913/1082352‑1082646
                                                                                                                                                    |                                                                                                                                                  
cGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCttt                                                                                                                                                    <  4:356618/149‑1 (MQ=255)
  gATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAg                                                                                                                                                  >  1:684529/1‑149 (MQ=255)
   aTAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                                                                                                                 <  2:684529/149‑1 (MQ=255)
    tAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGcc                                                                                                                                                <  2:434844/149‑1 (MQ=255)
        cAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaa                                                                                                                                            <  2:97341/149‑1 (MQ=255)
          gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTAAAGCCAggg                                                                                                                                            <  4:453768/147‑3 (MQ=255)
          gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTAAAGCCAggg                                                                                                                                            >  3:453768/1‑145 (MQ=255)
              gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                                                                                                                 <  3:350102/138‑1 (MQ=255)
              gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                                                                                                                 >  4:350102/1‑138 (MQ=255)
               aaCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                                                                     <  2:250648/149‑1 (MQ=255)
                        ccATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTa                                                                                                                            >  2:972345/1‑149 (MQ=255)
                           ttACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAg                                                                                                                         >  4:238621/1‑149 (MQ=255)
                           ttACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAg                                                                                                                         >  1:999980/1‑149 (MQ=255)
                               ggTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTAAAGCCAggg                                                                                                                                            >  4:800965/1‑124 (MQ=255)
                               ggTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTAAAGCCAggg                                                                                                                                            <  3:800965/126‑3 (MQ=255)
                                 ttAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                                                                                                                   >  3:630640/1‑149 (MQ=255)
                                         gggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           <  1:1262775/149‑1 (MQ=255)
                                         gggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           >  1:532093/1‑149 (MQ=255)
                                          ggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTa                                                                                                                                                   >  4:858776/1‑108 (MQ=255)
                                          ggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                                                                                                                   >  3:1083680/1‑140 (MQ=255)
                                           gAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGata                                                                                                         <  1:972345/149‑1 (MQ=255)
                                               acacAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGAc                                                                                                     >  4:1158850/1‑149 (MQ=255)
                                                 acaGTTCATTACTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGACGTGATGGTGGTGGAAGAGCGGGAAAATATCTCCCCTGGCTATAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCt                                                                                                   >  4:1004542/1‑149 (MQ=255)
                                                   aGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGc                                                                                                 >  2:611924/1‑149 (MQ=255)
                                                     ttCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           >  4:503604/1‑137 (MQ=255)
                                                     ttCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           <  3:503604/137‑1 (MQ=255)
                                                       cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTAAAGCCAggg                                                                                                                                            >  1:1402048/1‑100 (MQ=255)
                                                        attattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAAt                                                                                                                                         >  4:625498/1‑104 (MQ=255)
                                                        attattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                                                                                                      >  2:1098048/1‑107 (MQ=255)
                                                        attattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                                                                                                      <  1:1098048/107‑1 (MQ=255)
                                                        attattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                                                            >  1:315564/1‑149 (MQ=255)
                                                        attataCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAAt                                                                                                                                         <  3:625498/104‑1 (MQ=255)
                                                           attCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                                                                         >  1:1340137/1‑149 (MQ=255)
                                                           attCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                                                                         <  1:207289/149‑1 (MQ=255)
                                                            ttCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGata                                                                                                         >  2:1236334/1‑132 (MQ=255)
                                                            ttCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGata                                                                                                         <  1:1236334/132‑1 (MQ=255)
                                                            ttCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgctg                                                                                        >  3:946885/1‑149 (MQ=255)
                                                             tCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGa                                                                                       >  2:1200245/1‑149 (MQ=255)
                                                              cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGAc                                                                                      <  1:611924/149‑1 (MQ=255)
                                                               aGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTAAAGCCAggg                                                                                                                                            <  3:567566/94‑3 (MQ=255)
                                                                gAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTa                                                                                    >  4:1052036/1‑149 (MQ=255)
                                                                    cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGa                                                                                                                                       <  4:1133733/94‑1 (MQ=255)
                                                                    cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAccc                                                                                <  1:860495/149‑1 (MQ=255)
                                                                       aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                                                                     >  2:568733/1‑93 (MQ=255)
                                                                        gCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                                                            <  4:630640/149‑1 (MQ=255)
                                                                            gAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGa                                                                        <  2:329021/149‑1 (MQ=255)
                                                                             aaGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGaa                                                                       <  1:1200245/149‑1 (MQ=255)
                                                                             aaGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGaa                                                                       <  3:1158850/149‑1 (MQ=255)
                                                                             aaGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGaa                                                                       <  2:1340137/149‑1 (MQ=255)
                                                                                gcgcTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGtt                                                                    >  1:96890/1‑149 (MQ=255)
                                                                                 cgcTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTg                                                                   <  4:93194/149‑1 (MQ=255)
                                                                                      gAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           >  3:1401309/1‑104 (MQ=255)
                                                                                      gAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           <  4:1401309/104‑1 (MQ=255)
                                                                                        gTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATAc                                                                                                             >  4:485144/1‑100 (MQ=255)
                                                                                        gTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATAc                                                                                                             <  3:485144/100‑1 (MQ=255)
                                                                                           ggAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGt                                                                                              >  1:150565/1‑112 (MQ=255)
                                                                                           ggAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                                                            >  1:331305/1‑114 (MQ=255)
                                                                                           ggAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                                                            <  2:331305/114‑1 (MQ=255)
                                                                                           ggAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAg                                                         >  1:585705/1‑149 (MQ=255)
                                                                                               aTCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           <  4:58161/95‑1 (MQ=255)
                                                                                               aTCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                                                           >  3:58161/1‑95 (MQ=255)
                                                                                               aTCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa                                                     >  3:899618/1‑149 (MQ=255)
                                                                                                 ccTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTgg                                                                                                                    >  3:1150119/1‑84 (MQ=255)
                                                                                                 ccTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAgg                                                   <  3:783496/149‑1 (MQ=255)
                                                                                                  cTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGTCTGGCGAATACGATATGa                                                                                                      >  3:1095320/1‑97 (MQ=255)
                                                                                                  cTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGTCTGGCGAATACGATATGa                                                                                                      <  4:1095320/97‑1 (MQ=255)
                                                                                                   tCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                                                            <  3:1434916/122‑1 (MQ=255)
                                                                                                   tCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCa                                                 >  3:1023818/1‑149 (MQ=255)
                                                                                                    cAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAg                                                    >  2:1294236/1‑145 (MQ=255)
                                                                                                    cAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAg                                                    <  1:1294236/145‑1 (MQ=255)
                                                                                                           gTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGt                                                                     <  1:682859/121‑1 (MQ=255)
                                                                                                           gTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGa                                         >  4:692201/1‑149 (MQ=255)
                                                                                                             gATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAccc                                                                                <  2:66268/108‑1 (MQ=255)
                                                                                                                      ggAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa                                                          <  3:23204/121‑1 (MQ=255)
                                                                                                                      ggAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa                                                          >  4:23204/1‑121 (MQ=255)
                                                                                                                       gAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTg                                                                                                                     >  1:1036437/1‑61 (MQ=255)
                                                                                                                       gAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTg                                                                                                                     <  2:1036437/61‑1 (MQ=255)
                                                                                                                       gAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAgg                                                   >  4:553441/1‑127 (MQ=255)
                                                                                                                              gggAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAttt                                                                                                                             <  4:1191104/46‑1 (MQ=255)
                                                                                                                              gggAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAttt                                                                                                                             >  3:1191104/1‑46 (MQ=255)
                                                                                                                              gggAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa                                                     >  2:1347205/1‑118 (MQ=255)
                                                                                                                              gggAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATgcg                      >  2:127261/1‑149 (MQ=255)
                                                                                                                               ggAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATgcgc                     <  4:142818/149‑1 (MQ=255)
                                                                                                                               ggAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATgcgc                     <  3:280136/149‑1 (MQ=255)
                                                                                                                                gAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa                                                          >  2:852430/1‑111 (MQ=255)
                                                                                                                                   aaTATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                                                                         <  1:1406755/77‑1 (MQ=255)
                                                                                                                                   aaTATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                                                                         >  2:1406755/1‑77 (MQ=255)
                                                                                                                                     tatCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTaa               >  2:201128/1‑149 (MQ=255)
                                                                                                                                      atCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGt                                                       <  2:145691/108‑1 (MQ=255)
                                                                                                                                      atCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGt                                                       >  1:145691/1‑108 (MQ=255)
                                                                                                                                        cGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATc                                                                >  1:486513/1‑97 (MQ=255)
                                                                                                                                        cGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATc                                                                <  2:486513/97‑1 (MQ=255)
                                                                                                                                          ccccTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTAAGTCtt          >  4:770699/1‑149 (MQ=255)
                                                                                                                                             cTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAAcgg                                             <  4:1373800/111‑1 (MQ=255)
                                                                                                                                             cTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAAcgg                                             >  3:1373800/1‑111 (MQ=255)
                                                                                                                                             cTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTAAGTCTtggt       >  4:371833/1‑149 (MQ=255)
                                                                                                                                              tGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa                                                     <  1:1039359/102‑1 (MQ=255)
                                                                                                                                              tGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTAAGTCTtggtg      >  2:259008/1‑149 (MQ=255)
                                                                                                                                                  tttAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTAAGTCTTGGTGGTGc  >  4:1026244/1‑149 (MQ=255)
                                                                                                                                                  tttAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTAAGTCTTGGTGGTGc  <  1:127261/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGGCGGATGCGGCGCAAAGTGATGCGCTGTTAAGTCTTGGTGGTGC  >  NC_000913/1082352‑1082646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: